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Abstract of the Dissertation 

Dilepton spectra in p + p and Au + Au 
collisions at RHIC 

by 
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Doctor of Philosophy 

in 

Physics 

Stony Brook University 
2008 



Recent experimental results show that a strongly coupled quark- 
gluon plasma (sQGP) is created in heavy ion collisions at the Rela- 
tivistic Heavy Ion Collider (RHIC) at Brookhaven National Labo- 
ratory (BNL). Electromagnetic radiation, i. e., photons and lepton 
pairs, are penetrating probes that allow investigating the full time 
evolution and dynamics of the matter produced, as they do not 
undergo strong interaction in the final state. 

This work presents measurements of electron-positron pairs from 
p + p collisions at y/s = 200 GeV collected during the 2005 RHIC 
run and compares them to results from Au + Au collisions at ^s^n 
= 200 GeV taken in 2004 with the PHENIX detector. The invari- 
ant mass distribution of e + e~ pairs in p + p is consistent with the 
expected contributions from Dalitz decays of light hadrons, dielec- 
tron decays of vector mesons and correlated charm production, 
which have been measured in the same experiment. The charm 
and bottom cross section extracted from the measured dielectron 
yield are = 544 ± 39(stat.) ± 142(syst.) ± 200 (model) /xb and 
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a b i = 3.9±2.4(stat.)l 2 (syst.) fib, respectively. The dielectron con- 
tinuum measurement in p + p provides a crucial baseline for the 
modification of the dielectron continuum observed in Au + Au. 

In min. bias Au + Au collisions the yield of dielectrons in the 
low mass region (150 < m ee < 750 MeV/c 2 ) is enhanced by a 
factor of 4.0 ± 0.3(stat.) ± 1.5(syst.) ± 0.8(model) compared to 
the known hadronic sources. The centrality dependence of this en- 
hancement suggests emission from in-medium scattering processes. 
The excess dominates the yield in the transverse momentum re- 
gion below 1 GeV/c and shows significantly lower (p?) than the 
expected sources. The low p? enhancement is currently not un- 
derstood by any theoretical model of heavy ion collisions. The en- 
hancement extends to larger transverse momenta (pr > 1 GeV/c) 
where it is also observed in p + p and explained by virtual direct 
photons. The p + p measurement serves as an important test to 
pQCD calculations of direct photon production from hard scat- 
tering processes in this momentum range. An excess with an in- 
verse slope of T eS = 221 ± 23(stat.) ± 18(syst.) MeV is observed 
in central Au + Au collisions above the binary scaled direct pho- 
ton yield in p + p. This can be qualitatively explained by hydro- 
dynamical models including thermal photon radiation with initial 
temperatures of 300 < T init < 600 MeV and formation times of 
0.12 < t < 0.6 fm/c. 
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Chapter 1 
Introduction 



Relativistic heavy ion collisions, i. e., collisions in which the projectile energies 
are much larger than their rest masses, have been an excellent tool to study 
nuclear matter under extreme conditions and explore the phase diagram of 
strongly interacting matter for new states of matter. At the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) gold ions are 
brought to collisions at center of mass energies of up to y/sJ^N = 200 GeV. 
Before RHIC, the Alternating Gradient Synchrotron (AGS) at BNL and the 
Super Proton Synchrotron (SPS) have been colliding heavy ions at energies 
up to v /J/viv ~ 5 GeV (AGS) and y/s NN « 17 GeV and therefore studying 
nuclear matter at lower temperatures and higher net baryon densities. With 
even lower energies the heavy ion synchrotron "Schwerlonen Synchrotron" 
(SIS) provides heavy ion collisions at beam energies up to 2 AGeV. The Large 
Hadron Collider (LHC) will achieve even higher temperatures in lead on lead 
collisions at ^s^n = 5.5 TeV. 

At energy densities of e ~ 1 GeV/fm 3 a phase transition from hadronic 
matter to deconfmed quarks and gluons, the so called quark-gluon plasma 
(QGP), is predicted [1]. Calculations of lattice quantum chromodynamics 
(1QCD) predict that such an energy density is reached at a temperature of 
T m 170 MeV 1 ~ 10 12 K [2]. This state may have existed in the early universe 
a few tens of microseconds after the big-bang and may still exist in the core 
of neutron stars at densities of 10 18 kg/m 3 = 0.6 fm~ 3 which exceeds normal 
nuclear matter density by more than a factor of four. The same calculations 
also show that, at vanishing net baryon density, with the transition into the 
deconfmed phase, chiral symmetry, which is spontaneously broken in vacuum 
due to non- vanishing effective quark masses, is restored [2]. 



1 Natural units fi, = c = fcs = l are employed throughout this thesis except where noted 
otherwise 
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Figure 1.1: Schematic phase diagram of QCD matter as function of temperature 
T and baryonic chemical potential /U#. The measured chemical freeze out points 
for SIS, AGS, SPS, and RHIC energies are shown as points [5]. The dashed line 
denotes the thermal freeze out. The existence and exact location of a critical point 
are unknown. 



In Fig. 1.1 a schematic phase diagram of QCD matter [3-5] is shown as 
function of temperature T and baryon chemical potential /xg. Below the phase 
boundary quarks and gluons are confined to hadrons which form a interacting 
hadron gas. Above the phase transition, matter of asymptotically free quarks 
and gluons with partonic degrees of freedom is created. The exact details such 
as the order of the phase transition are unknown. Lattice QCD calculations 
indicate a cross over at vanishing net baryon density, while at large net baryon 
densities a first order phase transition is expected. If such a scenario is real- 
ized, a critical point, below which the transition becomes a cross over, would 
terminate the first order phase transition. 

Also shown are for different collision energies, the points at which the 
fireball created in the collision, reaches its chemical freeze out, i. e. the end of 
inelastic collisions and therefore the production of new particles has stopped. 
The thermal freeze out refers to the end of elastic collision, i. e. the end of 
momentum transfer between particles, as the mean free path of the particles 
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has become larger than the size of the system. The experimental points of 
chemical freeze out are calculated based on the final multiplicities of particle 
species [6, 7]. 

Experimental results from all four RHIC experiments (BRAHMS, PHENIX, 
PHOBOS, and STAR) give clear evidence that a new state of matter has been 
created in Au + Au collisions at y/sNN = 200 GeV [8-11]. These results include 
that a very high energy density (e) ~ 15 GeV/fm 3 is achieved, as indicated by 
the large energy loss of light hadrons [12, 13], and heavy quarks [14-16]. Fur- 
thermore, a large elliptic flow is developed on a partonic level including heavy 
quarks which is evidence for an early thermalization of the medium[14, 17-19]. 

Colorless, electromagnetic probes such as photons and dileptons (electron- 
positron or muon pairs) are particular useful probes to study the the quark- 
gluon plasma, as they, once emitted, do not undergo strong final state interac- 
tions with the hadronic medium. Dileptons are created via various processes 
during all stages of the collision. There is an approximate time ordering in 
the invariant mass of the lepton pair; the earlier produced the larger its mass. 
In the earliest stage, Drell-Yan annihilation between incoming qq pairs cre- 
ates dileptons dominating the invariant mass region well above 3 GeV/c 2 . As 
the formed hot and dense medium rapidly thermalizes, a significant contri- 
bution of dileptons from qq annihilation is expected, but in contrast to the 
DY pairs, these pairs should follow a thermal distribution in the intermediate 
mass region (IMR) 1 < ran < 3 GeV/c 2 . The major background source for 
the thermal radiation in the IMR comes from semi-leptonic decays of open 
charm. Pairs of cc quarks created in the initial hard scattering independently 
hadronize to DD mesons, which consist of one charm and one light quark, 
therefore often referred to as "open charm" mesons. These mesons inherit the 
strong initial back-to-back correlation of the charm quarks. Their individual 
weak decay D — > Klv\ leads to a continuum of lepton pairs dominating the 
IMR. Medium modifications of charm quarks, which leads to a significant flow 
and suppression at large transverse momenta of single electrons from open 
charm decays [15], may also be reflected in the invariant mass distribution of 
dileptons from open charm decays. After hadronization the main contribution 
is expected from annihilation of pions and other hadrons. Of particular in- 
terest are the two body annihilation of pions through the light vector mesons 
p, U, and <p, as they decay directly into lepton pairs, whose invariant mass 
reflects the mass of the vector meson at the time of the decay. With the short 
life time of the p of r p = 1.3 fm and a life time of the hadronic medium of 
pa 10 fm, this channel provides an excellent probe to the effects of chiral sym- 
metry restoration [20]. Once the chemical freeze out is reached, the hadronic 
resonances and Dalitz decays of light pseudo-scalar mesons as 7r°, rj contribute 
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to the low mass region (LMR) below 1 GeV/c 2 . 

A schematic view of all contributions is shown in Fig. 1.2 and is supple- 
mented with a realistic "cocktail" of hadron decays as expected to be measured 
with PHENIX in p + p collisions shown in Fig. 1.3. This cocktail includes the 
effects of detector resolution and acceptance. 

This thesis is about the analysis of the dielectron continuum in the low and 
intermediate mass region in p + p collisions at y/s = 200 GeV recorded in 2005 
with the PHENIX detector and the comparison to the dielectron continuum 
measured in Au + Au collisions taken in 2004. As part of this effort the 
cross sections of the u and in p + p collisions mesons have been measured. 
The e + e~ pairs in the intermediate mass region allowed the extraction of the 
total charm cross section in p + p [21], which is in very good agreement with 
previous measurements via single electrons [16]. The high mass region (HMR) 
m ee > 3 GeV gave access to the first measurement of the bottom cross section 
in p + p collisions at \fs = 200 GeV. 

In the continuum of the low mass region a contribution from internal con- 
versions of direct virtual photons has been extracted and provided an im- 
portant reference for the measurement of direct virtual photons in Au + Au 
collisions [22] . Furthermore, the measurement of the full dielectron continuum 
over a wide range in mass (0 < m ee < 8 GeV/c 2 ) and px (0 < pt < 5GeV/c) 
in p + p provides an important baseline for the interpretation of the Au + Au 
result [23]. 

This thesis is structured in the following way: In Chapter 1 an introduction 
into the physics of the dielectron continuum is given. After a brief discussion of 
the theoretical background in Section 1.1 and the different dilepton sources 1.2 
an overview over recent results at other dilepton experiments is given in Sec- 
tion 1.3. In Chapter 2 a description of the PHENIX detector is given. The 
analysis details are presented in Chapter 3. The various results of this analysis 
are shown and discussed in Chapter 4. The thesis ends with a summary and 
an outlook for further analyses of the dielectron continuum in Chapter 5. 
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Figure 1.2: Expected sources of dilepton production as function of invariant mass 
in ultra-relativistic heavy ion collisions. 
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Figure 1.3: Expected sources of dilepton production as function of invariant mass 
in p + p collisions at y/s = 200 GeV. 
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Figure 1.4: Examples of Feynman diagrams for the interaction between two 
charged particles, and the Compton scattering of a photon on an electron. Time is 
advancing from the left to the right. 

1.1 Quantum Chromodynamics 

The theory of quantum electrodynamics (QED) is probably the best under- 
stood and tested theory in the field of physics. On very short time scales, 
which can be estimated by the uncertainty principle (AEAt m H), electrons 
are allowed to emit and absorb a virtual photon with non-zero rest mass. The 
electromagnetic interaction between two charged particles is described by the 
exchange of such virtual photons. As the photon is charge neutral it does not 
interact with itself. The strength of the electromagnetic interaction is given 
by the coupling constant g e = yAna which depends on the fine structure 
constant a = e 2 / [hdireo) ~ 1/137. QED predicts a logarithmic increase of 
the interaction strength with increasing momentum transfers. The interaction 
between two charged particles, annihilation or scattering, can be represented 
by the Feynman diagram shown in Fig. 1.4. Also shown are the two diagrams 
for the Compton scattering of a photon from an electron. 

Higher order processes are responsible for changes in the strength of the 
electromagnetic interaction. As shown in Fig. 1.5 the virtual photon exchanged 
between two charges can create an e + e~ pair which subsequently annihilate 
each other. This virtual e + e~ pair acts as an electric dipole, effectively screen- 
ing the two interacting charges. Such vacuum polarization at short distances 
(large momenta) reduces the effective charge of the electron as well as the fine 
structure constant: 



\q 2 \ is the square of the momentum transfer, m the mass of the charged particle 
and Eq. (1.1) an approximation for large \q 2 \/m. Other higher order processes 
are responsible for the electron's self energy and its anomalous magnetic mo- 
ment (g s — 2) 7^ 1 and are shown in Fig. 1.5. 




(1.1) 
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Figure 1.6: Exemplary Feynman diagram for the strong interaction between quarks 

In an analog way, in the Standard Model quantum chromodynamics (QCD) 
describes the strong interaction between quarks via the exchange of color 
charges 2 carried by gluons. Fig. 1.6 shows the the lowest order Feynman 
diagrams for quark-antiquark annihilation and quark-gluon scattering. 

As gluons themselves carry a non-zero color charge, they can interact with 
themselves. This allows gluons not only to split into virtual quark-antiquark 
pairs, but also to split in pairs of gluon as shown in Fig. 1.7. This leads to 
an important difference from the electromagnetic interaction. The coupling 
strength of the strong interaction a s increases with increasing distance of two 
quarks. 

The momentum transfer dependence of the coupling strength can be writ- 
ten in analogy to Eq. (1.1) as: 

«s(|g 2 |)= f-^r + -^(11^-2^) In !*T) (1.2) 
\a s {/j, 2 ) 127T yu 2 J 

for {\q 2 \ ^> fi 2 ) with the number of colors N c = 3 and the number of flavors 



2 The necessity of color charges arises experimentally from the existence of baryons with 
three quarks of identical flavor, e. g. the A ++ consists of three u quarks. The Pauli exclusion 
principle demands an extra quantum number to allow such a configuration, as any two 
fcrmions must not occupy the same state but must at least differ in one quantum number. 
The introduction of color charges "red" , "blue" and "green" lifts the degeneracy of the three 
quarks. 
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Figure 1.7: Higher order processes of the strong interaction 

Nf = 6 (in the Standard Model). As for large distances or small momentum 
transfers the the coupling increases, known as anti-screening, a s cannot be 
expanded around \q 2 \ = 0, quite in contrast to a. With increasing momentum 
transfer the coupling decreases leading to quasi-free quarks and gluons, known 
as asymptotic freedom [24-27]. This allows an expansion of a s around large 
momentum transfers and therefore a perturbative treatment as long as a s < 1. 
However, for small momentum transfers on the order of q ~ 1 GeV/c the 
perturbative treatment breaks down and both the effective coupling and the 
relevant degrees of freedom change rapidly with scale — at large distances the 
degrees of freedom are colorless objects of two (mesons) or three (baryons) 
confined quarks. This regime of "strong QCD" poses the biggest challenge to 
the theory of the strong interaction. 

The QCD Lagrangian which describes the dynamics of the strong interac- 
tion can be expressed as: 

£qcd = ?(^D M - M^G^G""* (1.3) 

where tp is a vector of Dirac spinors 3 representing the wave function of a spin- 
1/2 quark-fields. 7^ are the Dirac matrices and M = diag(m u , m^, m s , . . . , m t ) 
is composed of the bare quark masses m u ~ 3, ~ 7, m s ~ 100 MeV/c 2 and 
m c ~ 1.25, nib ~ 4.1, m t ~ 175 GeV/c 2 . D M = + ig^-A 1 ^ is the covariant 
derivative with the Gell-Mann matrices A* and the spin-1 gauge fields 
with color index i = 1,...8. This term describes the interaction of quarks 
with gluons, to which the eight gauge fields correspond to. The field strength 
tensor 

G%, = d,K - <9,A; + igf^At (1.4) 



3 Each quark can carry one of three colors, red, blue or green. Therefore the vector of 
Dirac spinors reads: ip 
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describes the interaction of gluons with other gluons, where g = y/Aira s is 
the coupling constant of the strong interaction and fijk's denote the structure 
constants of the SU(3) group [28]. 

1.1.1 Chiral Symmetry 

The Lagrangian is invariant under local SU(3) gauge transformations, i. e., 
it is invariant under arbitrary rotations in color space. In addition, the La- 
grangian exhibits a global symmetry U(l) that corresponds to the baryon 
number conservation. In the limit of vanishing quark masses (which for mo- 
mentum transfers of q ~ 1 GeV/c is a good approximation for the light quarks 
u and d, and to a lesser extent also for the strange quark) the Lagrangian re- 
veals another symmetry under global vector and axial vector transformations 
in the SU(3) flavor space which are defined as: 

A^ ' A^ 

ip —> e~ ia v~2~i/j and ip — > e~ ta ' A ~'y 5 'ilj. (1.5) 
The conserved Noether currents associated with these symmetries are 

A* A* 
3v,i = irf-^b and 3% = ^7 M y7s^- (1.6) 

The quark-spinors in Eq. (1.3) can be decomposed into a left- and a right- 
handed component 

^=^(1^75)^ (1.7) 
which transform under the rotations defined in Eq. (1.5) as: 

ip L ^e- ia ^ip L and ip R ^ tp R (1.8) 

A* 

e- ia K-ip R and ^ L -> ip L , (1.9) 

which constitutes a global SU(3)l <8> SU(3) R chiral symmetry in flavor space. 
This symmetry conserves the handedness, i. e., the projection of the spin on 
the momentum direction, of a quark. Therefore the two kinds of quarks, left- 
and right-handed, do not mix dynamically. 

Spontaneous Breaking of Chiral Symmetry 

The observation of the mass splitting of chiral partners, e.g., the p and a\ 
mesons have a mass difference of ~ 550 MeV/c 2 (m p = 776 MeV/c 2 and 
m ai = 1230 MeV/c 2 [29]), implies a spontaneous breaking of chiral symmetry 
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due to a non-vanishing vacuum expectation value of the quark condensate 
{ipip) 7^ 0. While the vector current jv = Jl + 3r is still conserved the axial- 
vector symmetry (ja = Jl—Jr) is spontaneously broken, which means that the 
axial- vector charge Q k A = J d 3 xip^ ^75^ still commutes with the Hamiltonian 
but the ground state has a non-zero expectation value: Q\\0) 7^ 0. An often 
used analogy is a ferromagnet below the Curie temperature, in which the 
ground state of aligned spins breaks the rotational symmetry. The choice of 
a particular ground state leads to a spontaneous breaking of chiral symmetry 
which according to the Goldstone theorem [30] results in the appearance of 
eight massless Goldstone bosons (ti ± , 7r°, K ± , K°, K°, and rj) 

Such a system can be visualized by the following potential 

V = --n 2 (a*a + 7r + 7r) + ^V t(7 + ^nf (1.10) 

which is shown for /z 2 > and fi 2 < in Fig. 1.8 [31]. For /1 2 > the potential 
has a rotational symmetric ground state at a = n = 0. For /x 2 < 0, the state 
a = 7T = is a local maximum and therefore unstable. An infinite number 
of degenerate ground states lie on the circle a 2 + tt 2 = fi 2 /v 2 , which are not 
rotational symmetric. Therefore, picking a ground state spontaneously breaks 
the symmetry. However, effects of the symmetry are still present, as excitations 
around the symmetry axis, i. e. axial- vector rotations, do not cost energy. 
With the identification of the fields 7? = i^jXj^^p and a = these rotations 
correspond to the massless Goldstone bosons. In contrast, radial excitations 
along the a field do cost energy and correspond to massive particles. 

The strength of the symmetry breaking is measured by the vacuum expec- 
tation value of such Goldstone bosons which is, e. g., for pions: 

(o\jl k (x)Mp))^iS jk Up^. (1.11) 

with the measured pion decay constant / w = 93 MeV, which is the order 
parameter of the chiral symmetry. A second order parameter is the vacuum 
expectation value or quark condensate 

(^) = (0|^^ + <Ml|0> = (0\uu + dd\0). (1.12) 

The relation between quark condensate (ipip) and f n is given by the Gell- 
Mann-Oakes-Renner relation (GOR) [32]: 

m lfl = -m(ipip) 

= -2m(qq) (1.13) 
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where m = \{m u + m^) and (qq) = (uu) = (dd). 
Explicit Breaking of Chiral Symmetry 

The small but non-vanishing mass of the pseudoscalar mesons is explained by 
the explicit breaking of the chiral symmetry due the finite quark masses. 

The finite quark masses create a contribution —fn(ipip) in the QCD La- 
grangian in Eq. (1.3). This leads to a Mexican hat potential which is slightly 
tilted towards the positive a direction which breaks the symmetry of the po- 
tential: 

I i 
V = --fi 2 {a j a + ttV) + -v 2 ((J ] (T + vrV) 2 - ]^m\a (1.14) 

which is illustrated in Fig. 1.9. With the choice of fi 2 ' jv 2 = f 2 the tilting leads 
to a minimum at a = f n and axial-vector currents are no longer conserved, i. e. 
the pseudoscalar mesons, which were the Goldstone bosons of the spontaneous 
symmetry breaking, acquire a finite mass which is related to the current quark 
mass by Eq. (1.13). With m k 5 MeV/c 2 the vacuum value for the quark 
condensate becomes (qq) ~ —(254 MeV) 3 . 

In-Medium Quark Condensate 

So far only vacuum properties of hadronic matter have been discussed. Of par- 
ticular interest are the dynamics in the presence of a hot and dense medium as 
it is expected to be created in relativistic heavy-ion collisions. At high pressure 
or temperature hadronic matter is expected to undergo a phase transition to 
deconfined quarks and gluons accompanied by the melting of the quark con- 
densate. Already before reaching deconfinement, chiral symmetry is partially 
restored. 

The expected modifications can be derived in the limit of low temperatures 
T and density starting from the grand canonical partition function for a hadron 
gas in volume V in contact with a heat bath 

Z(V, T, n q ) = Tr ( e - (H -^ N)/T ) (1.15) 

with the Hamiltonian of the system H, the quark chemical potential and 
the quark number generator N. 

The expectation value of the quark condensate is given by the thermal 
average 

((qq)) = Z- 1 J2( n \^\ n ) e ~ iEn ~' iq)/T i ( L16 ) 

n 

where the sum is carried out over all eigenstates of H with the corresponding 
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Figure 1.8: Mexican hat potential as defined in Eq. (1.10). The case /i 2 > is 
shown in (a) has a symmetric ground state at a = tt = 0. For fi 2 < an infinite 
number of degenerate ground states lie on the circle around a 2 + tt 2 = [i 2 /v 2 , while 
a = 7r = is a local maximum. The choice of a ground state breaks the symmetry 
spontaneously. 
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V(G, 71=0) 



V(o, 71=0) 




(a) Symmetric Potential 



(b) Tilted Potential 



Figure 1.9: Mexican Hat Potential at it = 0. In case of (b) the potential is tilted 
along the a field breaking explicitly the chiral symmetry 



eigenvalues E n . For a non-interacting hadron gas the leading order result to 
Eq. (1.16) with respect to the vacuum expectation value for low densities is: 



(qq) 
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+ C(T 8 ). 



(1.18) 



The scale A g = 470 ±110 MeV is fixed from pion scattering data. To first order 
the quark condensate at low densities decreases quadratically with the tem- 
perature. Analog, the result for finite nuclear densities and low temperatures 
is a linear decrease with density: 



N) 



(qq) 



1.19) 



with baryon density /x/v = 3/x g . At nuclear matter densities (po = 0.17 fm -3 ) 
the quark condensate has dropped by 35% with a experimentally determined 
nucleon a-term of £ ^ ~ 45 MeV. Based on these calculations one can derive a 
prediction of the scalar quark condensate (qq) as function of temperature and 
density, which is shown in Fig. 1.10 [33]. 
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Figure 1.10: The expectation value of the chiral condensate as function of temper- 
ature T and nuclear matter density p as calculated with the Nambu-Jona-Lasinio 
model [33]. Figure based on [34]. 



Based on connections between hadron masses and the quark condensate 
various scenarios for the change of meson masses with decreasing quark con- 
densates are proposed. Brown and Rho predict in [35] a dropping of the p 
mass by 15-20% at normal nuclear matter densities and a universal Brown- 
Rho (BR) scaling of the in-medium vector meson masses with the in-medium 
pion decay constant /*: 

«W> fftY , ft m *rr m* m* m* 

^ZZlL = U-M and ^ = — = — = — = — (1.20) 

(ipip) \UJ Jtt rn a m N m p 

Other models explain medium modification of vector mesons by interactions 
with surrounding hadrons in the hot and dense medium. For an excellent 
review see, e. g., Ref. [3]. 



1.2 Dileptons 

Dileptons are produced during all stages of the collisions and carry a variety 
of signals. In this Section the various sources are discussed. 
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Figure 1.11: Feynman diagram for the production of a dilepton pair via the Drell- 
Yan process. 



In a collision of two nuclei, the Drell-Yan process is the annihilation of a quark 
in one nucleus with a sea antiquark from the other nucleus into a virtual 
photon which subsequently converts into a lepton pair. This process is shown 
in Fig. 1.11 for two nuclei A and B and has particular important at large 
invariant masses (m« > 3 GeV/c 2 ). At high energies the invariant mass of the 
lepton pair nin is given by the product of the quark momenta Xiy/s and £2\/2: 



with a/s the center of mass energy of the incoming nuclei. The cross section of 
Drell-Yan dilepton production can be calculated to leading order using parton 
distribution functions obtained from deep inelastic lepton-nucleon collisions. 
An additional factor K « 2 accounts for higher order corrections in a s is 
necessary to describe experimental dilepton data. 

One can show that in collisions of two equal nuclei with mass number A 
(e. g., for gold: Aj± u = 197), the contribution of dileptons from the Drell-Yan 
process scales with respect to the nucleon-nucleon CclSG clS A 4 / 3 [36]. In a 
thermalized partonic medium the same process, i. e., qq — > e + e~ , is possible 
between quarks and antiquarks with thermal momentum distributions. Theo- 
retical models predict such a contribution to be dominant in the intermediate 
mass region. 

1.2.2 Open Charm and Bottom 

Another contribution to the dilepton continuum are semi-leptonic decays of 
charm and bottom mesons. Heavy quark pairs (QQ = cc or bb, respectively) 



1.2.1 Drell-Yan 



rriu = X1X2S 
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Figure 1.12: Leading order processes of heavy quark production. 

are produced in nucleon-nucleon collisions in inelastic hard-scattering processes 
between constituent quarks of the two nucleons (qq — > g* — > QQ). The lowest 
order diagrams for heavy quark production are shown in Fig. 1.12. In addition 
to the annihilation of two light quarks into a virtual gluon, heavy quarks can 
be produced by interactions of a gluon of one nucleon with a gluon in the other 
nucleon (gg —> g* — ► QQ). 

Open charm and bottom mesons are formed out of a heavy quark and 
a light antiquark (u, d, or s) or a heavy antiquark and a light quark, e. g. 
D + = cd, B + = lib. After production of a heavy quark pair, the quarks 
undergo fragmentation and can form a D + D~ pair. Due to the large mass of 
the c quark, the D + D~ pair preserves most of the initial correlation of the QQ 
pair. The charmed (or bottom) mesons then can undergo weak decays into, 
e. g., D + — > K The total branching ratio of semi-leptonic decays of a D 
meson is on the order of 10%. The semi-leptonic decays of both the D + and 
the D~ results in the creation of a dilepton pair. 

Perturbative QCD calculations to leading order are not able to fully de- 
scribe charm production in nucleon-nucleon collisions. Similar to the Drell-Yan 
process, a K factor has to be employed to correct for the difference to the mea- 
sured charm cross section. Higher order calculations such as next-to-leading 
order (NLO), for which exemplary processes are shown in Fig. 1.13, and fixed- 
order-plus-next-to-leading order (FONLL) are in agreement with the D meson 
cross section measured by CDF in pp collisions at \fs =1.96 TeV [37] as well as 
with measurements of single electrons [16] and muons [38] from semi-leptonic 
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Figure 1.13: Examples of next-to-leading order processes contributing to the pro- 
duction of heavy quark pairs. 



charm decays performed by PHENIX in p + p collisions at y/s = 200 GeV. 
However, the theoretical uncertainties are large and the data prefer larger cross 
sections within these uncertainties. 

Shown in Fig. 1.14 is the invariant differential cross sections of electrons 
from heavy flavor decays [16] in comparison to a FONLL calculation [39]. The 
shape of the FONLL is in good agreement with the data, overall the data are 
a factor of ~ 1.7 higher than the calculation, as visible in the bottom panel of 
the same figure, which shows the ratio of data to FONLL. 

In the absence of medium modifications, the production of particles by 
inelastic hard scattering processes in heavy ion collision is given by the pro- 
duction cross section in p + p collisions scaled by the nuclear overlap fac- 
tor Tab, which is the integral over the product of the thickness functions of 
the two colliding nuclei A and B in the geometric overlap region (or Taa in 
case of two identical nuclei A). There is a simple relation between the nu- 
clear overlap factor, the inelastic p + p cross section and the average number 
of binary nucleon-nucleon collisions occurring in a nucleus-nucleus collisions: 

TaA = (Ncou) /(7pp. 

In order to quantify the deviation of the yield measured in a nucleus-nucleus 
collision from the binary scaled p + p expectation the nuclear modification 
factor Raa is defined as: 

r aa = dNAA = dNAA (i 22) 

T A a da p P ( Neon ) dNpp 

were dNAA (dN pp ) is the differential particle yield measured in the nucleus- 
nucleus (p + p) collision, and da pp = a pp dN pp . 

In Au + Au collisions at ^/snn = 200 GeV the single electron spectra 
from heavy flavor decays are strongly modified with respect to the expected 
yield form a hard probe [15]. Their nuclear modification factor is shown in 
Fig. 1.15. At high pt single electrons from heavy flavor decays are suppressed 
with respect the binary scaled yield of electrons in p+p collisions. Furthermore, 
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p T (GeV/c) 

Figure 1.14: (a) Shown is the differential cross section of electrons from heavy 
flavor decays in p + p collisions at y/s = 200GeV [16]. The data are compared to the 
central value of a FONLL calculation [39]. (b) Ratio of data to FONLL calculation. 
The upper (lower) curve show the upper (lower) theoretical limit of the FONLL 
calculation. 
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Figure 1.15: (a) Shown is the nuclear modification factor Raa of electrons from 
heavy flavor decays, (b) Elliptic flow of electrons form heavy flavor decays. A 
comparison to various models of medium effects is shown [15]. 

they exhibit significant elliptic flow, measured as the second Fourier coefficient 
V2 in the expansion of the azimuthal particle distribution, shown in the bottom 
panel of Fig. 1.15. Both these effects indicate that charm quarks are subject to 
significant energy loss and may thermalize in the medium created in Au + Au 
collisions at ^s^n = 200 GeV. This should leave an imprint also in the e + e~ 
pair spectrum in the IMR, where dielectrons from open charm decays are the 
dominant contribution. 

1.2.3 Direct Photons 

Direct photons are another important electromagnetic probe of the medium 
created in heavy ion collisions. As discussed below, direct photons are not 
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Figure 1.16: Feynman diagrams for direct photon production to leading order. The 
top row represents the annihilation process, the bottom two graphs the Compton 
process. 

only created as massless real photons, but also as virtual photons with non- 
zero mass, which internally convert into a dilepton pair. 

Direct photons are produced by inelastic scattering processes between the 
incoming partons. The lowest order processes as shown in Fig. 1.16 are quark- 
antiquark annihilation into a gluon and a photon and quark-gluon Compton 
scattering producing a quark and a photon. 

The production of direct photons in p + p collisions at = 200 GeV has 
been measured by PHENIX [40]. The measured direct photon cross section 
is in excellent agreement with NLO pQCD calculations [41-44] as shown in 
Fig. 1.17. The three curves which are shown correspond to three different 
choices of the momentum scale \i = O.bpx, l.Opr, 2.0pr- There are actually 
three different scales involved: the renormalization scale [lr, the factorization 
scale [ip, and the fragmentation scale p! F . The fragmentation scale is included, 
because the calculation includes direct photons from parton fragmentation into 
photons. All three scales are set to a common value [i = [lr = [lf = ^'f- 

In Au+ Au collisions, direct photons have been measured at high in vari- 
ous centrality bins [45] . The nuclear modification factor Raa for direct photons 
with pt > 6 GeV/c is shown in Fig. 1.18 as function of the participating nu- 
cleons A part . The production of direct photons is consistent for all centralities 
Taa scaled direct photons cross section measured in p + p, i. e., Raa = 1, i- e., 
the assumption that direct photons with transverse momenta above ~ 6 GeV/c 
are produced by inelastic scattering processes between the incoming partons 
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Figure 1.17: Direct photon section measured in p + p collisions at y/s = 
200GeV [40]. The data are compared with a NLO pQCD calculation [41-44]. The 
bottom panel shows the ratio of data to theory. 
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Figure 1.18: Nuclear modification factor of direct photons and 7r°'s with px > 
6.0 GeV/c as function of iVp art [45]. 



and escape the medium unperturbed by final state interactions. This observa- 
tion is supporting evidence that the observed n suppression [46], also shown 
in Fig. 1.18, is indeed a final state effect due to energy loss of the partons in 
the hot and dense medium. At very high px a deviation from one is observed 
which may, while still under investigation and despite large uncertainties, be 
attributed to isospin differences between the proton and a gold nucleon [47]. 

As direct photons carry the information about the momentum distribution 
of the partons involved in their production, the momentum distribution of di- 
rect "thermal" photons produced by partons in a thermalized medium directly 
reflects the temperature of the medium. 

In Au + Au collisions at RHIC energies, thermal photons are predicted to 
be the dominant source of direct photons in 1 < p T < 3 GeV/c [48]. This is 
illustrated in Fig. 1.19 which shows the predicted contributions to the total 
direct photon spectrum; from the initial hard scattering dominating the direct 
photon yield for px > 3 GeV/c, and the relatively soft thermal emission of 
from the hadron gas being the main contributor at px < 1 GeV/c, leaving 
a window for QGP radiation at 1 < px < 3 GeV/c. As initial conditions a 
formation time of To = 0.33 fm and a initial temperature T = 370 MeV were 
used. 

However, the inclusive photon yield is dominated by a large background of 
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Figure 1.19: Direct photon emission from central Au+Au collisions at ^/s/vtv = 200 
GeV as function of transverse momentum. The short dashed line shows the predicted 
contribution of photons form the initial hard scattering processes, the dashed-dotted 
line thermal radiation from the QGP and the long-dashed line emission from the 
hadron gas. [48]. 
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hadron decays such as ir° — > 77 for p? < 5 GeV/c, limiting the measurement 
of direct photons at low with the Electromagnetic Calorimeter (EMCal) in 
PHENIX, which is based on a statistical subtraction of the decay photon back- 
ground. The uncertainties in the knowledge of the decay background can be 
reduced by avoiding its explicit measurement, but rather tagging of the decay 
photons [49]. Furthermore, to circumvent the limitations due to the energy 
resolution at low photon energies, the excellent capabilities of the PHENIX 
detector to measure electrons have been utilized by measuring photons via 
their external conversion into e + e~ pairs, a method discussed in Appendix B. 



In contrast to massless real photons, virtual photons bring an additional 
observable, their invariant mass, which as discussed in the following brings 
the advantage of measuring thermal photons in a better signal to background 
region than real photons. In general any source of real photons, e. g. the ones 
shown in Fig. 1.16, can also create a virtual photon which is emitted as e + e~ 
pair. The relation between photon production and the associated e + e~ pair 
production can be written as [50]: 



d 2 n e 
dm P 




Sdn^ (1.23) 



with m ee being the invariant mass of the e + e~ pair and m e = 511 keV/c 2 the 
mass of the electron. The process dependent factor S aproaches 1 as m — > or 
m <C pr- For n° and 77 decays, S is given by S = |F(m 2 e )| 2 (l — m^/m^) 3 [51] 
where is the hadron mass and F(m 2 e ) the form factor. For e + e~ pair 
masses approaching nih, the factor S goes to zero. While the measurement of 
real thermal photons suffers from a large background of hadron decays (80% 
of the background comes from 7r° decays), measuring virtual photons allows 
to select a mass range m ee > m n o = 135 MeV/c 2 , in which the signal to 
background ratio is improved by a factor five, thus allowing a 10% signal of 
direct photons to be observed as a 50% excess of e + e~ pairs. 



For quark-gluon Compton scattering the exact relation between the photon 
production process and the e + e~ pair process can be calculated as [52]: 

d 3 n ee 2a 1 / 4m 2 / 2m 2 \ , 2u ,dn~, 

(1.24) 



dm ee dt 




25 




Figure 1.20: Vector dominance model in dilepton production via ir + ir annihila- 
tion. 



in which the Mandelstam variables s, t, and u are defined as: 

S = (Pg+Pgf 
t = (Pq~ P~,Y 
U = {Pq-P'af 



(1.25a) 
1.25b) 
(1.25c) 



with the four-momentum vectors of the incoming quark p q , incoming gluon p g , 
the outgoing real or virtual photon p g , and the outgoing quark p' . Comparing 
Eq. (1.24) with Eq. (1.23), the process specific S factor for the quark-gluon 
Compton process can be identified as: 



5 = 1 + 



2u 



t 2 + s 



2 m ee 



1 - 



2x 



1.26) 



(x + Vl + x 2 )(3x 2 + 1 + 2x Vl + x 2 ) 
with x = m ee /pT which for px 3> m ee simplifies to S ~ 1. 



1.2.4 Medium Modifications of Vector Mesons 

At later stages of the collision, after the medium has expanded and cooled down 
below the critical temperature, quarks and gluons are confined to hadrons in a 
hadron gas. The dominant dilepton production during this stage is expected 
from pion and kaon annihilation and scattering between other hadrons. The 
cross section for irn scattering is dynamically enhanced through the formation 
of light vector mesons p, u, and <fi. As these vector mesons carry the same 
quantum numbers as a photon, they couple according to the vector dominance 
model (VDM) directly to a lepton pair: ti + ix~ —> p —> 7* — > [53]. Thus 
the invariant mass of the lepton pair directly reflects the mass distribution 
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Figure 1.21: In-medium spectral function of the vector mesons for different tem- 
peratures and densities [54, 55]. 

of the vector meson at the time of its decay. As the lifetime of the p meson 
(r ~ 1.3 fm) is shorter than the lifetime of the hadron gas ~ 10 fm, most of the 
p mesons will decay in-medium imprinting any medium effect on the lepton 
pair. The invariant mass distribution of e + e~ pairs from p decays is given by 
the p spectral function (modulo a thermal Bose factor fs). Schematically p 
spectral function as function of its energy u and 3-momentum q, defined as 
the imaginary part of the propagator, can be presented in the form of: 



with the p the pole mass m p . The p self-energy S p includes the summation of 
all scattering amplitudes which includes in-medium interactions with baryons 
and mesons which broaden the spectral function. 

As an example, a prediction of the spectral functions of p, to, and <fi at 
RHIC energies by R. Rapp [54, 55] is shown in Fig. 1.21. Recent reviews of 
the role of the p meson in dilepton emission can be found, e. g., in Refs. [3, 
56-59]. In this calculation the p and u spectral function shows the strong 
broadening towards higher temperatures and densities, %. e., towards the phase 
boundary. The slight upward shift in mass is due to repulsive parts in the self 
energy, e. g., from baryonic particle- hole excitations. The does not show such 
strong modifications, but at the highest temperatures the width is significantly 
increased by a factor of ~ 7. 

This explains the special role of the light vector mesons, and their in- 
medium modifications, in dilepton measurements. In contrast, heavier vector 
mesons such as the J / ip or the T have a substantially longer lifetime and decay 
predominantly after freeze out. Also they carry information of modifications 
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to heavy quarks, e. g., suppression, but it is in the magnitude of their signal 
rather than their spectral shape, see e. g. Ref. [60]. 

After the thermal freeze out temperature is reached, the dominant sources 
of lepton pairs are resonance and Dalitz decays of light mesons with their 
vacuum properties, such as 7r°, 77, 7/. to, and </>. 

1.3 Previous Experimental Results 

A number of experiments have measured the dilepton continuum in heavy ion 
collisions under a variety of conditions. The chapter gives a brief summary of 
their results, but is by no means meant to be complete. Section 1.3.1 sum- 
marizes the recent results of the HADES experiment at the SIS accelerator 
at GSI. Section 1.3.2 presents the result of the NA45/CERES collaboration 
at the SPS accelerator at CERN. The measurement of the dimuon contin- 
uum by the NA60 collaboration is summarized in Section 1.3.3. The section 
concludes with the direct photon measurement of the WA98 collaboration in 
Section 1.3.4. 

1.3.1 HADES 

The High Acceptance DiElectron Spectrometer (HADES) is a fixed target ex- 
periment at SIS accelerator at GSI. It has measured the dielectron continuum 
up to m ee = 1 GeV/ c 2 in low energy collisions of light ions, such as C-C at 1 and 
2 AGeV [61, 62] as shown in Fig. 1.22. The yield of e + e~ pairs is compared to 
a cocktail of decays of 7r° — > ^e + e~, 77 — > 7e + e~, u — > 7i°e + e~, and uo — * e + e~ . 
While the mass region below 150 MeV/c 2 is well explained by the cocktail, an 
enhanced yield of e + e~ pairs in the mass region of 0.15-0.50 GeV/c 2 above 
the cocktail is observed. Also including in addition p and A resonance de- 
cays does not fully explain the measured yield. HADES also compared the 
yield of e + e~ pairs to the results of the DLS experiment [63] and finds them 
in good agreement [61, 64] challenging current theoretical descriptions of the 
dielectron production at these energies. The excess yield scales as the ir° yield 
with increasing beam energy, rather than the 77 yield [61]. From preliminary 
results [65] from elementary reactions (p — p, n — p) it has been shown that a 
isospin dependent Bremsstrahlung contribution, increased by a factor of ~ 4 
with respect to early theoretical calculations, can almost entirely explain the 
excess observed in C-C collisions at the same energy [66, 67]. This clearly 
shows the importance of a baseline measurement in elementary collisions at 
the same energy as the heavy-ion collisions. 
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(a) C-C at 1 AGeV [61] (b) C-C at 2 AGeV [62] 



Figure 1.22: Invariant mass spectrum of e + e~ pairs measured by HADES in C-C 
collisions at 1 ^4GeV (a) and 2 AGeV (b), respectively. The measured yields are 
compared to a cocktail of e + e~ pairs from decays of it , n and uj (cocktail A, solid 
line). A second cocktail, including A and p resonance decays is shown as dashed 
line (cocktail B). The bottom panels show the ratios of data and cocktail A. 
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1.3.2 NA45/CERES 



The NA45 experiment 4 better known as ChErenkov Ring Electron Spectrom- 
eter (CERES) was a fixed target experiment at the Super Proton Synchrotron 
(SPS) at CERN. It has measured the dielectron continuum in heavy-ion reac- 
tions at kinetic beam energies of 40-200 AGeV [68-70]. 

While the e + e~ pair yield in proton induced collisions such as p-Be (as 
shown in Fig. 1.23a) and p-Au at 450 AGeV is fully reproduced by a cocktail 
of hadron decays, an enhanced yield for m ee > 200 MeV/c 2 has been observed 
in S-Au collisions at 200 AGeV which is shown in Fig. 1.23b [68]. CERES has 
also measured a low mass enhancement of e + e~ pairs in Pb-Au collisions at 
40 and 158 AGeV. At the latter energy the pt and centrality dependence of 
the enhancement have been studied, which was found to be localized at low 
Pt and to increase stronger than linear with the charged particle density [70] 
suggesting binary annihilation processes such as n + 7r~ — > p — > 7* — > e + e~. 
Their results are compared to model calculations of three scenarios, (i) free 
7T7T annihilation in vacuum, (ii) assuming Brown- Rho scaling of the in medium 
p meson [71], and (iii) a modified in-medium p spectral function [3]. The 
models have employed a thermal fireball model to include the full time evo- 
lution of the system which includes the experimentally determined freeze out 
conditions (T,p B ) io = (115 MeV, 0.33p ) and initial conditions of (T,p B ) init = 
(190 MeV, 2.55po) ; as well as a finite pion chemical potential. For the dropping 
mass scenario a in-medium mass of: 



with C = 0.15, T c x = 200 MeV and a = 0.3 has been assumed. Clearly, 
the in vacuum im annihilation fails to describe the data, it gives too much 
yield at the u resonance peak and too little in the mass range 0.2 < m ee < 
0.7 MeV/c 2 , which leads to the conclusion that in-medium modifications must 
play a role. Within the experimental uncertainties both the dropping mass as 
the in-medium collisional broadening scenarios give a reasonable description 
of the enhanced dielectron yield. 



The NA60 experiment at SPS was created out of the NA50 experiment 5 [73] 
upgraded with a silicon- vertex tracker. NA50 had measured p + p~ pairs, rather 



4 named after its location in the North Area of the SPS, in contrast to the West Area 
5 itsclf an upgrade of the NA38 experiment [72] to study Pb-Pb collisions 




(1.28) 



1.3.3 NA60 
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m (GeV/c 2 ) m (GeV/c 2 ) 

(a) p-Be at 450 AGeV (b) S-Au at 200 AGeV 

Figure 1.23: Invariant mass spectrum of e + e~ pairs measured by CERES in p-Be 
collisions at 450 AGeV (a) and S-Au collisions at 200 AGeV (b), respectively. The 
measured yields are compared to a cocktail of e + e~ pairs from hadron decays; the 
various sources are shown individually. While the cocktail agrees with the measured 
yield of e + e~ pairs in p-Be collisions, an enhancement above the cocktail is observed 
in Pb-Au collisions [68]. 



than e + e~ pairs eliminating the major background source of 7r° Dalitz decays 
due to the larger muon mass (m^ = 105.66 MeV/c 2 ). For Pb-Pb collisions 
at 158 AGeV, they reported an enhanced dimuon yield in the intermediate 
mass region of 1.15 < m MM < 2.56 GeV/c 2 of up to a factor of « 1.65 above 
the expected sources from open charm decays and Drell-Yan [74]. They were 
able to rule out an enhanced Drell-Yan production, as the mass shape of the 
enhancement was much steeper then expected for Drell-Yan and concluded the 
enhancement was consistent with a charm production enhanced with respect 
to the one observed in p-A collisions, but did not excluded the contribution of 
thermal radiation. 

The upgrade with a silicon-vertex tracker allowed NA60 to determine the 
distance between the muon tracks and the collision vertex. Due to the long 
life-time of D mesons, their decays vertex would be offset from the collision 
vertex in contrast to a prompt source of n + fi~ pairs [75]. Fig. 1.25a shows 
the invariant mass spectrum of excess pairs compared to the expected 

shapes from open charm decays and Drell-Yan [76]. The shape is compatible 
with semi-leptonic decays of charmed mesons. Fig. 1.25b shows the weighted 
offset distribution of all fi + fi~ pairs. The weight in the pair offset from 

the vertex considers the momentum dependent resolution of the vertex mea- 
surement of a muon. The distribution of pairs is fitted with the expected 
shapes from open charm and prompt decays. The result is compatible with 
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Figure 1.24: Invariant mass spectrum of e + e~ pairs measured by CERES in Pb-Au 
collisions at 158 ^4GeV. The measured yield is compared to a cocktail of e + e~ pairs 
from free hadron decays and (i) without p decays (thin line), (ii) a vacuum p spectral 
function (thick dashed line) [71], (iii) a dropping in-medium p mass (thick dashed- 
dotted line) [3], (iv) a medium-modified p spectral function (thick solid line) [70]. 
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Figure 1.25: Shown in (a) is the mass distribution of the acceptance corrected ex- 
cess yield compared to the expected shapes from open charm and Drell-Yan decays. 
A fit of prompt and charm decays to the weighted offset distribution of fi + /-i~ pairs 
in the intermediate mass region is shown in (b) [76]. 



an open charm yield as extrapolated from the NA50 p-A result and a prompt 
excess with an enhancement factor of 2.4 over the Drell-Yan contribution. 

Furthermore, NA60 has observed an enhancement in the low mass dimuon 
continuum in In-In collisions at 158 AGeV shown in Fig. 1.26a [76-78]. The 
high precision of the data allowed them to subtract all sources of hadron decays 
(besides the contribution of p decays) and extract the enhancement. The two 
scenarios which were in reasonable agreement with the CERES data as shown 
in the previous chapter, an in-medium broadened p spectral function [3] and 
the dropping mass [71] are compared to the excess yield in Fig. 1.26b. While 
the broadening scenario is in good agreement, the Brown- Rho scaling does not 
explain the observed enhancement. 

The rriT spectra of the excess yield is shown in Fig. 1.27a for four slices in 
mass and the <p. They are fitted to an exponential and the inverse slope pa- 
rameter T e ff is shown in Fig. 1.27b together with the results of the same fitting 
procedure of the excess yield observed in the intermediate mass region. The 
inverse slope of the low mass enhancement follows closely the trend observed 
for the hadrons rj, u> and <fi which is consistent with the expected linear increase 
due to radial flow indicated by the solid line T eff = T fo + m((3 T ). The effective 
temperature of the vacuum p is much higher, which is interpreted as a later de- 
coupling of the p from the medium and therefore obtaining a larger radial flow. 
Above the <j) meson the effective temperature suddenly drops and stays con- 
stant with mass. This behaviour is consistent with an thermal emission from 
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Figure 1.26: Shown in (a) is the invariant mass spectrum of p + p~ pairs measured 
by NA60 in In-In collisions at 158 ^4GeV. The measured yield (open circles) is 
compared to individual sources of pairs from hadron decays without any p 

contribution. Their contributions are subtracted from the data, and the difference 
(solid triangles) attributed to p decays [76, 78]. In (b) the excess yield is compared 
to a number of theoretical predictions in [76]. 



the early partonic phase, before significant radial flow has developed [79], but 
hadronic scenarios are not excluded [80]. 

At very low ttit (ttlt — M < 0.2 GeV/c 2 ) a steepening of the tut spectra 
shown in Fig. 1.27a is observed for all four mass windows of the excess yield. 
This trend is opposite to the expectation for radial flow and is not seen for the 
(p meson. 

1.3.4 WA98 

The WA98 Collaboration has measured direct photons in Pb-Pb collisions 6 [81] 
based on a statistical subtraction of decay photons, which show an excess above 
the yield expected from the direct photon measurements in proton induced 
reaction. The result is shown in Fig. 1.28 in comparison to a prediction of 
direct photons [48], which includes a contribution of thermal emission from 
a QGP phase; rs 30% of the total thermal photon yield. This prediction is 
based on the same fireball model which successfully describes the low mass 
enhancement observed in the dilepton continuum in CERES and NA60 with 
a formation time of r = 1 fm and an initial temperature of Tj ~ 210 MeV. 



6 with an array of lead glass calorimeter which is now installed in PHENIX 
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(a) (b) 

Figure 1.27: Shown in (a) are the spectra of the excess yield of pairs in 

four mass windows in comparison to the <j> summed over all centralities excluding 
the most peripheral bin. In (b) the inverse slopes T e g are shown as function of the 
invariant mass. The inverse slopes are extracted from a fit of the mj- spectra in (a) 
(open symbols) and of narrower mass windows (filled triangles). The inverse slopes 
are compared to the slopes of hadrons (open circles) [76, 78]. 

The two recently published low px data points [82], extracted via photon HBT 
(Hanbury Brown and Twiss Effect [83, 84]) interferometry, are not described 
by current theoretical models which attribute this yield to the hadronic stage 
of the fireball [48, 85]. But the inclusion of soft Bremsstrahlung off irn and 
Ti K scattering improves the situation [86]. 
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Figure 1.28: Yield of direct photons in Pb-Pb collisions at 158 AGeV measured 
by WA98. The lowest to points are extracted from the strength of the two-photon 
correlation [82], the other points with a statistical subtraction method [81]. 
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1.4 RHIC 



The Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Labo- 
ratory (BNL) provides heavy-ion collisions at the highest energies currently 
available. It provides collisions of p + p, d + Au, Cu + Cu, and Au + Au 
with energies up to y / J/v r /v = 200 GeV, but can also vary the collision ener- 
gies to close the gap to SPS energies on the search for a critical point in the 
QCD phase diagram. Besides the study of the quark-gluon plasma, the po- 
larized p + p collisions are of particular interest for the study of the proton 
spin structure and allows to measure the gluon's contribution to the spin of 
the proton [87]. For the near future RHIC plans to provide polarized p + p 
collisions at y^s = 500 GeV, energetic enough to produce W bosons, which 
will allow to probe u and d quarks independently. 

The ions are accelerated in two intersecting rings, one clockwise, the other 
counter-clockwise and brought to collisions at the six ring intersections. Four 
of the intersections have been equipped with experiments. The BRAHMS 
experiment was designed to measure charged hadrons over a wide range of 
rapidity and transverse momentum. It completed its data taking program in 
2006. PHOBOS was equipped with subsystems to measure charge particle 
multiplicities over almost the entire solid angle, and in addition with two 
magnetic spectrometers providing particle identification in a narrow aperture. 
Phobos completed its running time at RHIC in 2005. The Solenoid Tracker 
At RHIC (STAR) with its large acceptance Time Projection Chamber (TPC) 
covering the full azimuth and \y\ < 1.5 is a multipurpose detector with focus 
on global event analyses, particle correlations and particle identification. 

Reviews of the results of the first three years of RHIC operation by all four 
collaborations can be found in their "white papers" [8-11]. 
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Chapter 2 

The PHENIX Experiment 



The PHENIX experiment 1 is a detector system which consists of four spec- 
trometer arms and two sets of global detectors. Of the four spectrometer 
arms, two, the so called central arms, are located at mid-rapidity covering 
each rj < 0.35 in pseudo-rapidity and ir/2 in azimuth. The two arms at for- 
ward rapidity, the muon arms, cover both 2ir in azimuth. The north arm covers 
the pseudo-rapidity range of 1.15 < rj < 2.44 while the south arms extends 
over the pseudo-rapidity range —1.15 > rj > —2.25. 

The detector setup is shown in Fig. 2.1; the upper panel shows the two 
central arms in a cut-away view perpendicular to the beam direction, while 
a view along the beam direction is displayed in the lower panel. The latter 
one also shows the two Muon Arms at forward rapidity. The azimuthal and 
rapidity coverage of each subsystem is summarized in Tab. 2.1. 

Two global detectors, the beam-beam counters (BBCs) and the zero-degree 
calorimeters (ZDCs), serve as event trigger and are responsible to measure 
global event parameters like the collision time, vertex and centrality. They are 
described in Section 2.1. 

The two central arms provide tracking and momentum measurements of 
charged particles over a large range in px from 0.2 GeV/c to 20 GeV/c, as 
well as particle identification; electrons are identified via the Ring Imaging 
Cherenkov Counter (RICH) and the energy-momentum matching measured 
in the Electromagnetic Calorimeter (EMCal) and the Drift Chamber (DC), 
respectively, photons via their electromagnetic showers in the EMCal and 
hadrons via their time of flight. All detector subsystem that have been used 
for the measurements presented in this thesis are described in Section 2.2. 

1 A variety of explanations for the name PHENIX exist, from the a posteriori definition 
as Pioneering High-Energy Nuclear Interaction experiment to the experiment which rose 
from the ashes of the original proposals for the RHIC experiments TALES, SPARC, OASIS, 
and DIMUON. 
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West Beam View East 




Figure 2.1: The PHENIX Detector configuration as of 2004. The upper figure 
shows the two central arms in a view perpendicular to the beam axis. The lower 
one shows a cut away view along the direction of the beam. 
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Table 2.1: Summary of the PHENIX detector subsystems [88]. 



Element 



Ar] 



Purpose and special 
features 



Magnet 
Central (CM) 
Muon (MMS) 
Muon (MMN) 
Beam-beam Counters 
(BBC) 
Zero-degree 
Calorimeter (ZDC) 
Drift Chambers (DC) 



Time Expansion 
Chambers (TEC) 
Ring Imaging 
Cherenkov Counter 
(RICH) 

Time of Flight (ToF) 
EMCal 

Lead-Scintillator 
(PbSc) 

Lead-Glass (PbGl) 



Muon Tracker (MuTr) 
MuTr South 
MuTr North 
Muon Identifier 
(MuID) 

MuID South 
MuID North 



±0.35 

-1.1 to -2.2 
1.1 to 2.4 
± (3.1 to 3.9) 

±2 mrad 

±0.35 



Pad Chambers (PC) ±0.35 



±0.35 
±0.35 

±0.35 

±0.35 

±0.35 



-1.15 to -2.25 
1.15 to 2.44 



-1.15 to -2.25 
1.15 to 2.44 
40" 



2tt 
2tt 
2tt 
2tt 

2tt 

2 x tt/2 

2 x vr/2 
tt/2 

2 x vr/2 
tt/4 

7r/2±7r/4 
tt/4 



2tt 
2tt 



2tt 
2tt 



Up to 1.15 Tm 
0.72 Tm for rj = 2 
0.72 Tm for ?] = 2 
Start timing, fast 
vertex 

Minimum bias trigger 

Good momentum and 
mass resolution 
a m /m = 1% at 
m = 1 GeV 
Pattern recognition, 
tracking for non-bend 
direction 

Pattern recognition, 
dE/dx 

Electron identification 



Hadron identification, 
a < 100 ps 

Electron and photon 
identification and 
energy measurement 
e^/ir^ separation at 
p > 1 GeV/c by EM 
shower and 

p < 0.35 GeV/c by ToF 
K ± /7r ± separation up 
to 1 GeV/c by ToF 
Tracking for muons 



Steel absorber and 
Iarocci tubes for 
/i/hadron separation 



For a description of the other central arm detectors, i. e. the Time Ex- 
pansion Chamber (TEC), the Time of Flight detector (ToF), as well as the 
muon arms, in which two detectors (MuTr and MuID) allow tracking and 
identification of muons at forward rapidity, one is referred to Refs. [89-91]. 

2.1 Global Detectors 

Global detectors are used to measure the event topology, i. e., the vertex posi- 
tion, the orientation of the reaction plane, and the centrality of the collision. 
Furthermore they measure the time at which a collision occurs. There are two 
pairs of detector systems installed on either side of the interaction point which 
are responsible for event selection and characterization in heavy ion collisions: 
two sets of beam beam counters cover 3.1 < \rj\ < 3.9 and two sets of zero 
degree calorimeters are installed at rj ~ ±6.9. 

2.1.1 Beam-Beam Counters 

The major tasks of the Beam-Beam Counters (BBC) [92, 93] are to serve as 
a trigger for collisions at the interaction point and to provide time and vertex 
information of the collision. The BBC comprises of two identical sets of 64 
hexagonal shaped Cherenkov counters as shown in Fig. 2.2, which are installed 
around the beam pipe at a distance of ±144 cm on the north and south side 
of the interaction point. 

Measuring the time difference between the BBC North and the BBC South, 
allows for the determination of the collision time as well as collision vertex: 

*o = 2^ BBCS + ^-bbcn) (2.1a) 

c 

^vertex = ^(^BBCS — £bBCn) (2.1b) 

where £bbcs and £bbcn are the average arrival time of particles in the BBC 
South and BBC North, respectively. The time resolution of the BBC is 52±4 ps 
(rms). This corresponds to a vertex position resolution of 1.1 cm. 

2.1.2 Zero Degree Calorimeter 

All four RHIC experiments are equipped with a pair of Zero Degree Calorime- 
ters (ZDC) [94] located at a distance of 18 m downstream of each interaction 
point behind the first accelerator "DX" dipole magnet as shown in Fig. 2.3. 
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(a) (b) 

Figure 2.2: Components of the BBC. (a) shows a photomultiplier tube (PMT) 
with a diameter of 25.4 mm and a 30 mm thick Quartz window mounted in front of 
it. (b) shows 64 PMTs assembled to one unit. 

Their task is to measure the energy of spectator neutrons, which did not par- 
ticipate in the collision and therefore carry still a large fraction of the beam 
momentum. While charged particles like spectator protons are deflected by 
the "DX" dipole magnet in front, neutrons hit the ZDC and create a hadronic 
shower. Neutral particles created within the heavy ion collision moving in 
forward direction have typically a much smaller energy. The ZDC consist of 
Cherenkov sampling hadronic calorimeter made of a tungsten alloy with a con- 
ical coverage of 21 mrad around the beam direction. The energy resolution of 
the ZDC is a E /E = 85%/VE © 9.1%. 

The total energy deposited by spectator neutrons can be used in anti- 
correlation with the total charge deposited in the BBC to determine the cen- 
trality of the collision as shown in Fig. 2.4. The centrality is a measure for 
how much the two colliding ions overlap. In addition to the centrality deter- 
mination the ZDC also serves as part of the minimum bias trigger in heavy 
ion collisions and provides timing information, but with a resolution of « 200 
ps it is less accurate than the BBC. 
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Figure 2.3: Top: Cut-away view of the ZDC in the plane perpendicular to the 
beam axis indicating deflection of protons and neutrons downstream of the "DX" 
Dipole magnet. Bottom: Plan view of the collision region [94]. 
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Figure 2.4: Correlation between the total energy deposited in the ZDC and the 
total charge measured in the BBC for Au + Au collisions at y/s^N = 200 GeV. The 
colored regions show the definition of centrality classes based on this correlation 
(black: 0-10%, red: 10-20%, green: 20-40%, blue: 40-60%, and magenta: 60- 
92%). Their boundaries are perpendicular to the centroid of the distribution. 
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2.2 Central Arm Detectors 



The two Central Arms consist out of several subsystems for charged particle 
tracking, momentum measurement and particle identification. Each arm cov- 
ers \r)\ < 0.35 in pseudo-rapidity and \<j)\ < tt/2 in azimuth. The central arm 
coordinate system is chosen with its origin at the nominal interaction point 
such that the z-axis is aligned with the beam direction pointing north, the 
x-axis pointing west and the y-axis upwards perpendicular to the two other 
axis. In this coordinate system the west arm covers — -^n < < j^n and the 
east arm < (f> < j^n. 

2.2.1 Central Magnet 

The transverse momentum of charged particles is determined by their bending 
in the magnetic field provided by the Central Magnet (CM) [95]. It consists 
out of an inner and an outer pair of concentric Helmholtz coils inside a steel 
yoke which provide a axially-symmetric magnetic field around the interaction 
point that is parallel to the beam direction as shown in Fig. 2.5. They cover 
a polar angle range of 70° < 6 < 110° which corresponds to a pseudo-rapidity 
range of \rj\ < 0.35. 

The two pairs of Helmholtz coils can be run with the same (++) or opposite 
(H — ) polarity. In the ++ configuration a total field integral of J Bdl = 1.5 Tm 
is achieved over the first 2 m from the interaction point, while the magnetic 
field in the region of the tracking devices (R > 2 m) is nearly zero in order 
to allow a tracking model which assumes straight tracks and to minimize the 
smearing of Cherenkov rings in the Ring Imaging Cherenkov Counter. 

The H — configuration leads to a cancellation of the field in the first 50 cm 
around the vertex to approximately zero field integral, which is used in combi- 
nation with the Hadron-Blind Detector (HBD) [96, 97], which was installed for 
the first time in 2007 inside this field free region and will be briefly discussed 
in Chapter 5. 

2.2.2 Drift Chambers 

A charged particle which traverses a gas-filled detector randomly ionizes the 
gas. The electrons from the primary ionization process are drifted in an electri- 
cal field towards an anode (sense) wire after a time proportional to the distance 
of the the track to the wire. Hits in subsequent anodes can be reconstructed 
to a track as described Section 3.2. 

Two Drift Chambers (DC) are installed in both Central Arms as the main 
tracking device for charged particles in PHENIX. They each consist out of a 
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Magnetic field lines for the two Central Magnet coils in combined (++) mode 

Figure 2.5: The PHENIX Central and Muon Magnets and their field lines shown 
in a cut-away view for adding (++) configuration [95]. Also shown are the two 
magnets in the Muon Tracker arms on either side of the Central Magnet. 

multiwire gas chamber located at a distance of 2.02 m to 2.40 m from the 
interaction point outside the magnetic field of the Central Magnet. The DC 
reconstructs the trajectory of charged particles in the r — <p plane in order to 
determine their transverse momentum p?. 

Both chambers extend over 2 m along the beam direction corresponding 
to Arj = ±0.35 in pseudo-rapidity; while the Drift Chamber installed in the 
Central Arm West covers — ^ix < 4> < j^tt in azimuth, the one in the East 
Arm covers < </> < j^ir. Their active volume is confined by Mylar windows 
and supported by a cylindrical shaped titanium frame as shown in Fig. 2.6. 
The detectors are filled with a gas mixture of 50% Argon and 50% Ethane. 

A design goal of the drift chamber was to measure the mass of the <j) meson 
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Figure 2.6: The PHENIX Drift Chamber Frame [89]. 

in the <p — > e + e~ decay channel with a resolution better than its natural width 
of 4.4 MeV/c 2 . In conjunction with the necessity to perform in a high particle 
multiplicity environment (as many as two hundred tracks in the central Au+Au 
collisions) this imposes the following requirements on the DC: 

• Single wire resolution better than 150 fim in r — <fi. 

• Single wire two track separation better than 1.5 mm. 

• Single wire efficiency better than 99%. 

• Spacial resolution in ^-direction better than 2 mm. 

Each Drift Chamber consists of 20 identical sectors covering 4.5°. As illus- 
trated in Fig. 2.7, each sector contains six different types of wire modules, XI, 
Ul, VI, X2, U2, and V2 stacked in radial direction. Every module contains, 
alternating in azimuth direction, four anode (sense) and four cathode planes. 



47 



The XI and X2 wires run parallel to the beam direction to track particles in 
the r — <p plane. Behind each X wire module two smaller U and V modules 
follow whose wires have a small stereo angle of 6° with respect to the X wires 
in order to measure the i-coordinate of the track. 

Each X module contains twelve sense wires separated by Potential (P) 
wires and surrounded by Back (B) and Gate (G) wires as shown in the left 
panel of Fig. 2.7 to shape the electrical field lines such that every sense wire is 
alternating sensitive to drift charges from only one side therefore limiting the 
left-right ambiguity to a region of ±2 mm. A calculation of the nominal drift 
field configuration due to this wire layout is shown in Fig. 2.8. The layout of the 
U,V-stereo modules is identical, but they contain only four sense wires. The 
stereo wires start in one sector and end on the other side of the Drift Chamber 
in the neighboring sector, as illustrated in the right panel of Fig. 2.7. 

For the pattern recognition to work with up to 500 tracks, each sense wire 
is separated in two halves at the center. Each half is read out independently. 
To electrically isolate the two halves, they are connected by a 100 /im thick 
Kapton strip. In total the Drift Chamber contains 6500 wires and therefore 
13000 read out channels. 

2.2.3 Pad Chambers 

The Pad Chambers are multiwire proportional chambers with a cathode pad 
readout that determine space points along the straight trajectory of charged 
particles to determine the polar angle 9 which allows to reconstruct the z- 
component of the momentum vector. 

The Central Arms are equipped with three layers of Pad Chambers in the 
West Arm and two layers in the East Arm, respectively. The first layer of Pad 
Chambers (PCI) is installed just behind the Drift Chambers, while the third 
layer (PC3) is situated right in front of the Electromagnetic Calorimeter. The 
second layer of Pad Chambers (PC2) is only present in the West Arm following 
the Ring Imaging Cherenkov Counter. 

Each PC contains a single layer of wires within a gas volume that is confined 
by two cathode planes. One cathode plane is solid copper, while the the other 
one is segmented into a fine array of pixels as shown in Fig. 2.9. The basic 
unit is a pad formed by nine non-neighboring pixels, which are read out by a 
common channel. Three pixels within a pad form a cell. For a valid hit, three 
neighboring pads must sense the avalanche. The interleaved design allows a 
fine position resolution of 1.7 mm in z direction and 2.5 mm in x and y and 
reduces at the same time the number of electronic channels by a factor of nine. 

Associating hits in PCI with tracks reconstructed in the DC is essential to 
determine the three dimensional momentum of a particle. 
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Wire orientations, top view 

V2 - planes 




Sector, side view 

Figure 2.7: Left: Cut-away view of the wire layout within one keystone of the 
Drift Chamber. Right: Plan view of the stereo wire orientation [89]. 
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Figure 2.8: Calculation of the drift lines for the nominal electrostatic field config- 
uration. Different wires are marked by letters: Back (B), Gate (G), Potential (P), 
and Sense (S) wires [89]. 
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Figure 2.9: Left: the pad and pixel geometry. Right: Interleaved pad design. [89]. 

2.2.4 Ring Imaging Cherenkov Counter 

A charged particle travelling in a medium with a velocity (3c that is greater 
than the speed of light in this medium, c n = c/n for a medium with refractive 
index n, emits Cherenkov radiation under angle cos 9c = l/(n/3). 

In each of the two central arms a Ring Imaging Cherenkov Counter (RICH) 
is installed between the inner and outer tracking detectors following the first 
layer of Pad Chambers [90]. Its main purpose is the separation of electrons 
from the large background of charged pions produced in heavy-ion collisions. In 
combination with the Electromagnetic Calorimeter it also provides information 
for an electron trigger in p + p collisions. Behind the entrance window with 
an area of 8.9 m 2 a volume of 40 m 3 is filled with CO2 as radiator gas, which 
has a refractive index of n — 1 = 410 x 10~ 6 [29], corresponding to a threshold 
velocity (3 t = l/n = 0.99590168 and a 7-factor of 7* = l/y/l-/3f = 34.932. 
This leads to a Cherenkov threshold of pt = m 7T / ~f t l3t = 4.87 GeV/c for charged 
pions (m n = 139.570 MeV/c 2 ), while electrons (m e = 0.511 MeV/c 2 ) exceed 
the Cherenkov threshold already with a momentum of pt = 0.018 GeV/c. 
Below the pion threshold the RICH has a hadron rejection of 10 4 to 1. 

A cut-away view of the RICH detector is shown in Fig. 2.10. The Cherenkov 
light is focused by two intersecting spherical mirrors with a total area of 20 m 2 
onto two arrays of 1280 photo-multiplier tubes (PMT) each which are located 
on either side of the entrance window. An average number of 10 photons per 
(3 ~ 1 particle are emitted under the angle of &c ~ 9 mrad. They are focused 
to a ring on the PMT array with an asymptotic radius of ~ 11.8 cm. The 
glass in front of the photo tube absorbs light with wave lengths below 200 nm. 
The mirror reflectivity is 83% at this wave length and rises to 90% at 250 nm. 

In p + p collisions the RICH also serves as Level-1 trigger on rare events 
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Figure 2.10: A cut-away view of one arm of the PHENIX RICH detector [90]. 

with electrons. The trigger is comprised of 64 non overlapping trigger tiles in 
each PMT array. Each trigger tile consists of 4(0) x 5(z) PMTs, an area which 
approximately corresponds to the size of a Cherenkov ring of a f3 ~ 1 particle. 

2.2.5 Electromagnetic Calorimeter 

High-energy electrons and photons lose energy in matter predominantly via 
Bremsstrahlung and e + e~ pair production, respectively. The amount of energy 
they lose is defined by the radiation length Xq, which is both (a) the mean 
length of traversed matter after which an electron has lost all but 1/e of its 
energy and (b) 7/9 of the mean free path for e + e~ pair production by a photon. 

The Electromagnetic Calorimeter (EMCal) measures the energy and posi- 
tion photons and electrons. Furthermore, it serves as trigger on rare events 
with high momentum photons. It comprises of eight sectors covering each 22.5° 
in azimuth and Ar] = ±0.35 in pseudo-rapidity. All four sectors of the West 
Arm and the two top sectors in the East Arm are shashlik type lead-scintillator 
(PbSc) sampling calorimeter. The two bottom sectors are lead-glass (PbGl) 
Cherenkov calorimeters, which had been used previously in the CERN exper- 
iment WA98 at the SPS. 

The PbSc calorimeter contains a total of 15,552 individual towers which 
are made of 66 sampling cells with alternating layers of 1.5 mm Pb and 4 mm 
scintillator (1.5%PT/0.01%POPOP) [98]. A module as shown in Fig. 2.11 
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Front 




Figure 2.11: View of a PbSc module showing the layers of Pb and scintillator, the 
wavelength shifting fibers and the phototubes attached to the back. 



comprises of four optically isolated towers which are read out individually. 
Each tower has measures 5.535 x 5.535 cm 2 across and has a length of 37.5 
cm, which corresponds to 18 Xq. 36 modules are held by a common support 
structure called super module. 18 super modules form a sector. The energy 
resolution of the PbSc Calorimeter is 

^ = ^2© 2.1%. (2.2) 



Each sector of the PbGl calorimeter comprises of 192 super modules (SM) 
which contain each 24 modules as shown in Fig. 2.12. Each module measures 
4x4 cm 2 across, is 40 cm long (14.3 X ), and read out with a photomultiplier 
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at its end. The energy resolution of the PbGl Calorimeter is 

^ = ^5 ©0.76%. (2.3) 

With a thickness of 18 Xq in the PbSc and 14.3 Xq in the PbGl, respec- 
tively, electrons and photons will deposit their energy within the calorimeter as 
electromagnetic shower of subsequent Bremsstrahlung and e + e~ pair creation. 
In the PbSc the electrons within the electromagnetic shower created in the 
Pb-layer produces scintillation light in the scintillator layers. The scintillation 
light is guided by wavelength shifting fibers to the phototubes located at the 
back of each tower. In the PbGl, which has a refractive index of n = 1.648, the 
electromagnetic showers is detected by Cherenkov light radiated by electrons 
in the shower. The Cherenkov light is read out at the end of the Calorimeter 
by photomultiplier tubes. 

In contrast to electrons and photons, the energy loss of hadrons in matter 
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Figure 2.13: Energy spectrum measured with the PbSc calorimeter, when exposed 
to pions, protons and electrons for incident energies of 0.5, 1, and 2 GeV [98]. The 
y-axis shows counts in arbitrary units. 



occurs primarily through ionization and atomic excitations. For typical hadron 
energies (0.1 < E < 10 GeV) the energy deposited in matter is nearly inde- 
pendent of the particle's energy, therefore these particles are called minimum 
ionizing particles, or mip's. Furthermore, PbSc has a nuclear interaction length 
Aj = 0.85 and PbGl A/ = 1.05, respectively. Therefore, only few hadrons will 
interact strongly and deposit a significant fraction of their energy. This leads 
to mip peaks in the energy spectrum as shown in Fig. 2.13 for charged pions 
and protons along with electrons for comparison. 

The EMCal serves as Level- 1 trigger for events with high momentum pho- 
tons, triggering when the energy deposited in an area of 4 x 4 overlapping 
towers surpasses a defined threshold. In addition the energy in a area of 2 x 2 
overlapping towers can be used in coincidence with the RICH trigger to trigger 
on events with electron candidates. 



2.2.6 Charged Particle Acceptance 

The central arm acceptance of charged tracks depends on their charge sign 
q, their transverse momentum and their azimuthal angle. The magnetic field 
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of the central magnets will bend a particle emitted under 0o and momentum 
Pt', its azimuthal angle measured at a distance r from the vertex increases 
inverse proportional to px- Therefore the relation between azimuthal angle at 
the vertex and at the DC (RICH) is given by: 

0dc = 0o + q— (2.4a) 
Pt 

, , . &RICH , n AU , 

0rich = 0o + q (2.4b) 

Pt 

An electron is accepted if its azimuthal angle is within the coverage of DC and 
RICH in one of the two central arms, i. e.: 

0min < 0DC — 0max 

&& 

min < 0RICH — 0max 

Fig. 2.14 shows the distribution of single electrons measured in p + p collisions 
in q/pr vs. 0o- The west arm coverage extends around 0^0 and the the 
east arm around ~ ir. The shapes can be described by the low px cut of 
0.2 GeV/c limits the distribution within \q/pr\ < 5 c/GeV shown as horizontal 
dashed black lines and the conditions defined in Eq. (2.5), with /cdc — 0.206 
rad GeV/c, fc RICH = 0.309 rad GeV/c, min = — ^7r rad and max = jq^ ra d 
for the west arm and m ; n = j^n rad, and max = y§tt rad for the east arm, 
respectively, whose boundaries are shown as solid black lines. 
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Figure 2.14: Single electron acceptance. The diagonal line represent the accep- 
tance limits due to RICH and DC as defined in Eq. (2.5). The dashed lines indicate 
the low p T cut off at 200 MeV/c. 
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Chapter 3 
Analysis 



In this chapter the analysis of the dielectron continuum in p + p as well as 
Au + Au collisions at a v /s/v r /v = 200GeV is presented. The analysis proce- 
dure of both data sets has large overlap so they are discussed in parallel, but 
emphasizing differences where they are present. 

The major steps of the analysis are outlined briefly here. The analysis 
begins with the event selection in Section 3.1, which is discussed for p + p in 
Section 3.1.1 and in Section 3.1.2 for Au + Au based on trigger and vertex in- 
formation. This is followed by a single electron analysis, including the tracking 
discussed in Section 3.2 and electron identification in Section 3.3. 

Single electrons are paired to electron-positron pairs as explained in Sec- 
tion 3.4. The subtraction of combinatorial and correlated background is dis- 
cussed in the Sections 3.4.3-3.4.5. The resulting invariant mass and pt spectra 
of e + e~ pairs are corrected for reconstruction and electron identification effi- 
ciencies as well as trigger efficiencies as presented in Section 3.5. The pt spec- 
tra of e + e~ pairs require an additional correction for the detector acceptance 
which is discussed in Section 3.5.4. 

Section 3.7 presents a calculation ( Cocktail) of the expected e + e~ pair yield 
from hadronic decays. The results will be compared to the experimental data 
in Chapter 4. 

3.1 Data Set and Event Selection 
3.1.1 p + p Collisions 

The analysis is based on the data set of p + p collisions at ^/s = 200 GeV 
which was collected during the run period in 2005. 

Two data samples are used for the measurement of the dielectron contin- 
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uum: a reference sample of events which were selected by the minimum bias 
trigger (MB) and a data set recorded with the single electron trigger (ERT). 
The MB trigger for p + p collisions requires that the BBC has registered hits 
in at least two photo tubes and an online determined collision vertex within 
30 cm: 

MB = (BBC > 2) n (|£ V c r tcx| < 30 cm) (3.1) 

The MB trigger cross section is obbc = 23.0 ± 2.2 mb corresponding to 
54.5 ± 6% of the inelastic p + p cross section a pp = 42.2 mb. Simulations, and 
data collected without requiring the BBC trigger, indicate that the triggered 
events include 79 ± 2% of events with particles in the central arm acceptance. 
This number coincides with the fraction of non-diffractive events triggered by 
the BBC from which it is concluded that for non-diffractive collisions the BBC 
trigger can have only little bias towards events with particles produced in the 
central arms. 

Events with a collision vertex far from the origin have a higher chance to 
create particles that hit the nose cones of the central magnet. This creates 
additional particles, e.g., from photon conversions which are reconstructed 
in the central arm detectors that do not originate from the actual collision. 
The Vertex distribution of events is shown in Fig. 3.1a. While this distribution 
is centered around z V ertcx = cm and has a full width at half maximum 
of w 30 cm, the number of electrons in the central arm peaks strongly for 
l^vertcxl > 25 cm, as shown in Fig. 3.1b. Therefore, an vertex cut of z vcrtcx < 
25 cm, with the vertex position determined offline by the BBC with better 
accuracy, is applied to avoid such contamination. Although with ~ 4 kHz the 
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Figure 3.1: Vertex distribution of events (a) and of electrons (b) in p + p collisions 
at y/s = 200 GeV. 

rate at which PHENIX can record data exceeds the capability of any other 
RHIC experiments, it would not be useful to only record events based on 
the MB trigger if one wants to study rare events, containing, e. g., electrons. 
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Therefore, PHENIX has a variety of Level-1 triggers to select particularly 
interesting events, e. g., events with a high particle in the central arm or an 
electron candidate. They are called ERT (EMCal RICH triggers), as they are 
based on the energy measurement in the EMCal, which in case of the electron 
trigger is matched to hits in the RICH. 

The high p? photon trigger is based on the energy measured within a 4 x 4 
neighboring EMCal towers. If this energy exceeds a set threshold, the trigger 
is fired. During the p + p run in 2005 there were three triggers with different 
energy thresholds: 

• ERTLLl_4x4a: E > 2.1 GeV 

• ERTLLl_4x4b: E > 2.8 GeV 

• ERTLLl_4x4c: E > 1.4 GeV 

In addition to the high pt photon trigger an electron trigger (ERTLL1_E) 
was active. This trigger requires a minimum deposited energy of 400 MeV in 
an overlapping tile of 2x2 EMCal towers matched to a hit in the RICH. The 
RICH hit is required within a trigger tile of 4x5 photo tubes. The location 
of the RICH trigger depends on the momentum of the trigger particle and 
is determined from a look-up table, assuming an electron to be the trigger 
particle, i. e. the momentum being equal to the energy deposited within the 
2x2 EMCal towers. Only events which are triggered in coincidence with the 
MB trigger (ERTLL1_E&BBCLL1) are considered for the analysis so that a 
cross section for MB collisions can be extracted. The efficiency of the ERT 
trigger is discussed in Section 3.5.2. When in the following the term "ERT 
trigger" is used, it refers to the ERTLL1_E&BBCLL1 trigger unless stated 
otherwise. 

The downscale factors on the MB trigger changes frequently to adjust for 
changing beam conditions. These downscale factors are recorded in a database 
and have to be considered when determining the total luminosity recorded. 

This analysis uses a MB data sample of 1.5 billion events and an ERT 
triggered sample of 270 million events. As the ERT trigger requires the coin- 
cidence of the MB trigger one can calculate the number of sampled MB events 
that the ERT data set of a given run corresponds to by as: 

^sampled _ jy/"^ . / scale _ down • ^jl^/A^g 1 " (3-2) 

Where Nmb is the number of events recorded with the MB trigger in this run 
and /scale-down is the scale-down factor for the MB trigger. The correction fac- 
tor N^rp/N^Q - is necessary in occasions when during the data reconstruction 
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a file segment of either the MB sample or the ERT sample is lost. In this case 
the number of ERT triggered events in the MB sample (A^b T ) does not equal 
the number of MB triggered events in the ERT sample (A^j^). This ratio is 
shown in Fig. 3.2a for all run numbers. Runs for which this ratio is > 2 or 
< 0.5 are rejected (5 runs), all others which are not equal to 1 are corrected 
(2 runs). 

The total number of MB events sampled by the ERT trigger is the sum 
of Eq. (3.2) over all runs: 

^sampled = g^gggg x jqIO 



3.1.2 Au + Au Collisions 

The analysis of the dielectron continuum in Au + Au at ^/snn = 200 GeV is 
performed on the data set recorded during the RHIC running time in 2004. 
Collisions were selected based on a minimum bias trigger, requiring at least 
two hits in each of the BBC and a coincident hit in one of the ZDCs. In 
addition the collision vertex had to be within 38 cm: 

MB = (BBCN > 2)n(BBCS > 2)n(ZDCS > 1UZDCN > l)n(|^ vcrtcx | < 38 cm) 

(3.3) 

92^ 3 ; % of the inelastic Au + Au cross section are selected by this trigger. For 
the same reason as in p + p a vertex cut is applied at z vcr tcx < 25 cm in the 
offline analysis. 

In contrast to protons, gold ions cannot be treated as point-like particles. 
They consist of 197 nucleons, 79 protons an 119 neutrons. When accelerated to 
an energy of 200 GeV per nucleon they are Lorentz contracted in longitudinal 
directions and look much like pan cakes with a radius of ~ 7 fm in the lab 
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Figure 3.3: Overlap of two gold nuclei for different impact parameters b. 



frame. Not every collision of two gold ions is head on, but their centers are 
shifted by a distance called the impact parameter b as illustrated in Fig. 3.3. 

The number of nucleons participating in a collision N paxt , depends on the 
impact parameter. Every nucleon that participates in a collision can undergo 
a number of binary collisions N co \\ with other nucleons. Both numbers can be 
calculated with a Glauber Monte Carlo simulation, which is a simple geometric 
model of the nuclei in which the nucleons are distributed following a Wood- 
Saxon potential. For the centrality classes used in this analysis the results of 
such a Glauber Monte Carlo simulations [99] are tabulated in Tab. 3.1. 
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Table 3.1: Average values of -/V par t, N co \\, and b for different centrality classes of 
Au + Au and p+p collisions at y/s~NN = 200 GeV [99]. Quoted errors are systematic 
uncertainties. Also listed are the number of events and the number of signal e + e _ 
pairs for each centrality class. 
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3.2 Charged Particle Tracking 

With the help of the Drift and Pad Chambers charged particles are tracked 
down to momenta of 200 MeV/c with a resolution of 

a p /p = 0.7% ®\%p/{GeV/c). (3.4) 

A track is reconstructed in the r — <p plane with six measurements in the XI 
and another six the X2 plane of DC over a radial distance of 20 cm. The drift 
time after which a hit is registered by the anode wire is measured with respect 
to the collision time to which is measured by the BBC. The working gas is 
chosen such that the drift velocity is uniform in the active region. This allows 
to convert the time after which a hit is measured by a particular wire into a 
distance at which the track passed this wire: 

x{t) = Vdrift • *■ (3-5) 

The drift velocity is obtained by measuring the time distribution of all hits in 
a DC arm. A typical distribution is shown in Fig. 3.4. The distribution has a 
characteristic peak which is about twice as high as the following plateau. This 
is due to tracks within a distance of 2 mm from the anode wire, the region 
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in which the back wires do not shield the field and the left-right ambiguity 
remains and therefore twice as many hits are registered as outside that region. 
The time interval between leading edge and trailing edge corresponds to a 
drift distance of 2 cm, which is the separation between the cathode and anode 
wires. The drift velocity is calibrated based on these time distributions for 
every run. The expected value for the DC working gas is: Udrift = 50 /mi/ns. 
Further calibrations are necessary to consider effects due to the signal width 
dependence of the arrival time, edge effects in the drift field, geometrical shifts 
between wires and the global alignment with respect to the collision vertex. 

Once the hit times are converted into distances a track candidate is re- 
constructed with a combinatorial Hough transformation over all possible hit 
combinations [100]. This is done for in a 2-dimensional space of the coordi- 
nates and a. The angle is the azimuthal angle of a track candidate at the 
DC reference radius of 223 cm, as shown in Fig. 3.5. The inclination a of the 
track candidate with respect to an infinite momentum (i. e., straight) track at 
angle <fi is inverse proportional to the momentum and its sign depends on the 
charge of the particle. Any combination of hits that have a local maximum 
and surpass the threshold criterion are considered a track candidate. After 
this Hough transformation further steps are performed to remove background 
tracks. First, hits are tested on their association to a track. A fit is performed 
that weights hits according to their distance from the track, such that the 
weight goes to zero for hits far away from the track. Then it is required that 
each hit can only belong to one track. Any track with less than eight hits is 
rejected. 

In addition the hits are measured in the UV wires as well as PCI hits 
are matched to track candidates to reconstruct the polar angle of the track 
candidate. If there are more than one PCI hit that can be associated to a 
track, the hit associated with the most UV hits is used. The following bits are 
used to determine the track quality: 

(1) XI used 

1 (2) X2 used 

2 (4) UV found 

3 (8) UV unique 

4 (16) PCI found 

5 (32) PCI unique 

Valid patterns include: 

49,50,51 1 1 x x PCI found/unique, no UV 

61,62,63 1 1 1 1 x x PCI found/unique, UV found/unique 
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17,18,19 1 x x PCI found/ ambiguous, no UVs 

21,22,23 1 1 x x PCI f ound/ ambiguous , UV found but tied 

29,30,31 1 1 1 x x PCI f ound/ ambiguous , UV found w/ one best choice 

The highest quality tracks can have is 63, i. e., it is reconstructed based 
on hits in the XI and X2 planes, has a unique PCI and UV hit. In addition 
to these, tracks with quality 31, i. e., requiring hits in XI and X2 plane and 
a unique UV hit, but only a ambiguous PCI hit, as well as with quality 51, 
%. e., demanding hits in XI and X2 plane and a unique PCI hit, but without 
a matching UV hit, are considered in this analysis. 

The measurement of a allows to reconstruct the transverse momentum of 
charged particles. There is a close relationship between a and the magnetic 
field integral along the trajectory of a particle: 



where K = 206 mrad GeV/c is the total magnetic field integral. With the 
measurements of a and <fi at the DC reference radius, one can determine the 
initial azimuthal angle at the collision vertex <p = <p — a + A0, with A</> = 



K 



(3.6) 



a ~ — 
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Figure 3.4: Example of a time distribution of hits in the Drift Chamber West arm 
summed over all wires (collected over ~ 15,000 p + p events). The x-axis displays 
time in units of ns. 




DC West Arm 



Figure 3.5: Definition of the coordinates <f> and a used in the Drift Chamber track 
reconstruction. 
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Figure 3.6: Ring mask on the RICH PMT array. -P croS s is the coordinate of the 
track projection around which the ring mask (shaded) is applied. The ring mask is 
limited by i2 min < r < # max 



3.3 Electron Identification 



Electrons emit an average number of 10 Cherenkov photons in the RHIC under 
an angle of cos(8) = l/(nf3), which are focused to rings in the PMT plane with 
an asymptotic radius of ps 5.4 cm. 

Based on the tracking information from Drift and Pad Chambers, a projec- 
tion of the track onto the PMT plane is calculated and stored in coordinates 
of the PMT plane: P C ross = ( Across, Across)- As illustrated in Fig. 3.6, around 
this point the number of PMT hits n with a distance between -R m i n = 3.3 cm 
and Rraax. — 8-4 cm are counted. Based on the pulse height measured in a 
PMT the number of photo electrons (npe i ) is determined. The ring center 
is reconstructed with no, the location Ri of the hit PMT, and the number of 
photo-electrons in this area (npe = ^ npej as: 

-'•'center — v-^ V 0, ' I 

22i n P e * 

with coordinates -Reenter = (-^center, ^center)- The distance to the track projection 
(displacement) is calculated as: 

disp = \/( Across ^center) 2 (0cross ^center) 2 - (3-8) 

In addition the quality of the measured ring shape is expressed in terms of the 
difference to the expected ring shape: 
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with R = 5.9 cm. 

The distance between the track projection onto the EMCal and the position 
of the associated EMCal cluster is calculated in azimuthal and z direction, 
emcdphLe and emcdz_e, respectively. These matching variables in units of 
radians and cm are converted into units of one standard deviation a by fitting 
the distribution to a Gaussian and stored as emcsdphLe and emcsdz_e. As the 
cluster position in the EMCal is particle species dependent due to the different 
shower shapes, this calibration is done for electron candidates. A circular 5 a 
cut on the distance is applied in p + p while in Au + Au a 3 a matching is 
required. 

In comparison to their momentum of p > 200 MeV/c, electrons have a 
negligible mass (m e = 511 keV/c 2 ). Therefore the energy they deposit in 
the EMCal must match the momentum (E = ^ p 1 + m 2 — p). In contrast, 
hadrons only deposit a fraction of their energy in the EMCal which leads to a 
measured energies which are smaller than their momenta. 

Fig. 3.7a shows the E/p distribution of charged tracks in min. bias p + p 
collisions in comparison with electron candidates, i. e., charged tracks fulfilling 
all but the E/p elD cut. While the distribution of all charged tracks shows no 
clear peak due to electrons at E/p = 1, requiring the elD cut greatly improves 
the signal to background ratio. The width of the peak of « 14% is dominated 
by the energy resolution of the EMCal given in Eqs. (2.2, 2.3). 

The reconstructed energy does not match the momentum for all electrons, 
as in some cases their shower overlaps with a photon shower which leads to 
a larger energy and causes the tail at E/p > 1. Electrons from off- vertex 
decays or late conversions have a misreconstructed momentum, as the tracking 
algorithm assumes all tracks to originate from the collision vertex. Off-vertex 
decays traverse less magnetic field integral and are therefore bent less, which 
leads to a larger reconstructed momentum and E/p < 1. In the momentum 
range of p < 4 GeV/c, relevant for this analysis, the momentum resolution of 
the drift chamber Eq. (3.4) is better than the energy resolution of the EMCal 
Eqs. (2.2, 2.3). Therefore, in the e + e~ pair analysis the invariant mass and pt 
of the e + e~ pair is calculated using the momentum information rather than 
the energy Eq. (3.10). As a consequence, electrons with a mismeasured energy 
can be kept, while electrons with misreconstructed momenta are removed with 
a cut on E/p > 0.5 in p + p. In Au + Au the momentum dependence of the 
E/p distribution is corrected and the matching is stored in units of a standard 
deviation dep = (E/p — l)/a E / p . This parameterization is determined by a 
Gaussian fit of the E/p distribution of electron candidates as function of px- 
A value of dep > —2 is required for tracks to pass the electron identification 
cut. 
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Figure 3.7: E/p distribution in min. biasp + p (a) and min. bias Au + Au (b) for 
charged tracks (black) , tracks after applying the RICH cuts (blue) and contribution 
of hadrons randomly associated to hits in the RICH (red). 



Table 3.2: Electron identification cuts for p + p and Au + Au. 



elD cut p + p Au + Au 

quality 63||31||51 

n > 1 > 2 

disp < 10 cm < 5 cm 

X 2 /npe < 15 cm 2 < 10 cm 2 

E/p > 0.5 

dep > —2 a 

y/ emcsdphLe 2 + emcsdz_e 2 < 5 a < 3 a 



In addition, some background remains underneath the peak which is at- 
tributed to random associations of RICH hits to hadrons. This contribution 
can be estimated by matching tracks in the north side of the drift chamber to 
hits in the south side of the RICH and vice versa. This is done by swapping 
the RICH sides during the track reconstruction. In contrast to a single elec- 
tron analysis, this hadron contamination is unreducable in the pair analysis. 
Fig. 3.7b shows the same distributions for min. bias Au + Au collisions, in 
which the hadron contamination is much larger. In the region 0.7 < E/p < 1.2 
the contribution to the sample of electron candidates is about 1.6% in p + p 
and 16% in Au + Au collisions All electron identification cuts are summarized 
in Tab. 3.2. 
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3.4 Pair Analysis 



In an event the source of any electron or positron is unknown and therefore 
all electrons and positrons are combined to foreground pairs, like sign FG ++ , 

FG and unlike sign FG + _. In the p + p data set it is checked that at 

least one of the tracks in each pair is associated to a hit of the ERT trigger. 
From the single track information, the invariant mass m ee and the transverse 
momentum px of the e + e~ pair are calculated as: 

m le = {.V++V-) 2 

= {E + + E.f-{p + +p.f (3.10) 

PT = (Px,+ +Px,-) 2 + (Py,+ + Py-f (3.H) 

with E± = \Jp± + ml , m e = 511 keV/c 2 , and the 3-momentum vector p± as 
measured with the drift chamber. 

The foreground pairs can be generally separated into physical and unphys- 
ical pairs. With physical pairs it is referred to an e + e~ pair which originates 
from the same parent, which can be either a hadron, e.g., u — > e + e~ or a 
photon. The photon can either be virtual and convert internally (such a vir- 
tual photon appears for instance in the it Dalitz decay: tt° — > 77* — > 7e + e~), 
or real and convert in detector material [e.g., the two photon decay of a tt°: 
7t° — > 77 — > 7e + e~). The latter type of photon conversions, also referred to as 
external conversion, are for reasons that will be explained in Section 3.4.2 not 
reconstructed with their real mass m = but contaminate the invariant mass 
spectrum of e + e~ pairs up to m ee ps 0.3 GeV and need to be removed from 
the sample of physical pairs in order to extract the physics signal of hadron 
decays and internal photon conversions. By definition, physical pairs can only 
exist as unlike sign pairs. However, the chance to reconstruct both the electron 
and the positron of such pairs is very much reduced by the limited azimuthal 
acceptance and the low momentum cut off at p < 200 MeV/c. 

The other, much larger contribution to the foreground, are unphysical pairs, 
i. e., electrons and positrons that are not from the same hadron, but just a 
result of combining all electrons in the event. Unphysical pairs can either be 
of unlike- or like-sign charge combinations. Most of these pairs are completely 
uncorrelated combinations of electrons, which can be described with a com- 
binatorial background generated with an event mixing technique. However, 
some fraction of the unphysical pairs is correlated, either due to correlations 
early on in the decay history, i. e., the pair does not share the same parent, but 
the same "grand parent" particle, in which case they are referred to as cross 
pairs or because they were produced within the same jet, therefore denoted jet 
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pairs. Another source of correlated pairs are ghost pairs, which appear when 
tracks share the same detector information. 

This classification scheme has ignored one very important source of e + e~ 
pairs, i. e., e + e~ pairs from semi-leptonic open charm and bottom decays cor- 
related through flavor conservation, discussed in Section 1.2.2. They do not 
originate from the same parent, however, they are considered physical pairs. 

In this Section all the details of the various background subtractions are 
discussed. Beforehand, a schematic overview over the various background 
sources is given as orientation. Noted in parenthesis are the Section number 
in which the particular background is discussed and a keyword referring to the 
method with which this background is removed. 

All Pairs: 

• physical pairs 

— hadron and open charm decays, internal conversions (Signal) 

— external photon conversions (4>v cut, Section 3.4.2) 

• unphysical pairs 

— uncorrelated (Mixed Events, Section 3.4.3) 

— correlated 

o Detector Ghosts (Pair Cuts, Section 3.4.1) 

o Jet and Cross Pairs (Simulations, Section 3.4.4) 

3.4.1 Pair Cuts 

In order to reproduce the shape of the combinatorial background of uncorre- 
lated e + e~ pairs with mixed events, it needs to be assured that any detector 
correlations are removed from the sample of e + e~ pairs in real events. Such 
correlations can occur, when two electron candidates share detector informa- 
tion, e. g. have overlapping rings in the RICH, share the same hit in PCI 
or have overlapping cluster in the EMCal. For instance, two electrons which 
traverse the RICH along parallel trajectories would emit Cherenkov photons 
which are focused on the same ring, as any two parallel photons are focused 
by the spherical mirror onto the same focal point. Therefore, a RICH hit due 
to an electron would be matched also to any charged track crossing the RICH 
parallel to the electron and therefore fulfill all elD cuts on RICH variables, 
even if the charged track itself does not emit Cherenkov photons. This effect 
is illustrated in Fig. 3.8. 
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Figure 3.8: Overlapping rings in the RICH PMT plane due to two parallel tracks. 



While an uncorrelated hadron contamination would be reproduced by mixed 
event, such a ghost pair is highly correlated in its geometry and therefore in 
the invariant mass spectrum. They have typical small opening angle and are 
therefore reconstructed with a small invariant mass. As like-sign tracks are 
bent by the magnetic field bents in the same direction, in contrast to unlike- 
sign pars which are bent in opposite directions, they have even smaller masses 
than unlike-sign pairs. 

Naturally, these ghost pairs do not exist in mixed events and therefore 
have to be rejected in real events. They are rejected with an event cut which 
rejects events (real and mixed) if there are two tracks which overlap in RICH, 
EMCal, or PCI. On overlaps in the RICH a cut is applied based on the 
distance between the ring centers associated to the two electron candidates of 
a pair. The distance between two tracks % and j is calculated as: 



^center — \l (^center ^center) 2 /Renter (^center Reenter) 2 /°" 2 > ccntor (3.12) 

with o~ ZccntcI = 3.6 cm and cr<£ centor = 10 mrad being the rms of the distributions 
of ring distances (without any pair cut) in z and directions, respectively, 
as shown in Fig. 3.9. A cut of A ccntcr > 10 (equivalent to 36 cm) is applied, 
which corresponds to twice the expected ring diameter (~16.8 cm). One may 
consider to cut only on one ring diameter as this would be enough to remove 
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Figure 3.9: Distance of the ring centers associated to two tracks in A(p centcr (a) 

and Az C enter (t>). 



false positive ring association to charged hadrons, but in addition one observes 
a reduction in the reconstruction efficiency of two electrons which have close 
rings, as their rings may overlap and distort their shape, which increases, e. g., 
X 2 /npe . This effect only vanishes once the two ring are clearly separated 
by more than one ring diameter, i. e. their centers by more than two ring 
diameters. In the EMCal events are rejected if two electron candidates are 
associated to two clusters within a 3 x 3 tower region (which is the average size 
of an electromagnetic shower). In the PCI the cut is applied at Az < 0.5 cm 
and A<f> < 20 mrad, which is the size of a cell. The pair cut on RICH variables 
has the largest effect and the combined cuts remove ~ 0.08% of all events in 
Au + Au collisions and a fraction of pairs which varies from 4% in the most 
central to 2% in the most peripheral collisions. In p + p this cut removes only 
0.1% of all pairs. 

The effect of ghost pairs on the e + e~ pair distribution can be studied 
with the invariant mass distribution of like-sign pairs, i. e., e + e + and e~e~ 
pairs, which do not contain any physics signal. The real and mixed event 
distributions of like-sign pairs in MB Au + Au collisions is shown in black and 
red, respectively, in Fig. 3.10a. The mixed events clearly deviate in shape from 
the real event distribution. In the invariant mass distribution of e + e~ pair the 
ghosts are also clearly visible at m ~ 500 MeV/c 2 as shown in Fig. 3.10b, 
where the mixed events deviate from the real events. 

3.4.2 Photon Conversions 

Photon conversions contribute to the final spectrum of e + e~ pairs and should 
be rejected. Electrons from conversions that are created off vertex in detector 
material, e.g., the beam pipe made of Beryllium (0.3 X ), pass through less 
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Figure 3.10: Invariant mass distribution of like-sign (a) and unlike-sign pairs(b). 
All pairs are shown in black, ghost pairs in blue and all others in red. 

magnetic field integral than the PHENIX tracking algorithm assumes, which 
leads to a misreconstructed momentum. In addition the opening angle of the 
e + e~ pair is misreconstructed, as it is tracked back to the collision vertex, 
rather than to the conversion vertex, as illustrated in Fig. 3.11. Therefore, 
conversion pairs have an apparent mass that increases with the distance from 
the collision vertex and therefore with the misreconstructed opening angle. 

Photon conversions occur in the beam pipe at a radius of 4 cm which 
corresponds to the peak at m ee « 20 MeV/c 2 in the invariant mass spectrum 
shown in Fig. 3.12, in detector support structures (m ee « 80MeV/c 2 and m ee ~ 
125 MeV/c 2 ) as well as in the air before, and much reduced within, the installed 
Helium bag which contributes as a continuum out to m ee < 300 MeV/c 2 , 
corresponding to the DC entrance window, beyond which electrons do not 
bend anymore because the region is field-free. Therefore, electrons from photon 
conversions within the DC are straight and rejected with the high-p<r cut of 
20 GeV/c. As photons are massless, conversion pairs do not have an intrinsic 
opening angle, but are aligned perpendicular to the axial magnetic field. A cut 
on the orientation of e + e~ pairs is used to remove those conversion photons. To 
measure the orientation of the e + e~ pair plane with respect to the z direction, 
i. e., the magnetic field, the angle (fiv is defined [101]: 



P1+P2 

U — T= =TT 

\Pl+P2\ 


(3.13a) 


V = P\ X f> 2 


(3.13b) 


W = U X V 


(3.13c) 


U X Z 


(3.13d) 


Ua= 


it x z\ 
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Figure 3.11: Misreconstruction of the opening of a conversion pair. Shown are 
the photon and the conversion pair which is created in the beam pipe. The real 
opening angle between the original trajectories (dashed lines) is much smaller than 
the reconstructed opening angle, artificially created by tracking the particles back 
to the collision vertex (dotted lines). 



Here pi is the 3-momentum vector of the electron and f>2 the 3-momentum 
vector of the positron. The following cuts are applied: 

• m< 0.60 kk <p v > 0.1 rad 

• m < 0.03 kk <p v > 0.25 rad 

To remove beam pipe conversions a harder cut has to be applied, because 
the multiple scattering within Beryllium smears the orientation with respect 
to the magnetic field more than in air. In addition the resolution of the <fiv 
measurement improves with increasing conversion radius. This cut removes 
> 98% of all conversion pairs. 



As mentioned in the introduction to this Section, most of the e + e~ pairs are of 
uncorrelated origin. While the unlike-sign spectrum contains both physics sig- 
nal and background pairs, the like-sign spectrum does not contain any physics 
signal and therefore is a direct measure of the background. It can be shown (see 




(3.14) 



3.4.3 Event Mixing 
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Figure 3.12: Invariant mass distribution of all e + e~ pairs in real events (black), 
photon conversions (red) and remaining pairs (blue). The peak at 20 MeV/c 2 is 
caused by photon conversions in the beam pipe. The remaining shape is due to 
other detector support structures (see text). 
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Appendix A) that, as long as electrons and positrons are produced in pairs, the 
unlike-sign combinatorial background is the geometrical mean of the like-sign 
backgrounds, independent of efficiency and acceptance differences, i. e.: 



(BG+-) = 2y/(BG ++ )(BG—}. (3.15) 

In experiments with equal acceptance for like- and unlike-sign pairs, the distri- 
bution of like-sign pairs in real events can be directly related to the unlike-sign 

background through their geometrical mean BG + - = 2^FG ++ FG and 

subtracted from the unlike-sign pair spectrum to obtain the physics signal. 
Because of a different acceptance for like- and unlike-sign pairs a correction 
needs to be applied to use this method in PHENIX. 

Another possibility, which leads to a higher statistical accuracy, is the event 
mixing. Electrons from events with similar topology are combined to pairs to 
statistically approximate the combinatorial background. To ensure the same 
topology, events are categorized in pools according to their centrality, z-vertex 
and reaction plane. In Au + Au, mixed events were generated in ten centrality, 
six z-vertex and one reaction plane pools 1 , while five z-vertex pools were used 
for p + p. 

The ERT trigger used in the p + p analysis biases the single electron dis- 
tribution towards higher pr- Therefore, to generate the correct combinatorial 
background shape of e + e~ pairs, the mixed events must be generated from the 
minimum bias sample, but as in real events it is required that in every pair at 
least one of the electrons must fulfill the ERT trigger condition. 

The shape of the like-sign background from mixed events can be compared 
to the like-sign distribution in real events to determine in which region and how 
well the background shape is reproduced in the event mixing. The absolute 
normalization of unlike-sign combinatorial background is given by the geomet- 
rical mean of the observed positive and negative like-sign pairs 2y / BG ++ BG 

where BG ++ and BG are determined by integrating the mixed event dis- 
tributions after they have been normalized in a region where the mixed events 
reproduce the distribution from real events, which means it is enforced that: 

/ ,, BG-L+ + BG — 
where A denotes the normalization region. 



^^One reaction plane pool means no separation by reaction plane. However, it was tested 
that the reaction plane pooling has no effect on the mixed event distribution and can be 
neglected 
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Figure 3.13: Mixed event subtracted like-sign spectrum in data (left) and as 
calculated with Pythia (right) as function of invariant mass and px- The dashed 
region outlines the normalization BjI Q'(\ 3jS defined in Eq. (3.17). 



Background Normalization 



The like-sign spectra as function of mass and px after mixed event subtrac- 
tion are shown in the left panel of Fig. 3.13. The figure clearly shows that 
mixed and real events do not have the same distribution like-sign pairs, which 
is an indication of correlated pairs. However, in the region m ee ~ px the 
distributions are very similar. Since for a pair with opening angle 5 it is 
m ee = Pi£>2(l — cos 5), the condition m ee ~ p T roughly corresponds to an 
azimuthal opening angle of A<fr ~ ir/2. This was checked with a Pythia 
simulation of minimum bias events which is described in Section 3.4.4. The 
A0 distribution of like-sign pairs from these simulated events is shown in 
Fig. 3.14. This distribution is compared to the mixed event distribution also 
generated from the Pythia simulation. The mixed events agree well with the 
real events around A</> ~ it/2. Thus, the mixed event distribution in Pythia 
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Figure 3.14: A(p distribution of Pythia like-sign pairs in same (black) and mixed 
(red) events. 



is normalized in the region defined by: 

m ee > 300 MeV/c 2 && (3.17a) 

m T < 1.2 MeV/c 2 && (3.17b) 

p T /c - 1.5m ee < 200 MeV/c 2 && (3.17c) 

p T /c - 0.75m ee > MeV/c 2 . (3.17d) 

which is shown as the dashed contour in Fig. 3.13. The yield remaining after 
the subtraction such normalized mixed event is considered correlated and dis- 
cussed in Section 3.4.4. The same region is used to normalize the mixed events 
to the real events in data. The normalization to this region has a statistical 
accuracy of 2.8%. 



3.4.4 Correlated Background 

Fig. 3.15 shows the invariant mass distribution of like-sign pair of real and 
mixed events in p + p collisions. The pairs from mixed events do not reproduce 
the shape of like-sign pairs in real events which is an indication for correlated 
background whose two sources are discussed in the following. 
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Figure 3.15: Invariant mass distribution of all like-sign pairs in real events (black) 
and mixed events (red). The normalization is discussed in Section 3.4.3. A clear 
shape difference is observed at small and at large masses, indicating a source of 
correlated background. 



SO 



Cross Pairs 



The first correlation of like-sign pairs stems from hadron decays with two e + e~ 
pairs in the final state, e. g., for a tt° these are the double Dalitz decay (tt° — > 
e + e~e + e~), the Dalitz decay (ir° —> ^e + e~) in which the photon converts into 
an e + e~ pair in detector material, and the two photon decay (tt — > 77) in 
which both photons convert into two e + e" pairs 2 : 

7T° -> 7l72 

-> e^eT/eJe^ (3.18) 

The e + e~ pair with the same parent photon is considered physics signal (i e., 
e^~e7 and e^e^~), and in case of a real photon conversion in detector material 
removed with a cut (see Section 3.4.2). The "cross" combinations into two 
unlike-sign pairs {efe^ and e^ej") as well as into two like-sign pairs {efe^ and 
e i~ e 2~) are n °t purely combinatorial, but correlated via the ir° mass (m ee < 
m n o). Therefore, their contribution is present in the real event spectra of 
both unlike- and like-sign pairs, but not in mixed events. Fig. 3.15 shows the 
invariant mass distribution of like sign pairs in real (black) and mixed events 
(red). The peak visible in the real event distribution of like-sign pairs at a 
mass of m ee « 100 MeV/c 2 is due to cross pairs. The same cross combinations 
are possible for rj decays, with the additional decay channel of r/ — > 7r°7r 7r°. 
Therefore these correlated pairs extend up to a mass of m n m 550 MeV/c 2 . 

For every decay, contributions to unlike-sign signal U, unlike-sign cross 
C and like-sign cross pairs L can be calculated in terms of branching ratios 
{BRyry = 98.8%, -B-Roaiitz = 1-2%, and -B-Rdbi.Daiitz = 

3.14 x 10" 5 [29]) and 
conversion probability (P con v = 0.3%) relative to the 7r° multiplicity: 

• 7T° — > 77 — > e + e~e + e~: 

- U = 

- C = 2- BR jy ■ P 2 onv = 2 ■ 98.8% • 0.3% 2 = 1.76 x 10~ 5 

- L = 2 ■ BRy 7 ■ P 2 onv = 2 ■ 98.8% • 0.3% 2 = 1.76 x 10~ 5 

• 7T° — > 77* — > e + e~e + e~: 

- U = BR Dalitz = 1.2% 

- C = 2 ■ BR Dahtz ■ P conv = 2 ■ 1.2% • 0.3% = 7.2 x 10~ 5 

- L = 2 ■ Bifeaiitz ■ ^onv = 2 ■ 1.2% ■ 0.3% = 7.2 x 10- 5 



2 All these decays are in principle two photon decays, but with two, one or zero virtual 
photons 
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which gives a total of U = 1.2 x 10~ 2 unlike-sign signal pairs, C — 1.5 x 10~ 4 
unlike-sign cross pairs, and L — 1.5 x 10 -4 like-sign cross pairs per 7r°. 

The same calculation for the t] yields a total of U = 6 x 10~ 3 unlike-sign 
signal pairs, C = 6.7 x 10~ 4 unlike-sign cross pairs, and L = 6.7 x 10~ 4 like- 
sign cross pairs per r/. Here the branching ratios (i?i? 77 = 39.38%, i?i?Daiitz = 
0.6% [29]) have been used and the unknown double Dalitz branching ratio of 
the has been assumed to be the square of the single Dalitz: BR^u. Dalitz = 
B -^Dalitz = 1-23 x 10~ 5 . The three ir° decay of the eta contribute with an 
additional C = L = 6.09 x 10 -4 to the unlike- an like-sign cross pairs. 

All these contribution are simulated with a fast Monte Carlo (Exodus) of 
7T° and t] double Dalitz decay as well as t] — > 7r°7r 7r° decays and filtered into 
the PHENIX acceptance. The resulting unlike- and like-sign cross pairs are 
shown Fig. 3.16. 

Jet Pairs 

Besides cross pairs another source contributes to correlated background pairs. 
This contribution can be reproduced with minimum bias events from Pythia 
simulations [102-104]. Pythia is a computer program widely used as Monte 
Carlo simulation of high-energy elementary particle collisions, e. g. e + p or 
p + p. The particle interactions are calculated using leading order matrix 
elements and the Lund string fragmentation model [105]. Pythia 6.319 with 
CTEQ5L parton distribution function [106] has been used. The following hard 
QCD processes have been activated: 
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Here g denotes a gluon and a fermion with flavor i, i. e., u, d, s, b, t, e~, u e , 
fi~ , i>n, t~, or u T , and their corresponding antiparticles Flavors appearing 
in the initial carry indices i or j, while flavors created only in the final state 
are denoted k. A Gaussian width of 1.5 GeV for the primordial kx distribution 
is used (PARP(91) = 1.5), and the minimum parton p? has been set to 2 GeV 
(CKIN(3)=2.0). 

In this simulation the branching ratio of the ir° Dalitz decay was set to 
100% to increase the sample of e + e~ pairs. The resulting A0 distribution of 
like sign pairs from real (black) and mixed events (red) in Pythia is shown in 
Fig. 3.14. The mixed event like-sign pairs have been normalized to the like-sign 
pairs in real events in the range 70° < A0 < 110°, in which the distributions 
of like-sign pairs in real and mixed events agree well. The deviation in shape 
at A(p « can be attributed to electrons which originate from two n° decays 
which occurred within the same jet, while the correlation at A(f) « tt comes 
from electrons that are decay products of 7r°'s in back-to-back jets. 7r°'s within 
the same jet are close in <p which leads to a small opening angle of the e + e~ 
pair, which also causes these pairs to have a rather small mass and at the 
same time a large pt as their momentum vectors add constructively. e + e~ 
pairs correlated via back-to-back jet have large opening angles and therefore 
larger invariant masses, while their p? is smaller than the one of pairs in the 
same jet, as their momentum vectors cancel mostly. These correlations in 
m ee and pt are shown in Fig. 3.13 for like-sign pairs from Pythia after the 
subtraction of mixed events. 

Fig. 3.16 shows invariant mass spectra of like sign pairs (top) and unlike- 
sign pairs (bottom) in real events, the mixed event distribution normalized 
according to Eq. (3.17) and the distributions after mixed event subtraction. 
The remaining like-sign correlated pairs are fitted with the distribution of 
cross and jet pairs in the two mass regions 0.06 < m ee < 0.12 GeV/c 2 and 
0.80 < m ee < 2.00 GeV/c 2 : 

0.12 0.12 0.12 

— dm = A / — dm + B / — dm 3.19a 

WW like J \ dm J cross J WW jet 

0.06 0.06 0.06 J 

= *, (3,9b) 

jet 



V^/iike J \dmj 

0.8 0.8 



The resulting normalization factors are applied to the unlike-sign distributions 
of cross and jet pairs which are shown in Fig. 3.16. The unlike-sign signal 
results from subtracting the correlated background from all correlated pairs. 
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Figure 3.16: Raw dielectron spectra. The top panel shows like-sign pairs as mea- 
sured in the experiment, the combinatorial background from mixed events, the cor- 
related pair background obtained by subtracting the combinatorial background, and 
the individual contributions from cross and jet pairs to the correlated background 
(see text). The bottom panel shows the same distributions for unlike-sign pairs. 
The correlated background in both panels is normalized to the measured like-sign 
pairs remaining after subtracting the combinatorial background. 
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3.4.5 Like-sign Subtraction 

As a systematic check of the subtraction of combinatorial and correlated back- 
ground an independent, completely data driven method has been developed, 
that does not make any assumption about the background contributions. As 
stated at the beginning of the chapter, the like-sign spectrum is free of physics 
signal and therefore a measure of the background, whether correlated or uncor- 
rected. But the different acceptance of unlike- and like-sign pairs in PHENIX, 
needs to be considered before they can be subtracted from the real event spec- 
trum of e + e~ pairs. 

The acceptance difference between unlike- and like-sign pairs is preserved 

in mixed events, therefore the ratio of BG + _/(BG ++ + BG ) is the proper 

correction to convert the distribution of like-sign pairs in real events into the 
background distribution of e + e~ pairs. The signal of e + e~ pairs is calculated 
in m ee and px as follows: 



In the special case that the real event like-sign spectrum is perfectly described 
by the mixed events, this expression reduces to S — FG + - — BG + -. 

The background distributions of like- and unlike-sign pairs are shown in 
Fig. 3.4.5 together with the relative acceptance. 

Alternatively to using the like sign foreground, one can first subtract the 
mixed events (normalized as discussed in Section 3.4.3) and then subtract 
the mixed event subtracted like sign signal, corrected for the relative accep- 
tance, from the mixed event subtracted unlike-sign spectrum everywhere but 
in the normalization region 3 . To increase statistics at high p? and large in- 
variant mass, the Au + Au mixed events have been used as an alternative 
description of the relative acceptance. The relative acceptance for correlated 
and uncorrelated pairs may differ. This effect has been studied its effect by 
using Pythia foreground spectra to determine the relative acceptance, i. e., 

FG + ^/(FG ++ + FG ). In addition the mixed events from Pythia have 

been used to describe the relative acceptance. Fig. 3.17 shows a comparison 
of the subtracted unlike-sign spectra with the presented methods. 

There are two systematic uncertainties in the background subtraction, one 
due to the statistics available in the normalization region for the various back- 
ground sources: 

• mixed events: 2.43% 

3 Note that a subtraction including this region would be arithmetically identical to the 
previous method, which skips the mixed event subtraction. 




(3.20) 
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m e . e [GeV/c 2 ] 



(c) Relative Acceptance 

• cross pairs: 2.49% 

• jet pairs: 6.86% 

The systematic uncertainty on the measured signal is shown in Fig. 3.18. 
The other systematic uncertainty is in the shape of the subtracted correlated 
background, which is estimated from the RMS of the various subtraction meth- 
ods described above, which is shown in Fig. 3.19. For m ee < 2 GeV a system- 
atic uncertainty of 5% (10% above) has been assigned. 

3.4.6 Background Subtraction in Au + Au 

The like-sign distribution of real and mixed events in minimum bias Au + Au 
collisions are compared in Fig. 3.20. Their shapes agree not only in the region 
defined in Eq. (3.17), but for all masses above 550 MeV/c 2 . Qualitatively, this 
can be explained with the suppression of away-side jets observed in Au + Au 
collisions [107]. At small masses, below 200 MeV/c 2 a like-signal from cross 
pairs is observed. The ratio of mixed event subtracted to mixed events is fitted 
above 550 MeV/c 2 to a constant, which gives a result of -0.000259 ±0.000633 
well in agreement with zero and a x 2 7 1 DOF = 1.45. Fig. 3.21 shows the like- 
sign mass spectra of real and mixed events divided in centrality bins of 0-10%, 
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Figure 3.17: Comparison of background subtraction methods. 




Figure 3.18: Comparison of background subtraction methods. 
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(a) LMR (b) IMR 

Figure 3.19: RMS of all methods as function of m, 



Table 3.3: Fit parameters to ratios shown in Fig. 3.20 and Fig. 3.21. 



Centrality 


Constant 


X 2 /DOF 


X 2 test 


max. dist 


0-10% 


0.00063 ± 0.00088 


1.25 


0.58 


0.0013 


10-20% 


-0.0009 ± 0.0014 


1.42 


0.54 


0.0018 


20-40% 


-0.0024 ± 0.0018 


1.12 


0.56 


0.0037 


40-60% 


-0.0085 ± 0.0050 


1.42 


0.69 


0.0099 


60-92% 


-0.018 ± 0.016 


1.56 


1.02 


0.042 


00-92% 


-0.00026 ± 0.00063 


1.45 


0.51 


0.0011 



10-20%, 20-40%, 40-60% and 60-90%. The same fits as for the MB spectra 
are performed and summarized in Tab. 3.3, which also includes the results 
of a statistical x 2 test and the maximum deviation of a Kolmogorov-Smirnov 
test (also see inserts to Fig. 3.20). The maximum deviations of Kolmogorov- 
Smirnov test are small compared to the uncertainties of the absolute normal- 
ization of the mixed event background. The p-values corresponding the the 
X 2 test, which are greater than 0.999 for all centrality classes, confirm the 
hypothesis that the two distributions are compatible and can be accepted for 
any commonly used significance level. 

Finally, Fig. 3.22 shows the invariant mass spectra of like-sign pairs in real 
and mixed event for different px ranges: p? < 1 GeV/c, 1 < Pt < 2GeV/c, 
and px > 2 GeV/c. Again no deviation of the mixed event shape from the 
real event distribution is observed for m ee > 550 MeV/c 2 . 

The mixed event are normalized to the like-sign pairs in real events in the 
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Figure 3.20: Invariant mass distribution of like sign pairs in real events (open 
circles) and mixed events (red line). The mixed events to real events are normalized 
above 700 MeV/c 2 . The bottom panel shows the ratio of the mixed event subtracted 
spectrum and the mixed events. 
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Figure 3.21: Invariant mass distribution of like sign pairs in real events (open 
circles) and mixed events (red line) for different centrality bins. The mixed events 
to real events are normalized above 700 MeV/c 2 . The bottom panel shows the ratio 
of the mixed event subtracted spectrum and the mixed events. 
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Figure 3.22: Invariant mass distribution of like sign pairs in real events (open 
circles) and mixed events (red line) for different px bins. The mixed events to real 
events are normalized above 700 MeV/c 2 . The bottom panel shows the ratio of the 
mixed event subtracted spectrum and the mixed events. 
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mass range 0.7 < m ee < 4 GeV/c 2 : 

fo.7 BG++ + BG — 



Ioj + FG. 



(3.21) 



The lower limit of 0.7 GeV/c 2 is chosen conservatively to avoid any contami- 
nation from correlated background observed at low mass and provides at the 
same time sufficient statistical accuracy of 0.12%. The absolute normalization 
of the unlike-sign background is given by the geometrical mean of the nor- 
malized like-sign integrals: BG^ = ^BG ++ + BG An Additional 0.2% 

uncertainty, due to a correction of that normalization by 1.004 ±0.002 for the 
fact that the applied pair cuts removes a different fraction of like- than unlike- 
sign pairs, yields to a total systematic uncertainty of 0.25% on the background 
normalization, which is a conservative upper limit of the uncertainty. 

The correlated background, remaining after mixed event subtraction, due 
to cross pairs and a contribution of jet-pairs at small masses is subtracted 
analog to the p + p analysis. But, in contrast to p + p, there is no room 
for a correlated signal due to away-side jet correlations. Due to the known 
modifications of the away-side jet, the Au + Au data are indeed expected 
to differ from the Pythia simulation. Therefore, the near-side correlations, 
i.e., pairs with A0 < 7r/2 are separated from the away-side contribution with 
A0 > 7r/2 and the three components cross, near-, and aways-side pairs are 
fitted to the like-sign correlated background distribution: 

\<Wlikc \ dm J cross WW near V^Vaway 

As the mixed events have been normalized at m ee > 0.7 GeV/c 2 it follows 
that C = 0. The different background contributions are illustrated in Fig. 3.23. 
Here the away-side curve is shown for illustration only and is normalized as 
the near-side curve. The overall contribution of correlated pairs relative to the 
unlike-sign signal is small, and comparing different subtraction strategies leads 
to a systematic uncertainty of 10% on correlated background subtraction in 
the LMR. 
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Figure 3.23: Raw dielectron spectra in min. bias Au + Au. The top panel shows 
like-sign pairs as measured in the experiment, the combinatorial background from 
mixed events, the correlated pair background obtained by subtracting the combina- 
torial background, and the individual contributions from cross and near- and away- 
side jet pairs to the correlated background (see text). The bottom panel shows the 
same distributions for unlike-sign pairs. The correlated background in both pan- 
els is normalized to the measured like-sign pairs remaining after subtracting the 
combinatorial background. 
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3.5 Efficiency Correction 



The distribution of e + e~ pairs after subtraction of all background sources needs 
to be corrected for the pair efficiencies of track reconstruction (ej^), electron 
identification (s^d) an d (i n case of p + p) triggering by the ERT (^ert)- The 
pair efficiencies for track reconstruction and electron identification are the 
product of the single efficiencies, i. e., e££ = e+ co e- co and = e^e"^, 
while the trigger efficiency for a pair can be expressed in terms of the single 
trigger efficiency as: £ ERT = (1 — ^ertO^ert + (1 — £ ert) £: ert + £ ert £ ert- 
Tracking and electron identification efficiency include effects due to dead areas 
within the detector acceptance. In addition the efficiency of the pair cuts 
(£pair cuts) needs to be taken into account. The idea is that the invariant mass 
spectra are corrected "into the PHENIX acceptance" as: 

dN 1 1 A^ BB0 (3 23) 



dni ee N ev t Afn ee £p a ir ^bias 

with Am ee the bin width in GeV/c 2 , £ pair = e^ r e^e^ T e pSlir cuts is the pair 
efficiency, £bbc = 0.545 ± 0.06 is the BBC efficiency and £bias = 0.79 ± 0.02 is 
the BBC trigger bias. 

On the other hand, distributions of e + e~ pairs are corrected for the 
limited geometric acceptance (e^£ = £g eo £ geo) azimuth (0 < < 2n) 

and into one unit of rapidity (\y\ < 0.5). The spectra are shown as yield 
invariant under Lorentz transformation: 

1 d 2 N 1 1 N ee e BBC ^ 24 ^ 



2vr p T dp T dy 2np T N cvt Ap T e . e his 



-pair^geo ^mas 

This invariant yield can be converted into an invariant cross section by multi- 

PP 

d 3 a 1 d 2 N incl 
E df ~ 2np T dp T dy app (3 ' 25) 



plying with the inelastic p + p cross section <r™ el = 42.2 nb. 



3.5.1 Tracking and Electron Identification Efficiency 

The tracking and electron identification efficiency is calculated for single elec- 
trons based on a full Monte Carlo simulation of tt° decays. 450M events are 
generated containing 7r°'s flat in rapidity (\y\ < 0.5), azimuth (0 < < 2ir), 
vertex distribution (|z ve rtex| < 30 cm), and in p T (0 < p T < 25 GeV/c). To 
enhance statistics the Dalitz decay branching ratio is set to 100%. 

These events have been processed through a full GEANT simulation pro- 
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Figure 3.24: Single electron acceptance. The diagonal line represent the accep- 
tance limits due to RICH and DC as defined in Eq. (2.5). The dashed lines indicate 
the low pt cut off at 200 MeV/c. Dead areas are removed with fiducial cuts shown 
as red shaded areas. 



gram of the PHENIX detector [108] which includes the details of the detector 
response. The output is analyzed by the PHENIX event reconstruction chain. 
The same single electron cuts are applied in simulation as in real data. Every 
electron is weighted according to a realistic px distribution of the parent tt°. 
The electron identification efficiency is determined as function of the single 
electron px~. 

, _ dN±Jdp± 
£elD ~ dN±/dp± {6 b) 

dN^/dp^ is the px distribution of generated electrons that fall into active area 
within the PHENIX acceptance as parameterized in Eq. (2.5). The active area 
is parameterized from data as function of the single particle charge sign q, px 
and azimuthal angle 0o as shown in Fig. 3.24. Dead removed with 

fiducial cuts which are listed in Tab. 3.4 in the simulations. dN^/dp^ is the 
number of generated electrons in the same fiducial area that pass the electron 
identification cuts. Fig. 3.25 shows the single electron (positron) efficiency as 
function of px in p + p collisions. 

The correction for e + e~ pairs is determined with the hadron decay gener- 
ator Exodus. The single electron (positron) efficiency is parameterized by a 
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Figure 3.25: Track reconstruction and electron identification efficiency in p 
collisions (a)-(b) and in Au + Au (c)-(d) as function of the single electron px- 



Table 3.4: Fiducial cuts in the PHENIX acceptance. 



Dead Area 



A 


-O.88 + 4.560 o <g/pr 


&& q/pr> 


-0.28 


+ 4.560o 


B 


O.28 + 4.560 o <g/pr 


&& q/pr> 


0.58 


+ 4.560o 


C 


O.52 + 4.860 o <g/pr 


&& q/pr> 


0.88 


+ 4.860o 


D 


-2.68 + 4.560 o <g/pT 


&& q/pr> 


-2.38 


+ 4.560Q 


E 


-13.78 + 4.560 o <g/pT 


&& q/pr>~ 


-13.48 


+ 4.560o 
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Table 3.5: Parameterization of the single electron efficiency. 





V 


+ p 


Au - 


- Au 




electrons 


positrons 


electrons 


positrons 


a 


0.991 


-0.662 


0.824 


0.818 


b (GeV/c) 


-0.244 


-0.0400 


-0.168 


-0.182 


c (GeV/c) 2 


0.00963 


-0.00431 


8.12xl0~ 4 


8.88xl0~ 4 


d (GeV/c)" 1 


-0.0260 


-0.00349 


-0.0324 


-0.0301 


e 


-0.478 


0.417 


-0.340 


-0.236 


/ (GeV/c)- 1 


-2.01 


5.60xl0~ 9 


-3.14 


-3.15 



fit of the following functional form: 

f(pr) =a + b/p T + c/p 2 T + dp T + e {e+fPT) (3.27) 

The fit values for electrons and positrons in p + p and Au + Au collisions, 
respectively, are tabulated in Tab. 3.5. This parameterization has been imple- 
mented in Exodus and used as a weight to each electron. In addition e + e~ 
pairs are rejected if at least one of the electrons falls in one of the dead areas 
defined in Tab. 3.4. Then the efficiency to reconstruct and identify an e + e~ 
pair is determined as function of the pair px and the invariant mass. 

This method has been cross checked with a full Monte Carlo simulation of 
e + e~ pairs. About 1M e + e~ pairs have been generated, one half of them flat 
in mass (0-4 GeV/c 2 ) and px (0-4 GeV/c), the other half linearly falling in 
mass (0-1 GeV/c 2 ) and (0-1 GeV/c) to enhance the statistics at low mass 
and px where the efficiency varies the most. The pairs have been simulated 
flat in azimuth (0 < <fi < 2ir) and rapidity (\y\ < 0.5). Only pairs with 
both electron and positron in the ideal PHENIX acceptance were processed 
by GEANT, reconstructed and analysed with the same cuts as the data. As 
with the previous method, the efficiency is calculated as function of the e + e~ 
pair px and mass. 

The efficiency to reconstruct and identify an e + e~ pair is shown as function 
of its invariant mass in Fig. 3.26 averaged over all pt assuming the pt shape 
from the hadron decay cocktail. At high masses the pair efficiency of ~ 45% 
is the square of the plateau value of the single efficiency which consist out 
of the efficiency of the track reconstruction (90%), the electron identification 
(75%) and the dead areas (65%) within the acceptance. Towards smaller 
masses the shape of the pair efficiency is the result of the convolution of the 
single electron efficiency which steeply falls towards low px and the geometrical 
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Figure 3.26: Track reconstruction and electron identification efficiency in p + p 
collisions (a) and in Au + Au (b) as function of e + e~ pair invariant mass. 



acceptance which cuts the single electron spectrum at pr < 200 MeV/c. In 
the mass range 0.4 < m ee < 0.8 GeV/c 2 the pair efficiency drops as a result of 
the single efficiency decrease at low single electron momenta. But for masses 
below 0.4 GeV/c 2 the geometrical acceptance removes all pairs with itlt = 
\J p\ + m 2 e < 0.4 GeV/c 2 , which leads to a higher average px for e + e~ pairs 
in this mass window and therefore a rise of the efficiency. 

The effect of the pair cuts on legitimate signal pairs is shown separately in 
Fig. 3.26. The efficiency is due to legitimate pairs that accidentally fulfill one 
of the pair cut criteria, e. g., share the same ring in the RICH and are therefore 
removed. This efficiency is studied with mixed events in which any detector 
overlap is accidental by definition. The pair cut efficiency is calculated as: 

d 2 A^ pair cuts / d 2 N ■ 

_ " iv mix / a iv mix /„ 9 o\ 

^pair cuts — , , / , , [O.ZO) 

am ee apT / am ee apT 



where d 2 N^" cuts / (dm ee dpT) is the mixed event distribution of e + e pairs with 
pair cuts and d 2 N mix / (dm ee dpr) the one without pair cuts. 

The 0y removes not only photon conversions but also signal pairs which 
happen to be oriented in the magnetic field like conversions. This effect is stud- 
ied with the full GEANT Monte Carlo simulation of 7r°'s which was already 
used for the electron identification efficiency. Above the 7r° mass this is ex- 
tended with the help of hadron decay simulations with Exodus. To simulate 
the detector resolution, an empirical smearing of the magnitude and direction 
of the 3-momentum vector of the electrons was implemented to reproduce the 
<ftv distribution observed in data and the GEANT simulation. The Gaussian 
smearing of the azimuthal and polar angle of the 3-momentum vector p of a 
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Figure 3.27: Comparison of (j)y distributions for various mass regions in data 
(Au + Au), Exodus and the GEANT simulation. 



track with charge sign q — ±1 was parameterized as follows: 



( 0.16\ 

lMf = q- 0.002 ( 0.021 + -pyj (3.29a) 

I 0A6 

= 0.0023 J 5.1 + -pj (3.29b) 

^ = 0.0 (3.29c) 



0.54+^- (3.29d) 

The resulting <pv distribution is compared for different mass ranges to data 
and the full GEANT simulation in Fig. 3.27. The efficiency is calculated from 
the ratio of e + e~ pairs with and without applying the <fiy cut and is shown 
together with the other efficiencies in Fig. 3.26. 



3.5.2 ERT Efficiency 

The dielectron analysis in p + p requires that at least one electron in every 
event is associated with a hit in the ERT trigger. This biases the distribution 
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Figure 3.28: px distribution of single electrons (black) and single electrons which 
are associated with an ERT trigger hit (red). 



of e + e~ pairs and needs to be corrected to allow for a meaningful comparison. 
The trigger efficiency of a single electron is determined with the MB data 
sample. The Level-1 trigger decision is recorded also in the MB data set. It is 
possible that an event is not registered as ERT events, i. e., the Level-1 ERT 
trigger bit is not set, but contains an electron associated to a hit in the ERT 
tiles. Such events would lead to a miscounting of electrons which fired the 
trigger, therefore, only MB events are considered in which the ERT trigger bit 
is set. Within those the pt spectrum of electron candidates <^V"mb&&ert / ^Pt > 
i. e., tracks that fulfill the electron identification cuts, that can be associated 
to a fired ERT trigger tile is shown in Fig. 3.28. The distribution is compared 
to the pt spectrum of all electrons candidates in all MB events dN^ B /dp^ 
(i. e., the ERT trigger bit is not required). 

The trigger efficiency is defined as the ratio of the two as in Eq. (3.30), 
which is shown in Fig. 3.29 for each EMCal sector individually as function of 
the single electron px 



-ERT 



C ^MB&&EF 



/dp^ 



(3.30) 



The trigger efficiency rises steeply above the trigger threshold of 400 MeV 
and reaches half of it's plateau value at pr ~ 600 MeV/c 2 . The plateau value 
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itself strongly depends on the EMCal sector, due to the variation in the number 
of active trigger tiles. The trigger efficiencies are fit to two Fermi functions 4 , 
one below and the other above a single electron p? of 0.5 GeV/c: 

f(p T ) = 6(p T - 0.5) + — %z 6(pt + 0.5) (3.31) 

J ' e -(PT-P0)/k + 1 K ' e -(PT-Po)/k' _|_ I KyI ' I K I 

The fit results for every sector are tabulated in Tab. 3.6 

It is important to note, that to determine the ERT efficiency more stringent 
electron identification cuts need to be applied to reduce the hadron contam- 
ination in the electron sample to a negligible amount. This is necessary as 
any hadron in the electron sample lowers the apparent single electron trig- 
ger efficiency. But the contribution of hadrons to the e + e~ pair spectrum is 
subtracted with the combinatorial background, therefore, the hadron contam- 
ination in the pair spectrum is much reduced to a negligible fraction. It is the 
trigger efficiency of the electrons in this background subtracted sample, that 
are of interest. Therefore, the following elD cuts are applied: 

• 0.20 GeV < p T < 20.0 GeV 

4 for EMCal sector E0 it was sufficient to fit the ERT efficiency with one Fermi function 
over the full pr range 
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Table 3.6: Parameterization of the ERT efficiency for single electrons. 



Sector 




Po (GeV/c) 


k (GeV/c) 


4 


Po (GeV/c) 


k' (GeV/c) 


EO 


0.297 


0.326 


16.7 








El 


0.386 


0.331 


19.5 


0.553 


0.394 


9.61 


E2 


0.415 


0.425 


27.8 


0.790 


0.493 


8.15 


E3 


0.465 


0.418 


24.3 


0.753 


0.470 


9.92 


WO 


0.396 


0.439 


23.1 


0.675 


0.499 


9.08 


Wl 


0.460 


0.413 


26.6 


0.886 


0.471 


6.03 


W2 


0.303 


0.422 


24.3 


0.630 


0.505 


7.85 


W3 


66.7 


1.53 


6.32 


0.281 


0.562 


6.99 



• DC track quality == 63 1 1 51 1 1 31 

• n > 3 

• disp < 5 

• xV n Pe < 10.0 

• 0.8 < E/p < 1.2 

• \J emcsdphLe 2 + emcsdz_e 2 < 2a 

On the other hand, this introduces a bias towards electrons which are more 
likely to fire the ERT trigger, e. g., the harder no cut prefers electrons that emit 
more Cherenkov light in the RICH. This effect has been studied, by varying the 
electron identification cuts. The trigger efficiency is compared for the strong 
cID cuts listed above, the standard elD cuts used for the pair analysis and a 
third set of cuts: 

• 0.20 GeV <p T < 20.0 GeV 

• DC track quality == 63 || 51 || 31 

• n > 2 

• disp < 5 

• xV n Pe < 10-0 

• E/p > 0.8 
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• \J emcsdphLe 2 + emcsdz_e 2 < 3a 

which in its strength lies in between the two other sets. The resulting trigger 
efficiencies are shown in Fig. 3.31. The ratio with respect to the last set of elD 
cuts shows only a small effect for p T > 1.5 GeV/c, but for p T < 1.5 GeV/c 
differences up to 20% appear, which is assigned as systematic uncertainty on 
the trigger efficiency. 

The efficiency of an e + e~ pair is determined analog to the calculation of the 
electron identification efficiency with Exodus. e + e~ pairs from hadron decays 
are generated and the EMCal sector that an electron would hit is determined 
as function of its px and azimuthal emission. Then its probability to fire 
the ERT trigger is given by the trigger efficiency in this EMCal sector which 
has been parameterized from the fit results in Tab. 3.6. For an e + e~ pair to 
survive, at least one of the two electrons must have fired the ERT trigger. 

The ratio of e + e~ pairs surviving the ERT trigger condition and all pairs is 
determined in mass and px of the pair. For the pt inclusive mass spectrum a 
trigger efficiency is applied in mass only as shown in Fig. 3.32. The projection 
onto the mass axis is weighted by the pt distribution of e + e~ pairs from 
Exodus. The plateau at large invariant mass corresponds to the average 
single electron trigger efficiency at high p?- The shape is similar to the elD 
efficiency, but has a more pronounced minimum, as the trigger condition needs 
to be matched by only one of the electrons. 



Table 3.7: Parameterization of trigger efficiency as function of mass. 



a (c/GeV) 


b (c 2 /GeV 2 ) 


c (c 3 /GeV 3 ) 


d (c 4 /GeV 4 ) 




0.349 


0.861 


-2.47 


1.21 






e po (GeV/c) 


k (GeV/c) 




Po (GeV/c) 


Jfe' (GeV/c) 


0.260 1.14 


3.79 


0.585 


0.408 


7.24 



The mass range 1 < m ee < 3 GeV/c 2 has large statistical uncertainties due 
to the absence of hadron decays in this region. Therefore, and to smoothen 
other statistical fluctuations the result is parameterized by a fit. For m ee < 
400 MeV/c 2 a third order polynomial is fit. Above, the sum of two Fermi 
functions is fit to the ERT efficiency. The result is summarized in Tab. 3.7. 

For the differential analysis in mass and pt, the two dimensional trigger 
efficiency distribution is smoothened by fitting the pt distribution of e + e~ 
pairs in narrow mass slices in the LMR, the J/ip and ip'. The IMR is fitted 
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Figure 3.30: ERT trigger efficiency for single electrons with different elD cuts. 
The trigger efficiency for strong elD cuts is shown in (red), for medium elD cuts 
(blue), and for lose elD cuts in (black). 
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Figure 3.31: Shown are the ratios of ERT trigger efficiencies with different elD 
cuts. Strong divided by medium is shown in (red), lose divided by medium in (black). 
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Figure 3.32: ERT trigger efficiency for e + e pairs as function of m ee . The fits are 
summarized in Tab. 3.7. 

in slices of pr analog to the p T inclusive trigger efficiency in Fig. 3.32. The 
results are combined into a smoothed parameterization of the ERT efficiency 
as function of mass and Pt which is shown in Fig. 3.33. 

The uncertainty in the trigger efficiency due the electron identification 
cuts of 20% at low px affects e + e~ pairs with ttit < 1 GeV/c 2 , for which this 
uncertainty is assigned. Elsewhere a 5% systematic uncertainty is assigned. 
This includes uncertainties due in the precise location of dead trigger tiles, 
which has been studied by shuffling EMCal sectors in the simulation of the 
ERT pair efficiency. 

Figure 3.34 compares of the invariant mass spectra for the p + p data ob- 
tained with the MB and ERT data sets. The ERT data set has been corrected 
by ERT trigger efficiency and elD efficiency, while the MB is corrected only 
for the elD. The agreement between the two data sets confirms the assessed 
accuracy of the ERT correction. 

3.5.3 Occupancy Correction in Au + Au 

In Au + Au collisions there is also a finite efficiency loss for particle detec- 
tion due to the presence of other particles nearby. To get a quantitative un- 
derstanding of the multiplicity-dependent efficiency loss, single electrons and 
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Figure 3.33: ERT trigger efficiency for e + e pairs as function of m ee and px- 



positrons are generated with a full GEANT Monte Carlo simulation of the 
PHENIX detector and then embedded into data files containing detector hits 
from real Au + Au events. Next, these new files containing the embedded e^ 1 
are run through the entire reconstruction software to produce track candidates 
containing the variables used for the electron identification cuts. 

Since all the detectors used in the analysis are located after the pair has 
been opened by the magnetic field, the pair embedding efficiency is defined as 
the square of the single electron efficiency 

e embed = ( 5 embcd) (3.32) 

/ # reconstructed e from embedded data 
\# reconstructed e ± from single track data 

where a reconstructed particle from embedded data has most of its detector 
hits associated with hits from the simulated particle. 

Tab. 3.8 lists the embedding efficiencies for the centrality classes used in 
the analysis. For the minimum bias the occupancy correction is weighted by 
the effective pair signal in each centrality class; since most of the yield is 
concentrated in the most central classes, the resulting pair efficiency is 0.81 
instead of 0.92 = 0.96 2 . A systematic uncertainty of ~ 3% has been assigned 
by calculating the occupancy efficiency with a data-driven method [109]. 
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Figure 3.34: Invariant mass spectra for p + p data with the MB (hollow) and ERT 
(solid) data sets. The agreement between the two data sets is excellent. 
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Table 3.8: Embedding efficiency. 



Centrality 


e ± 

embed 


pair 
embed 


0-10% 


0.86 


0.74 


10-20% 


0.91 


0.83 


20-40% 


0.93 


0.87 


40-60% 


0.97 


0.95 


60-92% 


0.99 


0.98 


0-92% 


0.96 


0.81 



3.5.4 Acceptance 

In addition to corrections for the track reconstruction, electron identification 
and ERT trigger efficiency, the pt spectra of e + e~ pairs are also corrected for 
the geometrical acceptance of the two central arm detectors e^J. It accounts 
for the fraction of e + e~ pairs over full azimuth < <p < 2tt and one unit of 
rapidity (\y\ < 0.5) of which at least one electron was not within the detector 
acceptance. 

The acceptance correction has been calculated based on a Monte Carlo sim- 
ulation of Dalitz decays of pseudo-scalar mesons (ir°, r], rf) and direct decays 
of vector mesons (p, uj, <f), J/ip, ip')- The rapidity distribution of all mesons is 
assumed to be flat, which is well justified by measurements of BRAHMS [110]. 
The pt distribution of these mesons is taken from PHENIX measurements (see 
Section 3.7), their polarization from the PDG book [29]. Vector mesons are 
assumed to be unpolarized. 

The acceptance correction is calculated as all other corrections in mass and 
Pt by dividing the number of e + e~ pairs within the detector acceptance, i. e., 
both tracks fulfill Eq. (2.5), by all pairs in < < 2tt and \y\ < 0.5. The 
systematic uncertainty on the polarization and due to the pt parameterization 
is studied with a simulation of unpolarized pairs with a flat mass and pt distri- 
bution, 5 . The resulting acceptance corrections are compared in Fig. 3.35 for 
a number of mass ranges as function of pt- The shapes of the two acceptance 
corrections are very similar and agrees within 5% for the lowest mass bin, and 
better at higher masses. 

In the intermediate mass region (1-3 GeV/c 2 ) charmed meson decays dom- 
inate the e + e~ invariant mass spectrum. Their contribution has been simu- 
lated with Pythia and normalized to the measured cc cross section measured 
in Ref. [16]. These e + e~ pairs do not share the same parent, both are produced 

5 which has already been used as cross check for the cID efficiency 
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Figure 3.35: Acceptance of e + e pairs as function of pt for various mass ranges: 
0-0.1, 0.1-0.2, 0.2-0.3, 0.3-0.4, 0.4-0.5, 0.5-0.6, 0.6-0.7, 0.7-0.8, 0.8-0.9, 0.9-1.0, 
1.0-1.05, and 2.9-3.3 GeV/c 2 (from left to right, top to bottom). The acceptance 
based on the hadronic cocktail is shown in black, the one based on white pairs in 
blue. 
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by semi-leptonic decays of charmed mesons, e. g. D, but their parents are cor- 
related as they contain each one of the c quarks which are always produced in 
pairs, cc — > DD — > e + e~ . For masses m ee < 0.5 GeV/c 2 the charm contribu- 
tion is negligible, above a systematic uncertainty of 5% has been added due 
to the uncertainty in the charm cross section of 30%. 



3.5.5 Bin Shift Correction 

For the px spectra of e + e~ pairs a further correction is applied to take into 
account the difference between the average pt of e + e~ pairs within a bin and 
the bin center due to the finite bin width. 

The bin shift correction of the extracted yield is determined by iteratively 
fitting the dN / dpT spectra with a power law and shifting the yield in every pt 
bin the fits by the relative difference between the yield at the bin center (p^ nt ) 
and at the average pt of the bin (pr)- Three iterations are performed, but 
already after the first iteration the bin shift correction changes insignificantly. 
As example, the fit results are shown in Fig. 3.36 for the u and <fi mesons yields 
(also see Section 4.1.2). Instead of correcting the bin center for the bin shift, 
the yield is adjusted by a factor: 



p«nt+A/2 



cent_ A In 

SWT) = ( 3 - 34 ) 

where A is the bin width. The bin shift correction is applied on the final 
binning and only after this correction is the yield dN/dpx converted into an 
invariant yield or cross section Ed 3 a/dp 3 . 

The correction for the u yield in p? bins of 200 MeV/c is ~ 3.28% below 
3 GeV/c and ~ 1.8% above. The bin shift correction for the (f) cross section is 
~ 2.2% in the full range of 0-5 GeV/c. 
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Figure 3.36: Yield of u (left) and <j) (right) before bin shift correction. The data 
are fit to a power-law to determine the bin sift correction. The first iteration is 
shown in red, the third iteration in green. Open points show the bin shift corrected 
bin centers. 



3.6 Systematic Uncertainties 

The various contribution to the total systematic uncertainty, which are sum- 
marized in Tab. 3.9, are described in the following, (i) The uncertainty on 
the e + e~ pair reconstruction, which is twice the uncertainty on the single elec- 
tron reconstruction efficiency [15, 16], includes uncertainties on the efficiency 
of electron identification cuts, geometrical acceptance, and run-by-run fluc- 
tuations; the uncertainty on (ii) the conversion rejection cut, (iii) the ERT 
and minimum bias trigger efficiencies. These uncertainties are included on the 
yield of e + e~ pairs in the invariant mass spectra. The uncertainties on the 
track reconstruction do not have a strong dependence on the e + e~ pair pj>, 
therefore they are assigned independent of px- Pair and conversion cuts are 
localized in mass (m ee < 600 MeV/c 2 ) and also rather p? independent. 

The largest contribution is the systematic uncertainty on the background 
subtraction. The mixed event normalization enters with 3% x B/S in p + p 
and in Au + Au with 0.25 x B/S. But the background-to-signal ratio, which is 
strongly mass and pt dependent, is two orders of magnitude worse for Au + Au 
collisions than in p+p as shown in Fig. 3.37. As the signal-to-background ratio 
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Au+Au 60-92% 

□ p+p 



(b) 



Figure 3.37: Signal-to-Background ratio for e + e~ pairs in min. bias Au + Au 
collisions with different px ranges (a) and for different centralities (b) which also 
includes e + e~ pairs in p + p collisions. 



shows a px dependence, it is assigned for each p? bin individually. 

For the pt spectra the uncertainty on the acceptance correction of 10%, 
and for m ee > 0.5 GeV/c 2 an additional 5% due to the charm contribution, 
are included in the systematic uncertainty. 
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Table 3.9: Systematic uncertainties of the dilepton yield due to different sources 
and indication of mass range where the error is applied. 



Syst. Uncert. component 


p + p 


Au + Au 


Mass Range (GeV/c 2 ) 


pair reconstruction 


14.4% 


13.4% 


0-4 


conversion rejection 


6.0% 


6.0% 


0-0.6 


pair cuts 


5.0% 


5.0% 


0.4-0.6 


occupancy efficiency 




3.0% 


0-4 


ERT efficiency 


5.0% 




TTlrp > 1 




20.0% 




rax < 1 


combinatorial background 


3.0%-B/S 


0.25%-B/S 


0-4 


correlated background ( "near-side" ) 


3.0% 


10.0% 


0-0.6 


correlated background ( "away-side" ) 


11.0% 




> 1 


centrality 




10.0% 


0-4 


acceptance correction 


10.0% 


10.0% 


0-4 


charm acceptance 


5.0% 


5.0% 


>0.5 



3.7 EXODUS Simulation of Hadron Decays 

Exodus is a fast Monte Carlo simulation of hadron decays. Key input is the 
rapidity density dN/dy of neutral and charged pions which is determined by 
a fit to PHENIX data [87] with a modified Hagedorn function: 

E % = A { e ~ {aPT+bpl) + Pt/Po) ^ (3-35) 

with A = 377 ± 60 mb GeV~ 2 , a = 0.356 ± 0.014 (GeV/c) _1 , b = 0.068 ± 
0.019 (GeV/cT 2 , p = 0.7 ± 0.02 GeV/c, and n = 8.25 ± 0.04. 

PHENIX cross section measurements of 7r°, 7r =1= , t], rf, uj, <p, and J/i/j are 
shown in Fig. 3.38. The fit of the charged and neutral pion data with Eq. (3.35) 
is also shown. 

To parameterize the other hadron spectra, the modified Hagedorn fit of the 
pion data is m T scaled, replacing p T in Eq. (3.35) by J (pr/c) 2 — m 2 + m 2 , 
where rrth is the mass of the hadron and fit with a free normalization factor 
to the hadron spectra. The good description of the data by the fits indicates 
that meson production in p + p collisions follows scaling. Accordingly, this 
observation is extended to other mesons for which no data are available. 

In order to extract the meson yield per inelastic p + p collision the fits are 
integrated over all px- Results, systematic uncertainties, and references to data 
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are given in Tab. 3.10. For the p meson Opja^ = 1.15 ± 0.15 is assumed, con- 
sistent with values found in jet fragmentation [29]. The 7/ yield is scaled to be 
consistent with jet fragmentation o n ijo n = 0.15 ±0.15 [29] which is consistent 
with preliminary measurements by PHENIX [111]. The ip' is adjusted to the 
value of (T^pr/aj/^ = 0.14 ± 0.03 [112] in agreement with preliminary PHENIX 
results [113]. For the r], lo, (ft, and J ftp the quoted uncertainties include those 
on the data as well as those using different shapes of the pt distributions. 
Specifically their spectra have been fitted with the functional form given in 
Eq. (3.35) with all parameters free and also an exponential distribution in rriT- 
For the p, 7/, and ip' the uncertainty is given by the uncertainty assumed for 
the cross section ratios. It is important to note that the dilepton spectra from 
meson decays are rather insensitive to the exact shape of the pr distribution. 



Table 3.10: Hadron rapidity densities in p + p collisions at ^fs = 200 GeV used in 
Exodus. 



dN 
dy 


y=0 


re\. uncert. 


meson/V 


data used 




1.065 ± 0.11 


10% 


1.0 


PHENIX [87, 114] 


V 


0.11 ± 0.03 


30% 


0.1032 


PHENIX [115] 


P 


0.089 ± 0.025 


28% 


0.0834 


jet fragmentation [29] 


LO 


0.078 ± 0.018 


23% 


0.0732 


PHENIX [116-121] 





0.009 ± 0.002 


24% 


0.0084 


PHENIX [116-118, 120, 122] 


if 


0.016 ± 0.005 


40% 


0.0127 


PHENIX [111] and [29] 


J/ij 


(1.77 ±0.27)xl0~ 5 


15% 


0.0000166 


PHENIX [123] 


if)' 


(2.5 ±0.7)xl0" 6 


27% 


0.0000023 


world average [112] 



Once the meson yields and px spectra are known the dilepton spectrum 
is given by decay kinematics and branching ratios, which are implemented in 
our decay generator Exodus following earlier work published in [16]. The 
branching ratios are taken from the compilation of particle properties in [29]. 
For the Dalitz decays 7r°, 77, rj' — > 7e + e~ and the decay uj — > 7r°e + e~ the 
Kroll-Wada expression [50] is used with electromagnetic transition form factors 
measured by the Lepton-G collaboration [51, 126]. For the decays of the vector 
mesons p, to, <p — > e + e~ the expression derived by Gounaris and Sakurai [127] 
is used, extending it to 2 GeV/c 2 , slightly beyond its validity range. For 
the J /if) and ip' — > e + e~ the same expression modified to include radiative 
corrections is used as discussed in [123]. All vector mesons are assumed to be 
unpolarized. For the Dalitz decays of which the third body is a photon, the 
angular distribution is sampled according to 1 + A cos 2 #cs distribution, where 
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Figure 3.38: Compilation of meson production cross sections in p + p collisions at 
y/s = 200 GeV. Shown are data for neutral [87] and charged pions [114], rj [115], 
kaons [114], u [116-121], </> [116-118, 120, 122], rf [111], and J/i/> [123]. The data are 
compared to the parameterization based on mj- scaling used in our hadron decay 
generator. 
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p T (GeV/c) 

Figure 3.39: Compilation of meson production cross sections in Au + Au collisions 
at ^/stvat = 200 GeV. Shown are data for neutral [46] and charged pions [124], 
r] [115], kaons [124], u [120], (j) [120, 125] and J /if} [60]. The data are compared to 
the parameterization based on tut scaling used in our hadron decay generator. 



116 



Table 3.11: Hadron rapidity densities in Au + Au collisions at y / S7v r /v = 200 GeV 
used in Exodus. 



rel. uncert. meson/V data used 



7T° 


97.72 ± 9.5 


10% 


1.0 


PHENIX [46, 124] 


V 


10.77 ± 3.2 


30% 


0.112 


PHENIX [115] 


P 


8.60 ± 2.6 


30% 


0.08981 


jet fragmentation [29] 




9.88 ± 3.0 


30% 


0.1032 


PHENIX [120] 





2.05 ± 0.6 


30% 


0.0214 


PHENIX [120, 125] 


if 


2.05 ± 2.05 


100% 


0.02146 


jet fragmentation [29] 


J/t/j 


(1.79 ± 0.26)xl0~ 5 


15% 


0.0000182 


PHENIX [60] 


t' 


(2.6 ±0.7)xl0" 6 


27% 


0.0000027 


PHENIX [113] and [112] 



#cs is the polar angle of the electrons in the Collins-Soper frame. 

The resulting systematic uncertainties depend on mass and range from 
10 to 25%. They are combined of the uncertainty of the measured 7r° yield 
and the meson-to-pion ratios which is the dominant uncertainty. They also 
include uncertainties on the measured electromagnetic transition form factors, 
in particular for the uj — > 7r°e + e~ decay whose contribution, however, is only 
visible in the mass range around 500 < m ee < 600 MeV/c 2 . 

The e + e~ pairs from hadron decays are filtered into the PHENIX accep- 
tance as defined in Eq. (2.5) and their momenta and angular distributions are 
smeared to take into account the detector resolution, as defined in Eq. (3.4) 
and Eq. (3.29). The resulting invariant mass spectra are shown in Fig. 3.40 for 
p+p collisions and in Fig. 3.41 for Au + Au collisions, respectively. The contri- 
bution of the various sources are shown individually and include a contribution 
from open charm decays which is discussed in Section 4.1.1. 
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Figure 3.40: Expected invariant mass spectrum of e + e pairs from hadron decays 
(cocktail) in Au + Au collisions at y/s = 200 GeV from Exodus. 
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Figure 3.41: Expected invariant mass spectrum of e + e pairs from hadron decays 
(cocktail) in Au + Au collisions at ^/snn = 200 GeV from Exodus. 
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Chapter 4 
Results 



This chapter summarizes the various results of the dielectron continuum anal- 
ysis. At first the invariant mass spectrum of e + e~ pairs in p + p collisions at 
y^s = 200 GeV is shown in Section 4.1 followed by the measurement of the 
charm cross section via the e + e~ pair yield in the intermediate mass region 
in Section 4.1.1. Section 4.1.2 presents the pt spectra of the u and mesons 
measured via their resonance decays into e + e~ pairs. The results of the di- 
electron continuum analysis of p + p collisions are concluded with p? spectra 
of the low mass continuum and the contributions of virtual direct photons 
in Section 4.1.4. While the analysis of p + p collisions is certainly the major 
achievement of this thesis, a detailed comparison to the dielectron continuum 
measured in Au+Au collisions is given. The results are presented in Section 4.2 
and compared to theoretical models in Section 4.3. 



4.1 The Dielectron Continuum in p + p Colli- 
sions 

The invariant mass distribution of e + e~ pairs calculated according to Eq. (3.23) 
is compared to the hadronic cocktail as shown in Fig. 4.1. The contribution 
from semileptonic heavy quark decays is discussed in the Sec. 4.1.1. The 
data are well described by the expected contributions from hadron decays and 
semi-leptonic decays of heavy flavored mesons over the entire mass range. The 
resonance peaks of u, <p, J/ip, and the ip' are well reproduced and their width 
in good agreement with the expected mass resolution. The ratio of data to 
cocktail is unity within the quoted uncertainties as shown in the bottom panel 
of Fig. 4.1. 
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Figure 4.1: Electron-positron pair yield per inelastic p + p collision as function of 
pair mass. Data show statistical (bars) and systematic (shades) errors separately. 
The yield per event can be converted to a cross section by multiplying with the 
inelastic p + p cross section of 42.2 mb. The data are compared to a cocktail of 
known sources. The contribution from hadron decays is independently normalized 
based on meson measurements in PHENIX, the systematic uncertainties are given 
by the error band. The contribution from open charm production is scaled to match 
the data {a cZ = 544 ± 39(stat.) ± 142(syst.) ± 200(model) fib). The inset shows 
the same data but focuses on the low mass region. The bottom panel shows the 
ratio of data to the cocktail of known sources. The systematic uncertainties of the 
data are shown as boxes, while the uncertainty on the cocktail is shown as band 
around 1. 
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4.1.1 Charm Cross Section 



Except for the vector meson peaks, the dilepton yield in the mass range above 
1.1 GeV/c 2 is dominated by semileptonic decays of D and B mesons correlated 
through flavor conservation. 

In order to extract the charm cross section, first the measured e + e~ pair 
yield is integrated in the mass range 1.1 < m ee < 2.5 GeV/c 2 and compared to 
the cocktail yield and the yield from semileptonic decays of open heavy quarks 
(charm and bottom) as well as Drell-Yan calculated with Pythia [102, 104]. 
Pythia 6.205 with cteq5l parton distribution function [106] has been used. 
Following earlier analyses [128, 129] a number of modifications to Pythia 
parameters have been made: 

• parp(91) = 1.5 (fcr), 

• mstp(33) = 1 (use common K factor) 

• PARP(31) = 3.5 (K factor), 

• mstp(32) = 4 (Q 2 scale) 

in addition for charm production: 

• msel = 4 (cc production) 

• pmas(4,1) = 1.25 (m c ), 
for bottom: 

• msel = 5 (bb production) 

• pmas(5,1) = 4.1 (m&), 
and for Drell-Yan: 

• MSEL = 11 (Z or 7* production) 

• PARP(31) = 1.8 (K factor) 

• CKIN(3) = 2.0 (min. parton pr). 

The following cross sections are assumed: 

• a C c = 567 ± 57(stat.) ± 193(syst.) pb [16] 

• a bl = 3.7 ±3.7 pb [130] 
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• (Toy =42 ±42 nb [131] 

The integrated yield of e + e~ pairs in 1.1 < m ee < 2.5 GeV/c 2 in the 
PHENIX acceptance are: 

• measured: y da ta = (4.53 ± 0.28(stat.) ± 0.97(syst.)) x 1(T 8 

• PYTHIA charm: Y C c = 4.08 x 10~ 8 

• PYTHIA bottom: Y bl = 2.62 x 1(T 9 

• Pythia Drell-Yan: Y DY = 2.66 x 1(T 10 

• Exodus cocktail: Y hadr = 3.20 x 10~ 9 



The contribution from the cocktail is essentially only due to the tail of 
the p, which is not known at this high mass. Therefore, an uncertainty of 
100% is assigned. However, since the cocktail contributes only 7% to the data 
in this mass range, the systematic uncertainty on the measured yield due to 
the p contribution is only 7%. Also the bottom cross section from FONLL 
has a 100% uncertainty, which for the same reason translates into only a 6% 
systematic uncertainty on the measured yield. 

The ratio of (Ydata — ^hadr) / + YbB + ^dy) = 0.96 can be used to con- 
vert the charm cross section measured in Ref. [16] which is used as input for 
PYTHIA: 

^ = 123 ± 12(stat.) ± 37(syst.) pb 
dy 

= 567 ± 57(stat.) ± 193(syst.) pb 

into what is measured in this analysis: 

(4.21 ± 0.28(stat.) ± 1.02(syst.)) x 10~ 9 

118.1 ± 8.38(stat.) ± 30.73(syst.) ± 39.5(model) /ib 

544 ± 39(stat.) ± 142(syst.) ± 200(model) pb 

This measurement is an important confirmation of the result reported in 
Ref. [16], in particular, as the STAR Collaboration has reported a charm 
cross section in p + p collisions [132, 133], which is a factor of two larger than 
observed by PHENIX. Currently, the source for this discrepancy is unknown, 
but if PHENIX was to overestimate its single electron identification efficiency 



y y 

1 data 1 hadr 



dy 
a,.-. 
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by a factor of two, this effect would enter quadratically into the e + e~ pair mea- 
surement and therefore lead to a significant difference between the dielectron 
and single electron results. 

The systematic and model uncertainties on this measurement are due to 
a number of sources. As the fraction of e + e~ pairs from correlated heavy 
quark decays at mid-rapidity depends on the dynamical correlation between 
the quarks, additional systematic uncertainties beyond the parameterization 
of the PHENIX acceptance (10%). The momentum of a parton inside the 
colliding hadrons are expected to have a finite "primordial" component kx 
which is perpendicular to the beam direction due to Fermi motion. In Pythia 
the 'primordial k^ parameter is used to define the width of the Gaussian 
distribution of the primordial parton k?- The value of kx affects the azimuthal 
correlation between c and c. While a nominal value of 1.5 GeV/c is used, kr 
has been varied between 1 and 3 GeV/c. The resulting e + e~ pair spectra are 
shown in Fig. 4.2 and the integrated yields are: 
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Due to the variation in the yield an uncertainty of 20% is assigned. 

The uncertainty of the relative abundance of charm quarks and the branch- 
ing ratios of semileptonic decays is 21% [16]. A total branching ratio of c — > e 
of 9.5% ± 1% was calculated from the the particle ratios D + /D° = 0.45 ± 0.1, 
D s /D° = 0.25 ± 0.1, and A c /D° = 0.1 ± 0.05 and the branching ratios ac- 
cording to [29]. The change in the distributions leads to an uncertainty of 
15% [16]. 

Furthermore the systematic uncertainty on the longitudinal correlation of 
the charm quarks (rapidity gap) has been studied by changing the parton 
distribution function (PDF). A PDF ff(x,Q 2 ) describes probability to find a 
parton with flavor i inside a beam particle p participating in a hard-scattering 
process with momentum transfer Q 2 . As the PDF can not be calculated with 
pQCD, PDFs are constructed by various groups by fits to experimental data 
for a fixed momentum transfer and are then extrapolated in the (x, Q 2 ) 
plane. 

The resulting spectra for five different PDFs (cteq5l [106], cteq4l [134], 
GRV94LO [135], GRV98LO [136], and MRST(c-g)[137]) are shown in Fig. 4.3. 
The variation in the e + e~ pair yield in the PHENIX acceptance is 11%. 
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• k T = 1 .0 GeV/c 
1 .5 GeV/c 
2.0 GeV/c 
k T = 2.5 GeV/c 




m ee (GeV/c 2 ) 



Figure 4.2: e + e pairs from semileptonic open charm decays for different values of 
k T : 1.0 GeV/c (Wue), 1.5 GeV/c (Wadfc), 2.0 GeV/c (yellow), 2.5 GeV/c (magenta), 
and 3.0 GeV/c (cyan). 
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Figure 4.3: e + e~ pairs from semileptonic open charm decays for different PDFs: 
CTEQ5L (black), cteq4l (blue), GRV94LO (red), GRV98LO (green), and MRST(c-g) 
(yellow). 
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Figure 4.4: e + e~ pairs from Drell-Yan in \y ee \ < 0.5 from Pythia (black) and a 
NLO pQCD calculation [138] for three scales: \x = q/2 (red), fx = q (green), and 
fi = 2q (blue). 

Drell-Yan 

As estimate of the Drell-Yan contribution to the e + e~ pair spectrum a NLO 
pQCD calculation by Werner Vogelsang has been used [138]. The invariant 
mass spectrum of e + e~ pairs in \y ee \ < 0.5 is shown in Fig. 4.4 using cteq6m 
PDF [139] with the scales ji = q/2, q, and 2q. 

The result is compared to a Pythia calculation which is shown in the 
same figure. While above 4 GeV/c 2 the two calculations disagree in shape, the 
integrated yields agree within 20% in the region of interest. 



Fitting 

In a second approach the charm and bottom cross section have been extracted 
with a simultaneous fit of the mass shapes from Pythia to the cocktail sub- 
tracted data shown in Fig. 4.5. 

The fit is performed in a mass range 1.1 < m ee < 7.0 GeV/c 2 , but excludes 
the mass range 2.5 < m ee < 3.9 GeV/c 2 to avoid influences from imperfections 
in the subtractions of J/ip and ip', with the mass spectra from Pythia for e + e~ 
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pairs from charm, bottom and Drell-Yan: 

dN/dm ee = po ■ dN C c/dm ee + p\ ■ dN bb /dm ee + dN DY / dm ee (4.1) 

where po and p\ denote the two fit parameters and dN C c/dm ee , dN bb /dm ee , 
and dNDY/dm ee are the mass distributions of simulated e + e~ pairs in the 
PHENIX acceptance from charm, bottom and Drell-Yan with the cross sections 
a C c = 567 fib [16], a bb = 3.7 fib [130] and o"dy = 42nb [131], respectively. The 
fit result is p = 0.914 ± 0.082 and p x = 1.064 ± 0.714. When scaling the 
Pythia Drell-Yan spectrum above 4GeV/c 2 to the NLO pQCD calculation 
by W. Vogelsang [138] the results are summarized in Tab. 4.1. 

Table 4.1: Fit results for scaling the Pythia Drell-Yan spectrum to the NLO 
pQCD calculation [138] above 4 GeV/c 2 . 
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The resulting charm and bottom cross sections are therefore: 

a c - c = 518 ± 47(stat.) ± 135(syst.) ± 190(model) fib 
a bl = 3.9 ± 2.4(stat.)jJ(syst.) fib 

The fit result for the charm cross section is in good agreement with yield 
integration method. Although the significance of the bottom cross section is 
limited by large uncertainties, it is the first measurement of bottom production 
at RHIC energies. A preliminary result of the bottom to charm ratio measured 
via electron-hadron correlations [140] leads to a bottom cross section of a bb = 
4.61 ± 1.31(stat.)i2!22( s y s ^-) A*b which is in good agreement with the dielectron 
result. 

Dependence on Collision Energy 

The yfs dependence of charm production in p + p collisions as predicted by a 
NLO pQCD calculations [141] is shown in Fig. 4.6. It is compared to exper- 
imental data from SPS and FNAL experiments at \fs ~ 20 GeV, data from 
PHENIX at -s/s = 200 GeV and results from the UA2 Collaboration at -s/s = 
630 GeV. The two data points from PHENIX are the cross sections measured 
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Figure 4.5: Electron-positron mass distributions from semileptonic decays of heavy 
flavor, obtained by subtracting the contribution from it , rj, u, p, <p, rj , J/tp and ip' 
mesons from the inclusive e + e~ pair yield. The arrows indicate upper limits (95% 
CL) in the mass regions where the charm contribution is smaller or comparable to 
the systematic uncertainties. For all data points statistical error bars and systematic 
uncertainty boxes, including data and model contributions, are shown. Also shown 
are expected contributions from charm, scaled to data, and bottom as well as Drell- 
Yan. 
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Figure 4.6: Shown is a NLO pQCD calculation of the total charm cross section 
as function of the collisions energy yfs. The central value is shown as solid line, 
the upper and lower ends of the uncertainty band are shown as dashed lines. The 
calculation is compared to experimental measurements of the charm cross section 
including the result of this thesis labeled dielectron. 



with single electrons from semi-leptonic heavy flavor decays [16] (labeled spec- 
tra) and the result presented in this thesis (labeled dielectron). Within the 
large theoretical uncertainties, NLO calculations describes the experimental 
data over the full range of y^i. 

For the same NLO calculation, the bottom cross section is compared to 
measurements as function of y/s in Fig. 4.7. The PHENIX data point labeled 
spectra is the result of the electron-hadron analysis [140] and the point labeled 
dielectron is the result of this thesis. NLO calculation. 

Transverse Momentum Dependence 

While it is known that Pythia simulations as used in this analysis fail to 
describe the single electron spectra at high px [16], it is of little concern here, 
as the IMR is dominated by e + e~ pairs from cc decays with small pt but 
large opening angle. Fig. 4.8 shows the measured yield of e + e~ pairs in the 
mass range 1.1 < m ee < 2.5 GeV/c 2 as function of pt- The data are cor- 
rected for electron identification and ERT trigger efficiency, but not for the 
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Figure 4.7: Shown is a NLO pQCD calculation of the total bottom cross section 
as function of the collisions energy y/s. The central value is shown as solid line, 
the upper and lower ends of the uncertainty band are shown as dashed lines. The 
calculation is compared to experimental measurements of the charm cross section 
including the result of this thesis labeled dielectron. 
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Figure 4.8: Measured pt distributions of e + e~ pairs in the mass region 1.1 < 
m ee < 2.8 GeV/c 2 in. The data are corrected for elD and ERT efficiency but not 
for the geometric acceptance. The yield is compared to the expected yield from open 
charm calculated with Pythia filtered into the PHENIX acceptance. No systematic 
uncertainties are shown. 

geometric acceptance. The measured distribution of e + e~ pairs is compared 
to a Pythia calculation filtered into the PHENIX acceptance. While at high 
Pt Pythia fails to describe the pt spectrum of e + e~ pairs, the low pt region 
which contains the dominant fraction of the yield is sufficiently well described. 
Very much analog to the single electron px- It will be interesting to study this 
Pt dependency in more detail and compare to NLO calculations. To leading 
order heavy quark pairs are produced back-to-back, thus their pair pt sums to 
zero. Deviations from this are due to the intrinsic kx carried by the incoming 
partons. In contrast, as shown in Fig. 1.13, next-to-leading order processes 
have three particles in the final state and therefore the distribution of heavy 
quark pairs should be more isotropic and contribute more at high p?. 

Leading Order vs. Next-to-Leading Order Calculations 

In Ref. [142] parton showers in Pythia have been used to approximate some 
of the real corrections to leading order hard scattering processes. For example, 
a final state heavy quark can radiate a gluon or an initial state gluon can split 
into a quark-antiquark pair. According to the number of how many heavy 
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Figure 4.9: Examples of pair creation, flavor excitation, and gluon splitting in 
Pythia. The thick lines correspond to hard scattering processes, the thin ones to 
parton showers [142]. 

quarks are produced in the LO hard scattering process the following processes 
are defined: 

pair creation: Two heavy quarks are created in the final state by quark- 
antiquark or gluon annihilation. These are the LO processes shown in 
Fig. 1.12 

flavor excitation: One heavy quark, which is assumed to come from a gluon 
splitting process, is put on mass shell by hard scattering off a parton 
from the opposite hadron. 

gluon splitting: No heavy quark is involved in the hard scattering, but a 
heavy quark pair is created in initial or final state showers from a gluon. 

The corresponding Feynman graphs are shown in Fig. 4.9 in which thick lines 
correspond to the hard scattering process and thin lines to the parton shower. 

The results of a Pythia calculation including parton showers tuned for 
charm production in Pb + Pb collisions at y/s = 5.5 TeV are shown in Fig. 4.10. 
The total charm cross section has been scaled to a NLO calculation which is 
shown in comparison. Transverse momentum and rapidity distributions of 
single charm quarks, as well as the invariant mass, pt and A<p of cc pair are 
shown. Besides the fundamental differences in the two calculations, the overall 
agreement is rather good. However, especially A<p distributions show significant 
differences. It is very interesting to note that the px distributions of the pair 
creation, which as the LO processes has two hard scattered quarks in the final 
state, shows a significantly softer p? distribution than flavor excitation and 
gluon splitting which have both a more isotropic distribution of the heavy 
quark pair. This underlines the importance of higher order processes at high 
Pt- 
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Figure 4.10: Charm production in Pb + Pb collisions at \fs = 5.5 TeV calculated 
by Pythia including parton showers {solid line). The individual components are 
pair production (dashed), flavor excitation (dotted), and gluon splitting (dot-dashed). 
The Pythia result is compared to a NLO calculation (triangles)[142\. 
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4.1.2 u> and (ft Cross Sections 



As discussed in Section 1.2.4 the in-medium resonance decays of the vector 
mesons uj and provide information about the hadronic phase of the medium 
created in heavy ion collisions at the time of decay. To measure potential 
medium-modifications to the line shape of the resonances, a baseline measure- 
ment of the line shapes in p + p collisions is crucial. Furthermore, the study 
of hadronic versus electromagnetic decays of the are particularly interest- 
ing. As the mass (m^ = 1019 MeV/c 2 ) is very close to twice the kaon mass 
(rriK = 493.7 MeV/c 2 ), the decay — * K + K~ is very sensitive to any changes 
of the mass due to, e. g., the restoration of chiral symmetry. 

To provide this baseline, the cross sections of u and have been mea- 
sured in p + p collisions via their resonance decays into e + e~ pairs. Fig. 4.11 
shows the raw invariant mass spectra of e + e~ pairs for various slices in px 
of 200 MeV/c up to 2 GeV/c and beyond in two bins from 2-3 GeV/c and 
3-5 GeV/c. They are shown together with the the sum of correlated and com- 
binatorial background and the signal after background subtraction. Already 
before background subtraction the uj and peaks are very pronounced. The 
figure also shows a fit to a Gaussian with the mean being fixed to the PDG 
value of the a and mass {m w = 782.65 MeV/c 2 , = 1019.460 MeV/c 2 [29]), 
respectively. 

The ijj and yields (listed in Tab. C.l) are extracted by counting e + e~ pairs 
in a mass region which corresponds to a 3a interval, i. e., 0.740-0.815 GeV/c 2 
and 0.965-1.065 GeV/c 2 , respectively. 

Alternatively, to determine the systematic uncertainty on the peak extrac- 
tion, the peaks are fitted with a Gaussian + n th order polynomial (n = 0, 1, 2) 
to allow for contributions due to the underlying continuum. The fits for the 
uj have been performed in the mass range 0.6-0.9 GeV/c 2 while for the the 
range 0.9-1.2 GeV/c 2 was chosen. The systematic uncertainties on the peak 
extraction have been determined out of the variation of the yield with the 
different fits. For the uj this uncertainty is 14%, while for the it is 20% for 
Pt > 1 GeV/c and 12% below. This uncertainty is added in quadrature to the 
ones listed in Tab. 3.9. 

Analog to Eq. (3.24) the invariant cross section of uj and mesons is 
calculated as 

E d^Vco^ _ _L 1 Nwrf gBBC ^incl / 4 2 ) 

dp* " 2tt p T N evt Ap T B " e pair eg£ e bias pp 1 ' 1 

with the branching ratios for the decays uj — > e + e~: B^ e = 7.18 ± 0.12% and 
-> e+e": Bf e = 2.97 ± 0.04%, respectively [29]. 
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Figure 4.11: Invariant mass spectrum of e + e pairs. Foreground (black), combinatorial background (red) and Signal (blue). 
Also shown are two Gaussian fits (orange) to u and </>, respectively. 



The resulting ui cross sections as function of px are shown in Fig. 4.12. In 
Fig. 4.14 it is compared to existing measurements at higher px via hadronic 
decays, i. e., u — > 7r°7r + 7T~ and ui — > 7r°7 [119-121]. In the px region in 
which both measurements exist, they are in good agreement. Also shown is a 
common fit to a nix scaled modified Hagedorn parameterization of charged and 
neutral pions as defined in Eq. (3.35) with only the normalization factor as free 
parameter, which has been used as parameterization of the u in the Exodus 
cocktail calculation as described in Section 3.7. In addition, alternative fits 
with a modified Hagedorn and all parameters free and a simple exponential in 
m T (£0 = Ae- {mT/T ° s) ) are shown. 

Fig. 4.13 shows the measured 4> cross section, which in Fig. 4.15 is compared 
to the result of the hadronic decay into two kaons [120, 122]. The agreement 
in the region where both measurements are available is not very good and 
currently not understood. A further investigation of the problem is ongoing 
and both results. As for the uj, the data are fit to the ray scaled modified 
Hagedorn parameterization of charged and neutral pions, a modified Hagedorn 
with all parameters left free and an exponential in mx- 

For both, ui and <fi cross sections, the fits with a modified Hagedorn and all 
parameters free results in a shape very comparable to the fit of charged and 
neutral pions. The exponential fit in my offers only a good description at low 
p T but fails for p T > 2 GeV/c. 

4.1.3 Low Mass pt Spectra 

To study the px dependency of the dielectron continuum the invariant mass 
spectra for e + e~ pairs in various px ranges is shown in Fig. 4.16. The results 
are again compared to the cocktail, which has been normalized to the measured 
e + e~ pair yield in m ee < 30 MeV/c 2 . A cocktail with absolute normalization 
agrees with the data in this region within the systematic uncertainty of 20%. 
Again, the p + p data are in good agreement with the cocktail over the full 
mass range and for all px bins, except for some small deviations at high px 
which are discussed in Section 4.1.4. 

In a next step the yield of e + e~ pairs in small mass ranges is analyzed 
as function of px- The yield is now corrected for the geometric acceptance 
of the PHENIX central arm detectors. The invariant yield as function of the 
e + e~ pair p T as defined in Eq. (3.24) is shown in Fig. 4.17 for e + e~ pairs in 
different mass ranges and compared to the cocktail. The px spectrum of e + e~ 
pairs with m ee < 0.1 GeV/c 2 , dominated by 7r° Dalitz decays, is in excellent 
agreement with the expectation of hadron decays. Also the px spectra of the 
other mass windows agree within the systematic uncertainties well with the 
cocktail for px < 2 GeV/c. Above a small excess, which was also visible in the 
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Figure 4.12: Invariant cross section of u vnp + p collisions as function of px- The 
stat. errors are drawn as bars, the syst. uncertainties are represented with the gray 
boxes. 
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Figure 4.13: Invariant cross section of <j) in p + p collisions as function of pr- The 
stat. errors are drawn as bars, the syst. uncertainties are represented with the gray 
boxes. 
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p + p -> co + X 




Figure 4.14: Invariant cross section of w in p + p collisions as function of pt 
measured in electromagnetic and hadronic decay channels. The results are fitted to 
the pion parameterization (blue), a modified Hagedorn (black) and an exponential 
in rriT (red). 




Figure 4.15: Invariant cross section of <f> in p + p collisions as function of pt 
measured in electromagnetic and hadronic decay channels. The results are fitted to 
the pion parametrization (blue), a modified Hagedorn (black) and an exponential in 
nix (red). 
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Figure 4.16: The e + e _ pair invariant mass distributions in p+p collisions. Thepx 
ranges are shown in the legend. The solid curves represent an estimate of hadronic 
sources. 
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high px mass spectra in Fig. 4.16, is observed. This excess is analyzed in the 
next Chapter. 

4.1.4 High p T Direct Photons 

As discussed in Section 1.2.3 any source of real direct photons can also produce 
a virtual photon which subsequently converts internally into an e + e~ pair. 
The mass distribution is given by Eq. (1.23) and is proportional to 1/m for 

p T > m ee . 

In this Chapter an attempt is made to analyze the small excess observed 
at high px in the low mass region of the dielectron continuum, which is shown 
more detailed again in Fig. 4.18, under the assumption that it is solely due 
to internal conversions of direct virtual photons. While at zero mass e + e~ 
pairs from hadron decays have the same shape in mass as internal conversions 
of direct photons, the suppression due to S when approaching the hadron 
mass changes the shape of e + e~ pairs from hadron decays. Therefore, one 
can fit, after an initial normalization to the mass range 0-30 MeV/c 2 , the 
two expected shapes for e + e~ pairs from hadronic decays / coc ktaii (as shown 
in Fig. 4.17) and from direct photons /direct m the mass range 80-300 MeV/c 2 
in which the 7r° contribution is severely suppressed. The only fit parameter is 
the relative fraction of direct photons r: 

f(m ee ) = (1 - r)/ cocktail (m ee ) + r/ dircct (m ee ). (4.3) 

As an example, the fit result is shown in Fig. 4.19 for e + e~ pairs with 
2.0 < p T < 2.5 GeV/c in p + p collisions. The quality of the fit result (e. g., 
X 2 /NDF = 11.8/10 for the lowest pt bin) indeed justifies the assumption that 
the observed excess is due to internal conversion of virtual photons. For each 
Pt bin, the fit is performed for several mass ranges m low < m ee < 300 MeV/c 2 . 
It is assumed that the form factor F{m 2 ee ) for direct photons is 1 independent of 
mass. This may not be the case for direct photons from parton fragmentation 
or, in case of heavy ion collisions, direct photon radiation from a hadron gas. 
If the form factor of direct photons is arbitrarily set to the one of the r] meson, 
the fraction of direct photons r would decrease by ~ 10%. The mass region 
above 300 MeV/c 2 is not used in the fits to avoid possible effects due to the 
form factor and to stay well within the region of m ee <C Pt, m which the 
approximation S = 1 is valid. 

The largest uncertainty in the fit is the particle composition of the cocktail, 
in particular the rj/n ratio, which is 0.48 ±0.03 at high p?, based on PHENIX 
measurements [13, 115]. This corresponds to a 7% uncertainty in the p + p 
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Figure 4.17: Shown is the invariant yield of e + e pairs in different mass ranges in 
p + p collisions. The mass ranges are defined in the legend. 
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Figure 4.18: The e + e~ pair invariant mass distributions in p+p collisions. Thepx 
ranges are shown in the legend. The solid curves represent an estimate of hadronic 
sources. This is a zoomed version of Fig. 4.16. 
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Figure 4.19: Invariant mass distribution of e + e~ pairs in p + p collisions with 
2.0 < p T < 2.5 GeV/c. The fit with Eq. (4.3) in the range 80 < m ee < 300 MeV/c 2 
is explained in the text. The black dashed curve shows f(m ee ) outside of the fit 
range. 
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cocktail for the mass range 100 < m ee < 300 MeV/c 2 . Other sources cause 
only a few percent uncertainty in this mass range, once the normalization of 
the cocktail is fixed to the mass range m ee < 30 MeV/ c 2 . As the r\ meson is the 
dominant hadronic source in this mass range the data has been fit with three 
components f(m ee ) = (1 -r -r r? )/ cocktai i(m ee ) + r/ direct (m ee ) + r r ,f r) (m ee ), with 
a constrain on r v such that the 77/71"° range is not changed beyond 0.48 ± 0.03. 
This alternative method results in values for r which are consistent within the 
statistical uncertainties. 

The fraction of direct photons r is shown in the left panel of Fig. 4.28 
for and compared to a NLO pQCD prediction [41, 138]. The curve shows 
the predicted direct photon cross section da^ LO (pT) divided by the inclusive 
photon cross section da^ cl . The three curves which are shown correspond to 
three different choices of the momentum scale fi = 0.5/>r, l.Opx, 2.0pp. For a 
perturbative treatment, the covered pr range of 1 < pr < 5 GeV/c is very 
low, and the theoretical uncertainties are quite large. Nevertheless, the NLO 
calculation seems in quite good agreement with the data. 

The fraction of direct photons r can be converted into a cross section of di- 
rect photons by multiplying with the inclusive photon spectrum. The inclusive 
photon cross section is determined for each p T bin by da™ cX = a 1 ™ 1 dN^ t& x 
dN° ocktail /dN™ cktail , where diV e data and diV e c e ocktail are the yields of e+e~ pairs in 
m < 30 MeV/c 2 for data and cocktail, respectively, and rfjv™ cktai1 is the yield 
of photons from the cocktail. The resulting direct photon spectrum for p + p is 
shown in Fig. 4.29. The pQCD calculation is consistent with the p + p data for 
Pt > 2 GeV/c within theoretical uncertainties. A similar good agreement is 
observed for n° production in p + p [87]. The p + p data can be well described 
by a modified power-law function: 



as shown by the dashed line in Fig. 4.29. 

In summary the dielectron continuum in p + p collisions at y/s = 200 GeV 
can be fully understood by known hadronic sources, correlated semi-leptonic 
decays of D and B mesons and direct photons produced in hard-scattering 
processes of the incoming partons. The measurement of "quasi-real" direct 
photons in a mass region in which 7r° Dalitz decays are suppressed has allowed 
the first measurement of direct photons at low p?. These results will be crucial 
for the understanding and interpretation of the observed dielectron yields in 
Au + Au collisions. 




(4.4) 
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4.2 Comparison to Au + Au Collisions 



The yield e + e~ pairs in the PHENIX acceptance measured in min. bias 
Au + Au collisions at ^snn = 200 GeV is shown in Fig. 4.20 as function 
of invariant mass. The yield is compared to the expected sources from hadron 
decays from EXODUS and open charm from PYTHIA. The open charm con- 
tribution is scaled to the measured charm cross section in p + p (a C c = 567 ± 
57(stat.) ± 193(syst.) fib) [16] multiplied by the average number of binary 
collisions: (A" coU ) = 257.8 ± 25.4 (see Tab. 3.1). 

The data are well described by the cocktail of hadron decays in the mass 
region below 150 MeV/c 2 . Due to the larger background to signal ratio the 
resonance peaks of uj and <fi are not as well measured as in p + p but within the 
systematic uncertainties in agreement with the expected yield. However, the 
yield of e + e~ pairs in the mass range 150 < m ee < 750 MeV/c 2 is significantly 
enhanced above the expected cocktail yield. The enhancement factor in this 
mass range, defined by the ratio of the measured yield to the expected yield 
from hadronic sources, is 4.0 ± 0.3(stat.) ± 1.5(syst.) ± 0.8(model). Here the 
first uncertainty is the statistical error, the second the systematic uncertainty 
on the data, and the last error the uncertainty on the cocktail. The last 
two uncertainties are very conservative estimates and not equivalent to one 
standard deviation (la). 

In a direct comparison to the measured dielectron continuum in p + p the 
result of min. bias Au + Au is shown in Fig. 4.21a. The invariant mass spectra 
of Au + Au collisions in five centrality bins (0-10%, 10-20%, 20-40%, 40- 
60% and 60-92%, see Tab. 3.1) are shown in Fig. 4.21b and compared to the 
individual cocktail expectations. The low mass enhancement is concentrated in 
the 20% most central collisions and becomes less significant in more peripheral 
ones. The yield in the IMR is consistent with the Pythia calculation for all 
centrality bins. In the following the centrality dependence of the IMR and 
LMR are further investigated. 

4.2.1 The Intermediate Mass Region 

The intermediate mass region which served in the p+p analysis as a window to 
extract the charm cross section is well described by the Pythia calculation of 
open charm pairs. This result is somewhat surprising as the measurement of 
single electron from semi-leptonic charm decays shows a strong modification 
of charm, i. e., a suppression of high px electrons and a significant elliptic 
flow [16]. Therefore, one would expect this modification to also influence the 
pair properties. An extreme scenario is shown in addition to the Pythia 
curve as dotted line in Fig. 4.20 assuming the initial correlation of the cc pair 
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Figure 4.20: Electron-positron pair yield as function of pair mass. Data show 
statistical (bars) and systematic (shades) errors separately. The data are compared 
to a expected sources from the decays of light hadron and correlated charm de- 
cays [23]. The charm contribution expected if the dynamical correlation of the cc 
pair is broken, is shown separately. The systematic uncertainty on the cocktail is 
not shown. 
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Figure 4.21: (a) shows the invariant mass spectrum in p+p and min. bias Au+ Au 
collisions, (b) shows the invariant mass spectra of five centrality classes in Au + Au 
collisions. The data are shown with statistical (bars) and systematic (shades) errors 
separately. The data are compared to the expected sources from the decays of light 
hadron calculated with Exodus and correlated charm decays based on Pythia. 



is completely broken, e + e~ pairs are generated with a single pr distribution 
following the measured spectrum of single electrons from semi-leptonic heavy- 
flavor decays [16], but a random azimuthal opening angle. This distribution 
is much softer, than the Pythia curve, which would leave room to other 
contributions such as thermal radiation via qq annihilation. 

If the yield in the IMR is dominated by open charm, it is expected to in- 
crease proportional to the number of binary collisions. The yield per number 
of binary collision N co n in the mass range 1.2 < m ee < 2.8 GeV/c 2 is shown 
as function of N paTt in Fig. 4.22. The data show no significant centrality de- 
pendence and are consistent with the expected yield calculated with Pythia. 
But the apparent scaling with the number of binary collisions may be a co- 
incidence of two counteracting effects: (i) the suppression of e + e~ pairs from 
open charm in the IMR due to modifications of charm which increases with 
N part and (ii) an additional contribution due to thermal radiation from qq an- 
nihilation which is expected to increase faster than proportional to N part . As 
discussed in Section 1.3.3, such a coincidence may have been observed at the 
SPS [74], where a prompt component has been suggested by NA60 [75]. 
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Figure 4.22: Dielectron yield per binary collision as function of iVp ar t for the mass 
range 1.2 < m ee < 2.8 GeV/c 2 . Statistical and systematic uncertainties are shown 
separately. Also shown are two bands corresponding to two different estimates of 
the contribution from charmed meson decays. The width of the bands reflect the 
uncertainties on the charm cross section only. 
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Figure 4.23: Dielectron yield per participating nucleoli pair as function of N pax t for 
the mass range 150 < m ee < 750 MeV/c 2 (a) and m ee < 100 MeV/c 2 (b). Statistical 
and systematic uncertainties are shown separately. Also shown are the expected 
cocktail yields. The band reflects the systematic uncertainty on the cocktail 



4.2.2 The Low Mass Enhancement 

The centrality dependence of the yield in the low mass region has been studied 
more quantitatively. The first range, m ee < 100 MeV/c 2 , is dominated by 
7T° Dalitz decays. As the ir° multiplicity roughly scales with the number of 
participating nucleons, the yield in this range should scale linearly with N part . 
Fig. 4.23 shows the yield in m ee < 100 MeV/c 2 divided by the number of 
participating nucleon pairs (A^ part /2) as function of A^ part . The yields in all 
centrality bins as well as p + p are in good agreement with the cocktail, and 
show the expected scaling. 

The top panel of Fig. 4.23 shows the yield per participating nucleon pair 
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for the mass region 150 < m ee < 750 MeV/c 2 , in which the enhancement is 
observed. While the agreement with the cocktail is good for p + p and the 
peripheral collisions, an increase of the yield is observed in the data above 
-Wpart ~ 200 that is stronger than linear. For the most central bin, the yield 
is a factor 7.6 ± 0.5(stat.) ± 1.5(syst.) ± 1.5 (model) above the cocktail. This 
stronger than linear rise of the yield is qualitative understood with the notion 
of qq or tttt annihilation as source of the enhancement. 

4.2.3 Direct Photons 

In this chapter it is tested whether the excess observed in the low mass region 
is compatible with a source of virtual direct photons. The analysis has been 
performed with the same method as for p + p which has been discussed in 
Section 4.1.4. Fig. 4.24 shows the yield of e + e~ pairs as function mass below 
1.2 GeV/c 2 for various p? ranges in min. bias Au + Au collisions. The low 
mass enhancement is concentrated at low transverse momenta: pt < 1 GeV/c. 
For pt > 1 GeV/c the enhancement, although reduced, is still quite significant 
(a factor rs 2 above the expected yield), but is quite different in shape. The 
low px mass spectrum divided in finer px bins as shown in Fig. 4.26. While 
the low pt enhancement has a "bump" like structure, the shape at higher pt 
is consistent with the 1 / m ee shape expected from internal conversions of direct 
virtual photons. 

Thus the excess of e + e~ pairs for p? > 1 GeV/c is treated as entirely due 
to internal conversions and as in p + p the mass spectra are fitted with the 
sum of cocktail and a direct photon component with a l/m ee dependence. An 
example for 1.0 < pr < 1.5 GeV/c is shown in Fig. 4.27. The uncertainty 
on the ?7 / 7r° ratio is a bit larger than in p + p, but the same central value is 
assumed: 0.48±0.08, which is the more conservative approach, since the actual 
PHENIX measurement of the ratio Au + Au is ~ 0.4 with an uncertainty of 
11% [13, 115]. This leads to a 17% uncertainty due to the t] contribution in 
the mass range 80 < m ee < 300 MeV/c 2 . As in p + p a fit with a variable rj 
contribution within a rj/ir ratio of 0.48 ± 0.08 gives consistent results within 
the statistical uncertainties. 

The fraction of direct photons r in min. bias Au + Au is shown as function 
of pt in the right panel of Fig. 4.28. It is compared to the NLO pQCD 
prediction of direct photons for p+p divided by the inclusive photon yield and 
scaled by the nuclear overlap factor Taa (Taa ^^ LO (pt)/^V™ c1 (pt)), which 
is also shown for the three theory scales fi = 0.5pT, l.Opr, 2.0p T . While the 
fraction of direct photons measured in p + p was consistent with the NLO 
prediction, the fraction measured in Au+Au is larger than than the calculation 
for p T < 4 GeV/c. 
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Figure 4.24: The e + e~ pair invariant mass distributions in Au + Au collisions. 
The pt ranges are shown in the legend. The solid curves represent an estimate of 
hadronic sources. 
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Figure 4.25: The e + e~ pair invariant mass distributions in Au + Au collisions. 
The pt ranges are shown in the legend. The solid curves represent an estimate of 
hadronic sources. This is a zoomed version of Fig. 4.24. 
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Figure 4.26: The e + e~ pair invariant mass distributions in Au + Au collisions for 
the low pt range. The solid curves represent an estimate of hadronic sources. 
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Figure 4.27: Invariant mass distribution of e + e pairs in Au + Au collisions with 
1.0 < p T < 1.5 GeV/c. The fit with Eq. (4.3) in the range 80 < m ee < 300 MeV/c 2 
is explained in the text. The black dashed curve shows f(m ee ) outside of the fit 
range. 
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Figure 4.28: The fraction of direct photons as function of pt in (a) p + p and 
(b) Au + Au collisions. The error bars and the error band represent statistical and 
systematic uncertainties, respectively. The curves show the direct photon prediction 
of a NLO pQCD calculation. 
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As for p+p, the fraction of direct photons r is converted into a direct photon 
yield by multiplying with the inclusive photon yield. The inclusive photon 
yield is determined for each p T bin by dN'™ 1 = dN^ ta x dN™ cHail /dN™ cHai \ 
where dN^ tS/ and dN™ cktai1 are the yields of e + e~ pairs in m < 30 MeV/c 2 
for data and cocktail, respectively, and ^jV^ ocktai1 is the yield of photons from 
the cocktail. The resulting direct photon spectrum for 0-20%, 20-40% and 
minimum bias Au + Au collisions is shown in Fig. 4.29 together with the direct 
photon cross section measured in p+p collisions. At high px the direct photon 
spectrum is extended by the EMCal measurement of direct photons [40, 45], 
which show good agreement with the NLO calculation. It is worth to point 
out the significant improvement of the systematic uncertainties by the internal 
conversion method with respect to the EMCal measurement 

The measured direct photon yield in Au + Au is above the p + p fit to 
Eq. (4.4) scaled by Taa for px < 2.5 GeV/c, which indicates the the direct 
photon yield at low px increases faster than the binary scaled p + p cross 
section. To quantify this excess, an exponential plus a modified power-law are 
fit to the Au + Au data. The power-law is fixed to the Taa scaled p + p fit: 



1wp T dpx dy 



(2 \ — n 
1 + &j (4.5) 



The inverse slope of the exponential T and the A are the only free parame- 
ters. The systematic uncertainties of this fit are estimated by changing the 
p + p fit parameters within the uncertainties. The results are summarized in 
Tab. 4.2, for which the fit result of A has been converted into a rapidity den- 
sity dN/dy(px > 1 GeV/c). For the most central collision an inverse slope of 
T = 221 ± 23 ± 18 MeV is found. If the direct photons in Au + Au collisions 
are of thermal origin, the inverse slop T is related to the initial temperature 
Ti n i t of the dense matter. In thermal models, T^t is 1.5-3 times T due to 
the space-time evolution [143]. Comparisons to a few models are shown in 
Section 4.3.1. 



4.2.4 pt of the Low Mass Enhancement 

As shown in Fig. 4.24, the dominant fraction of the low mass enhancement 
is localized at low pr- Clearly, in this region the approximation m ee px 
does not hold any longer, and consequently S ^ 1. Therefore, the excess in 
this mass and pr range cannot be analyzed rigorously under the assumption 
of internal conversions of direct photons as done at high i n the previous 
chapter. 

The invariant yield of e + e~ pairs as function of pr is calculated as defined 



156 




Figure 4.29: Invariant cross section (p + p) and invariant yield (Au + Au) of direct 
photons as function of px- The filled point are from this analysis and the open points 
are from [40, 45]. The three curves on the p + p data represent the NLO pQCD 
calculation of direct photons, the dashed curves show a power-law fit to the p + p 
data, scaled by Taa- The solid black curves are exponential plus Taa scaled p + p 
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Table 4.2: Summary of the fits with Eq. (4.5) to the direct photon spectra measured 
in Au + Au at t/iTnn = 200GeV. The first and second errors are statistical and 
systematical, respectively. 



Centrality 


dN/dy(p T > 1 GeV/c) 


T (MeV) 


XVDOF 


0-20% 


1.10 ± 0.20 ± 0.30 


221 ± 23 ± 18 


3.6/4 


20-40% 


0.52 ± 0.08 ± 0.14 


214 ± 20 ± 15 


5.2/3 


MB 


0.33 ± 0.04 ± 0.09 


224 ± 16 ± 19 


0.9/4 



in Eq. (3.24) and compared to the cocktail of hadronic and charmed meson 
decays for various mass ranges in Fig. 4.30. The systematic uncertainties 
are dominated by the background subtraction which contributes as (Ts/S = 
0.25%-B/S*. As the signal to background ratio improves towards high p^, the 
systematic uncertainties decrease. For the mass ranges below 400 MeV/c 2 the 
data are cut at low pt due to the single electron px cut at 200 MeV/c 2 . 

In the mass region m ee < 100 MeV/c 2 the Au + Au data are in agreement 
with the cocktail, as expected. In the higher mass bins, the data show a large 
excess both at low and high pj>. At high px this excess has been interpreted 
as direct photons in the previous chapter. The pt spectra of direct photons 
have been successfully described with the functional form given in Eq. (4.5) 
an the fit values in Tab. 4.2. The yield of direct photons can be converted into 
a yield of e + e~ pairs in a mass range mi < m ee < m 2 according to Eq. (1.23): 



ni2 



1 d 2 N PP 1 d 2 N^ [2a 1 4m 2 / 2m 2 



' ee x 'y 



1 + — ^ Sdm, 



2npT dpx dy 2np T dpx dy J 3ir m ee \j m 2 e \ m 2 e 



in i 



1 + 

Po 



■AA A s±pp 

in-2 



37r m ee \l m 2 ee \ m 2 ee 



Till 



As previously, the distributions of direct photons and cocktail are normal- 
ized to the data in the mass range m ee < 30 MeV/c 2 , where both have nearly 
identical shapes. The extension of the direct photon contribution to lower p? 
is done using the suppression factor S for quark-gluon Compton scattering 
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Figure 4.30: Shown is the invariant yield of e + e~ pairs in different mass ranges in 
min. bias Au + Au collisions. The mass ranges are defined in the legend. 
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Figure 4.31: Shown is the invariant yield of e + e~ pairs in different mass ranges 
in p + p (left) min. bias Au + Au collisions (right). The mass ranges are defined in 
the legend. The data are compared to the cocktail from hadronic decays (solid line) 
and the sum of hadronic cocktail and direct photons (dashed line). 



defined in Eq. (1.26): 

= 1 (1.26) 

(x + VI + x 2 )(3x 2 + 1 + 2x VI + x 2 ) 

with x = pT/m ee . This contribution is added to the cocktail and the sum is 
shown as dashed line in Fig. 4.31, which shows the pr spectra of e + e~ pairs in 
both p+p and Au + Au collisions. The excess at at pr > 1 GeV/c with respect 
to the hadronic cocktail agrees with the sum of cocktail and direct photons. 
For the mass range 100-300 MeV/c 2 this is of course redundant, as the excess 
in exactly this region of phase space had been used to measure this direct 
photon component, but also above 300 MeV/c 2 the excess at px > 1 GeV/c is 
in good agreement with a direct photon contribution. 

At low pt, however, a significant excess beyond the extrapolated direct 
photon contribution remains in the Au + Au data. This excess increases in 
strength towards low p?, despite the restriction on the single electrons to 
have at least a transverse momentum of px > 0.2 GeV/c to be fall within 
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the PHENIX acceptance. This trend is quite remarkable and opposite to the 
behavior of the cocktail component. It is also reminiscent of the steepening of 
the tut spectra of the dimuon excess yield measured by NA60 [78] (also see 
Section 1.3.3). To further study the low excess, the spectra are combined 
to a one mass range of 300 < m ee < 750 MeV/c 2 . 

In Ref. [124] PHENIX has measured that the differential cross section of 
charged hadrons (n ± , K ± , p, and p) in Au + Au collisions at ^/snn = 200 
GeV can be described by an exponential in my — m , where m is the particle 
mass, and = \/p T + m^. This is also known from lower beam energies for 
p + p, p + A, and A + A collisions. Such an exponential shape of the invariant 
yield can be parameterized as: 

1 d 2 N _ 1 ^expf m T-™o \ ^ ^ 

2nmT dm,T dy 2%T(T e ff + m ) \ T c g J 

where T e g denotes the inverse slope or effective temperature, and A is a nor- 
malization parameter which is related to the rapidity density dN/dy. It has 
been observed that the inverse slope increases with the particle mass, which 
is consistent with a hydrodynamical description of particle emission from a 
source with transverse flow [144-148]. Under such conditions the measured 
effective temperature T ef j is blue shifted with respected to the freeze out tem- 
perature T fo : 

T eS = T lo + m(p T ) 2 . (4.8) 

where (/3r) is a measure of the strength of the transverse flow. The inverse 
slopes measured for kaons are larger then 200 MeV for all centralities. Together 
with the effective temperatures measured for charged pions, protons, and anti- 
protons a centrality independent freeze out temperature of ~ 180 MeV has 
been extracted [124]. 

Following this approach to characterize the low pt excess, the px spectrum 
of e + e~ pairs with 300 < m ee < 750 MeV/c 2 is fitted with the sum of hadronic 
cocktail, direct photons and an my exponential as defined in Eq. (4.7). The 
direct photon contribution to the cocktail at low px is 10% at m ee = 0, which 
leads to equal contribution of cocktail and direct photon in the fitted mass 
range. The total contribution of cocktail and direct photons is fixed by the 
normalization to the data in m ee < 30 MeV/c 2 . 

The pt spectrum and the fit result are shown in Fig. 4.32, together with the 
individual components. The data are well reproduced by the fit which gives 
an inverse slope of T eff = 88 ± ll(stat.)t? 3 (syst.) MeV with a x 2 /DOF w 1. 
The systematic uncertainty accounts for the uncertainty on the data and the 
uncertainty on the cocktail normalization (20%). The yield attributed to the 
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exponential accounts for more than 50% of the total e + e~ yield. 

If the source of this excess is due to thermal radiation from the fireball, 
which is dominated by tttt annihilation, the excess yield in the LMR should 
exhibit a similar mass dependent effective temperature due to transverse flow 
as the hadrons. The effective temperature is significantly lower than for any 
of the identified hadron spectra [124]. In particular the kaon, with a mass 
comparable to the average invariant mass of the e + e~ pairs in the fitted mass 
region, has an effective temperature more than twice as large. This seems 
to advocate a source of dielectrons that does not exhibit the rise with mass 
typical for a radially expanding source. Also the temperature extracted from 
the direct photons at high p T (T cff = 221 ± 23(stat.) ± 18(syst.) MeV) is 
significantly larger than the one of the excess yield at low p?. 
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Figure 4.32: The figure shows pt spectra for the mass range 300 < m ee < 
750 MeV/c 2 . The spectrum is fitted to the sum of the cocktail expectation (thin 
line), direct photons and an exponential function in wit (thick line). The cocktail 
include contributions from charmed mesons. The dashed line shows the expectation 
from the cocktail plus the direct photon contribution. The exponential function is 
also shown. 
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Table 4.3: The initial Temperature Ti n ; t and formation time to for various models 
predictions of thermal photons in Au + Au collisions at y^sjvTV = 200 GeV. 





T init (MeV) 


To fm/c 


notes 


David d'Enterria et al. 


590 


0.15 


[143] 


Rasanen et al. 


580 


0.17 


[149, 153] 


Srivastava 


450-600 


0.2 


[150, 154] 


Turbideet al. 


370 


0.33 


[48] 


L. Liu et al. 


370 


0.6 


[152] 


J. Alam et al. 


300 


0.5 


[151] 



4.3 Model Comparisons 

In this section the enhanced production of direct photons and the low mass 
enhancement of the dielectron continuum observed in Au + Au collisions are 
compared to theoretical models 

4.3.1 Thermal Photons 

In Fig. 4.33 the direct photon spectrum measured in 0-20% Au + Au collisions 
is compared to a number of hydrodynamical models [48, 143, 149-152]. All 
of these models are in good agreement with the data, and with each other 
within a factor of two. This fact is quiet remarkable as the initial temperature 
and formation times of these models vary in the range of T m i t = 300 MeV, 
To = 0.6 fm/c to T init = 600, To = 0.6 fm/c. Tab. 4.3 lists the initial conditions 
of all the models mentioned here. They are also shown as pairs of (T init , To) in 
Fig. 4.34, which illustrates the apparent anti-correlation between the formation 
time To and the initial temperature T init which produces time averaged direct 
photon spectra all in agreement with the data. 

The general good agreement of the hydrodynamical models with the data 
confirms that the px range of 1-3 GeV/ c is dominated by thermal radiation. To 
disentangle the details of the initial condition and the hydrodynamic evolution 
of the system more detailed studies will be necessary. 

4.3.2 Low mass dileptons 

While the excess of direct photons produced in Au + Au collision in compar- 
ison to p + p seems to be qualitatively explained by thermal radiation from 
a medium with temperatures of 300-600 MeV, the situation for the low mass 
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Figure 4.33: Comparison of direct photon yield Au + Au to theoretical predic- 
tions [48, 143, 149-152]. Note that the cyan curve by Liu et al. has a pQCD 
contribution included, while all other model curves only show the thermal contribu- 
tion. For a comparison to the data above 3 GeV/c a pQCD contribution, e. g. the 
one shown by W. Vogelsang [41, 138] needs to be added. 
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Figure 4.34: Comparison of formation time To and initial temperature Tj n ;t of the 
thermal photon predictions. 



enhancement of dielectrons of the expected cocktail of hadron decays is quiet 
different. 

In Fig. 4.35 the LMR of the dielectron continuum in min. bias Au + Au 
collision at ^/s^n = 200 GeV is compared to a number of predictions in- 
volving different scenarios for in-medium modifications of the vector mesons 
by R. Rapp and H. van Hees [54], K. Dusling and I. Zahed [155], and E. L. 
Bratkovskaya and W. Cassing [156]. As discussed in Section 1.3, an expanding 
fireball model with a modified in-medium spectral function by R. Rapp suc- 
cessfully describes the dilepton spectra measured by CERES (see Fig. 1.24) 
and NA60 (see Fig. 1.26b) [80]. Also the other authors find good agreement of 
their models with NA60[156, 157]. Both, Rapp and Dusling include a contri- 
bution from qq annihilation during the QGP phase in their predictions which, 
however, is important only in the intermediate mass range. 

Predictions for different p spectral functions (vacuum, dropping mass, and 
collisional broadening) are filtered into the PHENIX acceptance and added 
to the hadronic cocktail which beforehand has its vacuum p contribution sub- 
tracted as this is yield is part of the prediction. The results by R. Rapp [54] are 
shown as blue lines in Fig. 4.35, the vacuum rho is shown as dashed, the drop- 
ping p as dotted and the broadening scenario (for which the p spectral function 
was shown in Fig. 1.21) as solid line. The predictions are in agreement with the 
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data within the systematic uncertainties for m ee > 0.5 MeV/c 2 , but clearly 
do not provide enough yield in the mass range 150 < m ee < 0.5 MeV/c 2 
to describe the data. All three scenarios have in common a sharp drop of 
the p spectral function below m ee < 2m n , something that is not observed in 
data. In data the enhancement stays at an almost constant level down to 
m ee « 100 MeV/c 2 . 

Also the hydrodynamic calculation by K. Dusling is able to describe the 
NA60 data [157] with dilepton rates following from a chiral virial expan- 
sion [158], does not describe the low mass enhancement below pa 600 MeV/c 2 , 
as shown with the red solid line in Fig. 4.35. 

The third theoretical calculation by E. L. Bratkovskaya [156] which is based 
on a microscopic HSD (Hadron String Dynamics) transport model [159] pro- 
vides two scenarios, one is a collisional broadening of the p meson, the other 
includes in addition a drop of the p mass. They are shown in Fig. 4.35 as solid 
green and dashed green line, respectively. As the other model calculations, 
also this does not account for the enhancement observed below 600 MeV/c 2 . 
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Figure 4.35: Invariant mass distribution of e + e~ pairs in min. bias Au + Au 
collisions compared to theoretical predictions including dropping mass and broad- 
ening scenarios by Rapp and Hees [54], Dusling and Zahed [155], Bratkovskaya and 
Cassing [156]. 
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Chapter 5 

Summary and Outlook 



This thesis has presented the analysis of the dielectron continuum in p + p 
collisions at y^s = 200 GeV which allowed to extract a variety of physics 
signals: The measurement of the total charm and bottom cross section in 
the intermediate and high mass region (cr C c = 544 ± 39(stat.) ± 142(syst.) ± 
200(model) fib and a b i = 3.9 ± 2.4(stat.)+2( s Y s t-) A 40 ), the differential cross 
section measurement of u and <fi mesons and the px dependence of the low 
mass dielectron continuum. A detailed study of the spectral shape in the 
mass region 100 < m ee < 300 MeV/c 2 resulted in the measurement of direct 
photons via internal conversions in the pt range of 1 < px < 5 GeV/c which 
is consistent to NLO pQCD prediction. 

These results provide an important baseline for any measurement of the 
dielectron continuum in heavy ion collisions. As part of this thesis a discussion 
of the dielectron continuum measured in Au + Au collisions has been given. 
The results can be summarized briefly as follows: a dielectron yield is ob- 
served in the mass range 150 < m ee < 750 which is enhanced by a factor of 
4.0 ± 0.3(stat.) ± 1.5(syst.) ± 0.8(model) in min. bias Au + Au collisions with 
respect to the expected contribution from hadron decays and semi-leptonic 
charm decays. This enhancement is localized in central collisions and is absent 
in p + p and peripheral collisions. It increases faster than linear with N paxt and 
is moreover localized at low px- The px dependence of the enhancement has 
an inverse slope of T e g = 88 ± ll(stat.)^ 3 (syst.) MeV and does not show any 
signs of radial flow. The concentration of the enhancement at low px is quali- 
tatively consistent with results from SPS experiments [70, 78], but the missing 
evidence for radial flow and the very small effective temperatures seem to indi- 
cate some differences to the observations at SPS. This becomes more evident 
when comparing the enhancement with current models of in-medium modifica- 
tions to the p meson which are able to explain the SPS results with collisional 
broadening of the p [80, 156, 157], but fail to give a complete description of the 
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enhancement observed at RHIC. The measurement of direct photons via inter- 
nal conversions at high pt shows an excess compared to the p+p measurement. 
This excess has an inverse slope of T e ff = 221 ± 23(stat.) ± 18(syst.) MeV in 
0-20% most central collisions and can be described by hydrodynamical models 
with an initial temperature in the range of 300 < T init < 600 MeV at formation 
times of 0.12 < tq < 0.6 fm/c as thermal radiation. 

With upcoming data taking periods and the expected increases in luminos- 
ity, one of the most interesting regions, besides a more quantitative analysis of 
the low mass enhancement, will be the intermediate mass region. Currently, 
the statistics available in Au + Au does not allow to distinguish thermal from 
open charm contributions, but with more statistics a study of the pt depen- 
dence of the IMR may provide new opportunities. This is equally true for 
the study of p + p collisions. The comparison of the measured pr spectrum 
of the IMR and leading-order Pythia calculations suggests that higher order 
processes will need to be considered at high px and invariant masses. With 
more statistics to come the high mass region beyond the J/ip will allow more 
precise measurements of the bottom cross section. The Silicon Vertex Detec- 
tor (VTX) upgrade [160] will accurately measure the collision vertex as well 
as secondary decay vertices. This will tremendously improve the open heavy 
flavour measurement, as it will allow to differentiate between prompt decays 
and the off-vertex decays of charmed and bottom meson. 

Returning to the low mass region, it may be possible (at least in p + p) to 
measure low momentum r] via their Dalitz decays by combining e + e~ pairs with 
100 < m ee < 500 with photons measured in the calorimeter. This would be 
very analog to the beam pipe conversion analysis presented in Appendix B.2. 
A low momentum r\ measurement would significantly reduce the uncertainties 
on the hadronic cocktail and may, together with a good understanding of the 
open charm background, lead the way to a measurement of the u form factor 
in p + p collisions at y/s^N = 200 GeV. 

The Hadron Blind Detector (HBD) upgrade [97] should significantly im- 
prove the signal-to-background ratio of any low mass dielectron analysis, by 
identifying e + e~ pairs with small opening angle as originating from 7r° Dalitz 
decays and photon conversion in a magnetic field free region surrounding the 
beam pipe. The HBD is a windowless proximity-focusing Cherenkov detector 
with a radiator length of 50 cm operating in pure CF 4 as radiator and de- 
tector gas. Cherenkov photons emitted by electrons create photoelectrons in 
a Csl photocathode, which are avalanched by a triple stack of Gas Electron 
Multipliers (GEM) and read out by PCB boards. The operation in a reverse 
bias mode together with a shielding mesh above the photocathode makes the 
HBD blind to minimum ionizing particles, i. e. hadron blind, as their charge is 
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collected on the mesh. e + e~ pairs from tt° Dalitz decays or photon conversions 
with small opening angle will produce overlapping avalanches which could be 
identified by their signal height compared to avalanches from single electrons. 
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Appendix A 

Background Normalization 



A.l Pairing of electrons and positrons 

In the following it is assumed that positrons and electrons are always pro- 
duced in pairs. Let N pairs be produced in a particular event and N is given 
by a probability distribution P(N). Of the N pairs only a fraction p is recon- 
structed, and then the number of reconstructed pairs n p is given by a binomial 
distribution B sampling out of N "events" with a probability e p . 

• Probability to get n p pairs from N true pairs: Lo(n p ) = B(n p , N,e p ) 

• with an average: (n p ) = e p N 

• and variance: a p = e p N(l — e p ) 

Of the remaining pairs one track is reconstructed with a probability e + 
or e_. For a given N and n p the number of additional single positive tracks 
n + and negative tracks n_ follow a multinomial distribution M with possible 
three possible outcomes for each of the N — n p unreconstructed pairs: no track, 
one + track or one — track. 

The probability to get n + and n_ single tracks from N true pairs with n p 
reconstructed pairs, i. e., from N' = (N — n p ) not fully reconstructed pairs is: 
uj(n + , n_) = M(n + , N', e + , e_) 

• with average: (n±) = e±N' 

• variance: a\ = e±N'(l — e±) 

• and covariance: cov(n + ,n„) = — N'£ + £:_ 
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In this case the number of unlike-sign pairs for a given A and n p is: 

N—rip N— n p 

(n+-) = n 2 + n p ^n + CL>(n + ) + n p ^n_cu(n_) + n + n_uj(n + , n_) 

+ + 
= n 2 + n p e + (N - n p ) + n p e^(N - n p ) + (n + n_) 

= n 2 + n p e + (N - n p ) + n p E-(N - + £ + £_(A - n p ) 2 

— e + e^(N — n p ) 
= n 2 + e + Nn p — E + n 2 + e_Nn p — E_n 2 + e + e_N 2 — 2e + e^Nn p 

+ e + E-n 2 — E+E-N + E + E-U p 
= (rip + e + {N - n p )) (rip + £_ (A - n p )) - e + e_{N - n p ). (A.l) 

Similarly one can calculate the number of like-sign pairs: 

2{n ++ ) = ^2(n p + n + )(n p + n + - l)w(n+) 
+ 

= n 2 p -n p + (n\) - (n+) + 2n p (n + ) 

= n 2 p -n p + e 2 + (N - n p ) 2 + e + (l - e + )(N - n p ) - e + (N - n p ) 

+ 2e + u p (N - n p ) 
= n 2 p -n p + e 2 + (N - n p ) 2 - e 2 + (N - n p ) + 2e+n p (N - n p ) (A.2) 

and 

2(n__) = n 2 p - n p + e 2 _(N - n p ) 2 - £ 2 _(N - n p ) + 2e_n p (A - n p ). (A.3) 

To obtain the expected number of like- and unlike-sign pairs for a fixed 
number of real pairs A it is averaged over all possible reconstructed pairs n p : 

n p 

= (1 - e+ - e_ + £+£-)(n 2 p ) + (e+A + e_A - 2e + £-N + E + £_)(n p ) 

+ £ + £^N 2 - £ + £_N 
= (1 - S+ - £_ + £ + £_)(£^ 2 + £ p (l - £ p )A) 

+ (e+A + £_A - 2£ + £_A + e + e_)e p N + e+e.A 2 - e + e_A 

= (£p & p E+ ^p^— £p^+^~ £p£+ £p£— 2E p E+E_ -\- 

• (A 2 - A) + E P N 
= (e p + e+(l - e p )){e p + £_(!- £ P ))(A 2 - A) + e p N. (A.4) 
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Now the like-sign background is calculated: 

2(N ++ ) = J2 2 ( n ++) B M 

= e 2 p N 2 + e p (l - s p )N - e p N + e 2 p e 2 + N 2 + e\e p {\ - s p )N - 2e 2 + e p N 2 
+ e 2 + N 2 - e 2 + N + e 2 + e p N + 2e + e p N 2 - 2e + e 2 p N 2 - 2e + e p (l - e p )N 
= e 2 JN 2 -N) + e 2 + e 2 JN 2 - N) + e 2 + e p N - 2e 2 + e p N 2 + e 2 + (N 2 - N) 



+ f 2 +p 2 p 2 V/V 2 - N\ -9f 2 fJN 2 - ATI +9p 

V 

F , F 

V 

2/ at2 



(e 2 + e 2 + + e 2 + e 2 )(N 2 - N) - 2e 2 + e p (N 2 - N) + 2e + e p (N 2 - N) 
1 



-2e + e 2 (N 2 -N) 



(N++) = - (e p + £+(1 - e P )) 2 (iV 2 - TV) (A.5) 
and 

(N-) = \{s p + - £ p )) 2 (iV 2 - iV). (A.6) 
Finally, it is averaged over all A" to get the foreground unlike-sign pairs: 

N 

= (e p + e + (l - e p ))(e p + e_(l - £ P ))((iV 2 ) - (iV» + e P (iV> 

= ( J BG + _) + (5). (A.7) 

The unlike-sign foreground consists of the sum of the unlike-sign back- 

ground BG + _ and the signal S = e p {N). Similarly the like-sign foreground is 
calculated as: 



(FG ++ ) = J2(N ++ )P(N) 



N 

= l(e P + e + (l-e P )) 2 ((N 2 )-(N}) 

= (BG ++ ) (A.8) 

and 

(FG__) = (BG-). (A.9) 
The like- sign foreground contains no signal. 
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So due to the fact that electrons and positrons are always created in pairs, 
the unlike-sign background is the geometric mean of the like-sign backgrounds, 
independent of the primary multiplicity distribution. 



(BG + -) = 2y/(BG ++ )(BG—) (A.10) 

Comparing the background to the product of the average track multiplici- 
ties one gets for a fixed n p . 



+ 



n p + (n+) 

n p + e + (N -n p ) (A.ir 



averaged over all possible n p : 

( N +) = X>+Mn*) 



= e p N + e + N - e + e p N 

= {e p + e + {l-e p ))N (A.12) 
or averaged over all possible N: 

(FG + ) = ^2(N + )P(N) 

N 

= (s p + e + (l-e p ))(N) (A.13) 

and thus: 

(FG + )(FG + ) = (e p + e + (l - e p ))(e p + e_(l - e p ))(N) 2 (A.14) 



or 



<« G +-> - l + ^M. (A.15) 



(TO + )(FG_) (NY 

So in general (BG + ^) ^ (FG + )(FG-), except for the special case that 
P(N) is a Poisson distribution. Note this is the opposite conclusion one derives 
in the case that the sources of + and — tracks are independent, i.e., + and 
— tracks are produced as singles and not as pairs like they are for muons. In 
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that case (FG + ){FG-) is the correct background normalization. 
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Appendix B 



Beam Pipe Conversions 



B.l Introduction 

In this Appendix another approach to measure direct photons at low is 
presented. It is the measurement of real photons which convert externally in 
the beam pipe, surrounding the collision vertex at a radius of 4 cm. These 
are the pairs which are removed with the <pv cut in the regular dielectron 
continuum and internal conversion analysis. The basic idea is to identify e + e~ 
pairs from beam pipe conversions and to combine these with photons in the 
EMCal to tag conversion photons from ir° decays. It is certainly not as strong 
as the virtual photon analysis but is nevertheless an interesting approach. 



B.2 Analysis 

To identify e + e~ pairs from photon conversions, a single electron identification 
cut is applied, which requires signals from at least two phototubes in the Ring 
Imaging Cherenkov Detector (RICH) matching to a reconstructed charged 
track in the Drift Chamber (DC). No further electron identification cuts were 
applied since the pair cuts (see Section B.2.1) to separate conversion photons 
from other e + e~ pairs are more efficient and powerful enough to provide a very 
clean photon conversion sample. 

The extracted photon conversions are tagged with photons reconstructed 
in the EMC to determine the contribution from ir° — > 77 decays (see Sec- 
tion B.2.2). 

All yields are measured as a function of px of the e + e~-pair, which makes 
a direct comparison of the inclusive photon yield, iV" lcl , and the tagged photon 
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yield, iVT tag , 



possible: 



AT 1 M = e e 
N f^ (pr) = e e 



+e - a e+e - 7 incl (pr) 



e 7 (pr) f l 77 (pr) 



(B.l) 
(B.2) 



The measured yield of inclusive photons depends on the reconstruction effi- 
ciency e e + e - and the PHENIX acceptance a e + e - of the conversion e + e~ pair. 
The tagged photon yield depends in addition on the efficiency to reconstruct 
the second photon in the EMC e 7 (pr) and on the conditional probability / to 
find it in the EMC acceptance, given that the e + e~ pair has been reconstructed 
already Here, c 7 (pt) is weighted with the pr distribution of the e + e~ pair. In 
the ratio N™ d /N* ta9 the e+e" pair reconstruction efficiency and acceptance 
correction factor cancel. 

A ratio of the hadronic decay photon yield, A^ adr , and the tagged photon 
yield from 7r° decays, A^° tag , is calculated with simulations. 



The comparison of the ratio in data and in simulations in a double ratio 
leads to an expression that is equivalent to the ratio of inclusive and decay 
photons as shown in Eq. (B.4). 



The only remaining factors are the reconstruction efficiency of the photon in 
the EMC, 6 7 (pt), and the conditional acceptance / in the simulation part of the 
double ratio, which have both been determined with Monte Carlo simulations 
(see Section B.2. 3). 

B.2.1 Photon Conversions 

Since the PHENIX tracking algorithm assumes the track to originate from the 
collision vertex, off-vertex conversion pairs are reconstructed with an artificial 
opening angle which leads to an invariant mass that is proportional to the 
radius at which the conversion occurs. 

Therefore, photon conversions that occur in the beam pipe material (Be, 
0.3 % radiation length) at a radius of 4 cm are reconstructed with an invariant 
mass of ~ 20 MeV/c 2 . Fig. B.l shows an invariant mass spectrum of e + e~ pairs 
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(Pr) = f NZ fa) 



(B.3) 




(B.4) 
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Figure B.l: Invariant mass of e + e pairs before (black) and after (red) applying 
cuts on the orientation of the e + e~ pair in the magnetic field. 



in the range 0-0.1 GeV/ c 2 . The peak from photon conversions in the beam pipe 
at 20 MeV/c 2 can be clearly separated from Dalitz decays 7r° — > 7e + e~, which 
dominate the spectrum below 10 MeV/c 2 , and combinatorial background pairs, 
whose contribution increases toward higher invariant masses. 

The photon conversion pairs, which have no intrinsic opening angle, can 
be distinguished from Dalitz decays and purely combinatorial pairs by cutting 
on the orientation of the e + e~ pair in the magnetic field. 

Fig. B.l shows the invariant mass spectra of e + e~ pairs before (black) 
and after (red) applying these pair cuts. The yield from integrating the mass 
region < 35 MeV/c 2 of the conversion peak is corrected for the remaining 
dependent contamination of ~ 15.0 ± 2.0 (syst) % due to combinatorial e + e~ 
pairs which has been determined with mixed events. 

B.2.2 Tagging of Decay Photons 

To reveal which of these conversion photons come from 7r° — > 77 decays, the 
e + e" pairs in the conversion peak are combined with photons which have been 
measured in the EMC, under loose cuts based on the time of flight and the 
shower profile for photons with a minimum px of 0.3 GeV/c, and their invariant 
mass is calculated (see Fig. B.2). 
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ye + e" inv. mass with: 0.8 GeV < e + e" pair-p T < 5.0 GeV 
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Figure B.2: Invariant mass of "fe + e~ triplets in same events (black) and normalized 
mixed events (red) for e + e~ pairs with 0.8 < pr < 5.0 GeV/c. The insert shows 
the invariant mass of "fe + e~ triplets after background subtraction for e + e~ pairs 
with 0.8 < px < 1.2 GeV/c. A fit with a Gaussian is drawn and the resulting 
parameters shown in the box in the upper right of the graph. 



The reconstruction efficiency e 7 (pr) of the loose photon has been estimated 
with a full GEANT simulation which embed simulated photons into real EMC 
data, therefore providing a combined information on the photon identification 
efficiency and occupancy effects. The overall efficiency is determined to be 
82 ± 1 % independent of pt beyond the minimum pt cut off. 

Conversion photons that are identified as decay products of 7r° can be 
tagged as tag . This signal has a large combinatorial background due to the 
high photon multiplicity in Au + Au collisions. 

The combinatorial background is reproduced with an event mixing method, 
which creates uncorrelated pairs of photons and e + e~ pairs from different 
events. The mixed event spectrum is normalized to the same event spectrum 
well outside the vr° mass region (0-100 MeV/c 2 , 170-250 MeV/c 2 ) and sub- 
tracted. 

The statistical error on the normalization factor is on the order of 0.2 % and 
depends only on the statistics in the same event spectrum in the normalization 
region. As an example, the resulting 7r° signal for e + e~ pairs with 0.8 < pr < 
1.2 GeV/c is shown as insert in Fig. B.2. 
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Mean and a are determined by a fit of the background subtracted data with 
a Gaussian. The data are also fitted to the sum of a second order polynomial 
and a Gaussian, to take into account the possibility that the shape is not 
completely described by the mixed event spectrum. The difference in the 
resulting mean and a is negligible. The mean and a obtained by the fit are 
then used to integrate the data in a region ±1.5 a around the mean, chosen 
to optimize the signal to background ratio. 

The statistical error on the extracted 7r° signal is given by: 

°l = E FG ^ + a E BG '^ + f 77 E BG '^) (B-5) 

i i \ i J 

With FG(i) and BG'{i) being the yields in bin i of invariant mass spectrum 
in same events and normalized mixed events, respectively, the summations are 
performed over the integration region. It is important to note that the last 
term in Eq. (B.5), is the square of the sum over the normalized background, 
and therefore, depends on the integration region and is not bin independent. 
Different integration regions have been used. Variations in the resulting yield 
have been used to set a systematic uncertainty on the yield extraction of 2.5 % 
independent of px- 

The loss of i\^"° tag due to the external conversion of the second photon 
is corrected by a factor 1 — p CO nv = 94 ± 2 %. In this factor p conv is the 
conversion probability due to the material budget between the vertex and the 
Pad Chamber 3 (PC3) in front of the EMC. 

B.2.3 Simulations 

The contribution of hadronic decays has been determined with a fast Monte 
Carlo simulation of ir° and 77 Dalitz decays. A parameterization of the n° spec- 
trum measured by PHENIX [46] has been used as input. The 77 distribution 

has been generated assuming my scaling (p? —> \Jpt^~ m ^ ~ m ^> ) of the 7i~° 
spectral shape and a normalization at highpr to 77/71"° = 0.45 ±0.04, according 
to PHENIX data [13]. The relative error of 9 % on the 77/71"° ratio is reduced 
by the branching ratio of the two photon decay and results in a 3 % error in 
the ratio N^ dl /Nf. 

The contamination due to neutral Kaons which decay before the beam pipe 
has been found negligible (~ 1 %) and has been folded into the systematic error 
on the simulations. 

The conditional probability / that the photon is reconstructed in the EMC 
once the e + e~ pair is reconstructed already was calculated with a fast Monte 
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Figure B.3: Direct photon excess, (a) The fraction of the direct photon component 
as a function of pj- in p+p. (b) Au+Au (min. bias) compared to other measurements 
of direct photons as described in the legend. The curves are from a NLO pQCD 
calculation [41, 138]. 

Carlo simulation of tt° — > r ye + e~ . The use of Dalitz decays is justified by the 
fact that the pt spectra of photons from n — > 77 are essentially identical to 
the e + e~ pair px spectrum from 7r° — > r ye + e~. 

After the e + e~ pair has been filtered in the detector acceptance, the con- 
ditional acceptance / for the second photon has been calculated taking dead 
areas of the detector into account. Uncertainties in calculating / are found to 
be 5 %, which is the largest source of systematic errors. The limited energy 
resolution of the EMC of <Je/E = 5% © 9%/^/E introduced an additional 
systematic error of ~ 1 % due to the pt cut at 0.3 GeV/c. 



B.3 Conclusions 

The right panel of Fig. B.3 shows the preliminary result for the double ratio 
7 mcl (p T ) /7 hadr (p T ) as in Eq. (B.4) for Minimum Bias Au + Au collisions at 
y/SNN — 200 GeV. The main sources of systematic errors arise from the uncer- 
tainties in the description of the detector active areas, in the peak extraction 
and in the assumptions of the tt° shape and give a final systematic error on the 
double ratio of ~ 7%. The result is compared to the conventional statistical 
subtraction of hadronic decay photons [161], and a result [49] which is based 
on the same tagging method, but instead of photons coming from conversions 
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in the beam pipe, the clean photon sample is determined by selecting EMC 
clusters with very strict photon identification cuts. 

In addition the result of the internal conversion analysis is shown as 1 +r rs 
7inci/7hadr in the right panel for Au + Au and in the left of Fig. B.3 for p + p 
collisions. While all results agree within their uncertainties, the improvement 
in statistical and systematic uncertainties with the internal conversion analysis 
is quite significant. Also the NLO pQCD predictions of W. Vogelsang [41, 138] 
are shown.. 
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Appendix C 
Data Tables 



Table C.l: ui and 4> yields. 



p T (GeV/c) (0.740-0.815 GeV/c 2 ) iV (0.965-1.065 GeV/c 2 ) 



0.10 


33.6 


± 


9.0 


30.4 


± 


9.5 


0.30 


153 


± 


19 


141 


± 


18 


0.50 


372 


± 


28 


185 


± 


19 


0.70 


508 


± 


29 


223 


± 


20 


0.90 


358 


± 


24 


193 


± 


18 


1.10 


261 


± 


20 


140 


± 


15 


1.30 


182 


± 


16 


93 


± 


12 


1.50 


134 


± 


13 


57.6 


± 


9.2 


1.70 


68.3 


± 


9.3 


40.5 


± 


7.4 


1.90 


65.3 


± 


8.9 


28.4 


± 


6.2 


2.50 


142 


± 


13 


91 


± 


10 


4.00 


51.4 


± 


7.4 


31.2 


± 


5.8 


5.50 


3.8 


± 


2.0 


0.9 


± 


1.0 
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Table C.2: Fit results of uj yield extraction. 



Yield 

Pt (GeV/c) Gauss Gauss + PolO Gauss + Poll Gauss + Pol2 



o.io 


34.3 


1 

± 


9.7 


4/ 


1 


"1 o 

12 


1 1 

41 


1 

± 


12 


45 


1 

± 


"1 f 

lb 


0.30 


145 


± 


21 


126 


± 


23 


125 


± 


23 


116 


± 


26 


0.50 


247 


± 


26 


212 


± 


29 


210 


± 


29 


186 


± 


38 


0.70 


484 


± 


37 


359 


± 


33 


363 


± 


33 


299 


± 


34 


0.90 


487 


± 


35 


330 


± 


32 


340 


± 


33 


289 


± 


35 


1.10 


301 


± 


23 


230 


± 


22 


234 


± 


22 


215 


± 


21 


1.30 


215 


± 


19 


164 


± 


19 


166 


± 


20 


141 


± 


20 


1.50 


171 


± 


20 


10 


± 


14 


104 


± 


15 


85 


± 


15 


1.70 


140 


± 


16 


83 


± 


15 


85 


± 


15 


78 


± 


19 


1.90 


59.6 


± 


8.8 


51.2 


± 


9.0 


49.6 


± 


9.1 


43.4 


± 


9.2 


2.50 


178 


± 


15 


137 


± 


15 


139 


± 


15 


117 


± 


16 


4.00 


47.9 


± 


7.4 


50.6 


± 


7.7 


43.0 


± 


7.6 


52 


± 


10 



Sigma (MeV/c 2 ) 

Pt (GeV/c) Gauss Gauss + PolO Gauss + Poll Gauss + Pol2 



0.10 


17.6 


± 


5.4 


22.5 


± 


6.2 


20.2 


± 


6.0 


21.6 


± 


6.9 


0.30 


8.2 


± 


1.6 


7.2 


± 


1.6 


7.1 


± 


1.6 


6.7 


± 


1.7 


0.50 


16.7 


± 


2.4 


14.3 


± 


2.7 


14.1 


± 


2.7 


12.2 


± 


3.6 


0.70 


10.2 


± 


1.5 


6.89 


± 


0.78 


7.00 


± 


0.81 


5.92 


± 


0.71 


0.90 


25.7 


± 


3.3 


16.1 


± 


1.8 


16.6 


± 


1.9 


14.5 


± 


1.9 


1.10 


15.6 


± 


1.7 


11.7 


± 


1.3 


11.9 


± 


1.3 


11.1 


± 


1.3 


1.30 


20.6 


± 


3.0 


14.7 


± 


2.1 


15.0 


± 


2.2 


13.0 


± 


2.0 


1.50 


35.1 


± 


9.4 


11.4 


± 


2.3 


12.0 


± 


2.4 


9.3 


± 


2.5 


1.70 


43.5 


± 


9.7 


21.2 


± 


4.2 


21.8 


± 


4.5 


20.4 


± 


4.9 


1.90 


14.4 


± 


2.6 


12.4 


± 


2.3 


12.1 


± 


2.2 


10.9 


± 


2.3 


2.50 


20.6 


± 


2.9 


14.7 


± 


1.7 


15.0 


± 


1.7 


13.1 


± 


1.6 


4.00 


21.2 


± 


2.2 


22.8 


± 


2.5 


20.5 


± 


2.9 


24.7 


± 


4.2 



Continued on next page 
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Table C.2: (continued) 
X 2 /ndf 



p T (UeV/cJ 


Gauss 


Lrauss + roll) 


(~ ^ l T ^ , > 1 1 

LraUSS + rOll 


Liauss + rOlZ 


0.10 


25.51/28 


19.02/27 


16.98/26 


16.72/25 


0.30 


34.32/28 


32.22/27 


31.84/26 


31.47/25 


0.50 


52.24/28 


46.99/27 


46.30/26 


45.18/25 


0.70 


87.42/28 


48.89/27 


44.51/26 


32.96/25 


0.90 


86.84/28 


37.16/27 


31.66/26 


26.28/25 


1.10 


94.51/28 


38.00/27 


31.98/26 


28.46/25 


1.30 


67.84/28 


42.57/27 


40.62/26 


34.94/25 


1.50 


90.59/28 


51.59/27 


45.66/26 


38.84/25 


1.70 


43.87/28 


26.07/27 


25.52/26 


25.22/25 


1.90 


39.52/28 


31.76/27 


30.16/26 


26.10/25 


2.50 


89.68/28 


42.62/27 


37.99/26 


27.94/25 


4.00 


47.46/28 


44.28/27 


26.22/26 


18.73/25 
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Table C.3: Fit results of 4> yield extraction. 



Yield 

Pt (GeV/c) Gauss Gauss + PolO Gauss + Poll Gauss + Pol2 



(J. 10 


on i 

29.1 


1 


7.9 


on A 


1 

± 


o o 

8.8 


29.8 


1 


8.7 


27.7 


1 


8.7 


0.30 


139 


± 


16 


120 


± 


16 


119 


± 


16 


124 


± 


17 


0.50 


159 


± 


17 


134 


± 


17 


132 


± 


18 


130 


± 


18 


0.70 


169 


± 


18 


150 


± 


18 


139 


± 


18 


134 


± 


18 


0.90 


206 


± 


18 


181 


± 


19 


168 


± 


19 


188 


± 


21 


1.10 


151 


± 


16 


139 


± 


18 


138 


± 


18 


142 


± 


21 


1.30 


117 


± 


15 


103 


± 


16 


101 


± 


17 


97 


± 


19 


1.50 


67 


± 


10 


71 


± 


12 


75 


± 


13 


100 


± 


31 


1.70 


25.4 


± 


5.7 


30.4 


± 


5.9 


28.2 


± 


6.0 


29 


± 


6.0 


1.90 


28.4 


± 


6.0 


35.1 


± 


6.5 


27.3 


± 


6.1 


27.3 


± 


6.2 


2.50 


86.5 


± 


9.9 


96 


± 


10 


83 


± 


10 


84 


± 


10 


4.00 


20.4 


± 


4.9 


25.8 


± 


5.1 


23.2 


± 


5.5 


21.4 


± 


5.6 



Sigma (MeV/c 2 ) 

Pt (GeV/c) Gauss Gauss + PolO Gauss + Poll Gauss + Pol2 



0.10 


9.1 


± 


3.8 


10.4 


± 


4.2 


9.4 


± 


4.1 


8.6 


± 


5.2 


0.30 


5.48 


± 


0.72 


4.79 


± 


0.68 


4.77 


± 


0.68 


4.92 


± 


0.71 


0.50 


13.1 


± 


1.3 


11.8 


± 


1.3 


11.7 


± 


1.3 


11.6 


± 


1.3 


0.70 


6.55 


± 


0.76 


5.97 


± 


0.71 


5.66 


± 


0.68 


5.54 


± 


0.68 


0.90 


15.9 


± 


1.4 


14.4 


± 


1.4 


13.7 


± 


1.4 


14.8 


± 


1.6 


1.10 


19.4 


± 


2.4 


18.0 


± 


2.4 


18.0 


± 


2.5 


18.4 


± 


2.7 


1.30 


23.7 


± 


3.6 


21.2 


± 


3.5 


20.8 


± 


3.6 


20.3 


± 


3.8 


1.50 


18.4 


± 


4.0 


19.5 


± 


4.7 


21.0 


± 


5.2 


29.5 


± 


9.4 


1.70 


5.5 


± 


1.9 


7.0 


± 


1.4 


6.4 


± 


1.7 


6.5 


± 


1.6 


1.90 


8.7 


± 


2.6 


16.5 


± 


7.5 


8.3 


± 


2.5 


8.3 


± 


2.6 


2.50 


14.3 


± 


2.0 


16.9 


± 


3.2 


13.6 


± 


1.9 


13.8 


± 


1.9 


4.00 


21.2 


± 


2.6 


26.8 


± 


2.4 


24.9 


± 


3.1 


25.0 


± 


3.7 



Continued on next page 
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Table C.3: (continued) 
X 2 /ndf 



p T {(jtev /c) 


Gauss 


Lrauss + rOlU 




LrauSS + rOlz 


0.10 


18.57/28 


17.03/27 


15.78/26 


15.13/25 


0.30 


49.46/28 


33.43/27 


33.37/26 


32.19/25 


0.50 


51.69/28 


33.32/27 


33.23/26 


32.99/25 


0.70 


27.39/28 


17.95/27 


12.29/26 


11.45/25 


0.90 


60.09/28 


47.10/27 


41.32/26 


31.12/25 


1.10 


32.43/28 


29.74/27 


29.73/26 


29.56/25 


1.30 


19.12/28 


16.33/27 


16.06/26 


15.94/25 


1.50 


52.02/28 


51.64/27 


48.58/26 


46.97/25 


1.70 


133.5/28 


65.33/27 


61.16/26 


60.53/25 


1.90 


80.72/28 


51.22/27 


36.73/26 


36.72/25 


2.50 


111.4/28 


79.21/27 


60.80/26 


50.34/25 


4.00 


57.31/28 


24.48/27 


23.07/26 


20.02/25 
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Table C.4: Invariant Cross Section of u. 



Pt (GeV/c) Ed 3 a/dp 3 (mbarn GeV 2 /c 3 ) stat. error syst. error 



0.10 


2.91 


7.7X10- 1 


7.4X10- 1 


0.30 


2.12 


2.7X10- 1 


5.4X10" 1 


0.50 


1.39 


l.OxlO" 1 


3.5X10" 1 


0.70 


8.56x10^ 


5.0xl0- 2 


2.2X10" 1 


0.90 


3.67X10" 1 


2.4xl0^ 2 


9.3xl0~ 2 


1.12 


1.49X10- 1 


l.lxlO^ 2 


3.8xl0~ 2 


1.38 


1.049X10" 1 


7.6xl0~ 3 


2.6xl0~ 2 


1.75 


2.20xl0- 2 


1.8xl0~ 3 


5.5xl0- 3 


2.50 


2.15xl0- 3 


2.0xl0~ 4 


5.4xl0~ 4 


3.50 


2.50xl0~ 4 


4.0xl0~ 5 


6.3xl0~ 5 


4.50 


3.1xl0~ 5 


l.lxl0~ 5 


7.7xl0- 6 


6.50 


3.7xl0- 7 


2.0xl0~ 7 


2.0xl0~ 7 



Table C.5: Invariant Cross Section of (j). 



p T (GeV/c) 


Ed 3 a/dp 3 (mbarn GeV 2 /c 3 ) 


stat. error 


syst. error 


0.10 


1.14X10" 1 


3.5xl0~ 2 


3.4xl0~ 2 


0.30 


1.79X10- 1 


2.2xl0~ 2 


5.2xl0~ 2 


0.50 


1.17X10- 1 


1.2xl0~ 2 


3.4xl0~ 2 


0.70 


8.47xl0~ 2 


7.6xl0~ 3 


2.5xl0~ 2 


0.90 


5.03xl0~ 2 


4.7xl0~ 3 


1.5xl0~ 2 


1.12 


2.09xl0- 2 


2.3xl0~ 3 


4.9xl0~ 3 


1.38 


1.45xl0~ 2 


1.5xl0~ 3 


3.4xl0~ 3 


1.75 


2.97xl0~ 3 


3.8xl0~ 4 


6.9xl0~ 4 


2.50 


4.90xl0" 4 


5.7xl0~ 5 


l.lxlO" 4 


3.50 


5.1xl0~ 5 


l.lxl0~ 5 


1.2xl0~ 5 


4.50 


6.6xl0~ 6 


2.7xl0~ 6 


1.5xl0~ 6 


6.50 


1.9xl0~ 8 


3.9xl0- 8 


4.0xl0~ 8 
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